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ABSTRACT 


This  program  was  designed  to  permit  the  ONR  sponsored  program  in  the  Materials  Research 
Laboratory  at  Penn  State,  now  a  component  in  the  Materials  Research  Institute,  to  draw  upon 
the  expertise  of  Professor  Jan  Fousek  an  internationally  respected  authority  on  domain 
processes  in  ferroelectrics.  The  program  has  sponsored  Dr,  Fousek  for  3  month  visits  to  Penn 
State  in  1999, 2000  and  2001,  Dr.  Fousek  is  deeply  respected  in  the  MRL  and  has  been  able  to 
interact  effectively  and  provide  advice  which  is  bearing  fruit  in  both  practical  and  fundamental 
studies  now  ongoing  and  also  contemplated  in  future  studies.  Topics  which  have  benefited 
most  from  the  interaction  are  briefly  discussed  in  the  report  and  several  publications  which 
were  greatly  facilitated  are  included  as  appendices  to  this  report. 

1.  INTRODUCTION 

Domain  wall  processes  are  widely  accepted  to  be  of  major  importance  in  facilitating  the 
performance  of  piezoelectric  materials  derived  from  oxygen  octahedron  type  ferroelectrics 
with  centric  prototypic  symmetrys  as  exemplified  in  the  most  widely  used  perovskites  structure 
material  families.  The  solid  solution  system  between  lead  zirconate  and  lead  titanate  (PZT)  is 
much  the  most  widely  used  piezoelectric  ceramic,  and  here  the  compositions  of  most  interest 
cluster  around  the  morphotropic  phase  boundary  (MPB)  between  ferroelectric  rhombohedral 
and  ferroelectric  tetragonal  forms.  Unfortunately  at  the  52%  Zr/48%  Ti  composition  of  this 
phase  boundary  PZT  has  proven  impossible  to  grow  in  single  crystal  form  however,  recently 
lead  zinc  niobate:  lead  titanate  and  lead  magnesium  niobate:  lead  titanate  have  been  shown  to 
have  pseudo-morphotropie  rhombohedral:  tetragonal  ferroelectric  phase  boundaries  and  both 
can  be  grown  (with  some  difficulty)  in  effective  single  crystal  form.  Now  with  control  of 
orientation  of  the  field  with  respect  to  the  crystal  axes,  the  concepts  of  Domain  Engineering 
and  Domain  Averaging  could  be  examined  with  Fousek,  the  huge  anisotropy  in  properties 
explained  and  massive  anhysteretic  strain  capability  rationalized. 

A  second  area  which  has  been  most  productive  was  the  suggestion  to  explore 
flexoelectricity  i.e.  the  generation  of  electric  field/polarization  by  an  elastic  strain  gradient. 

The  work  stemmed  from  discussions  with  Fousek  concerning  composites  with  oriented  cone 
shaped  inclusions  which  from  the  shape  must  be  piezoelectric.  We  concluded  the  origin  of  the 


charge  separation  must  be  the  strain  gradient  induced  by  the  cone  shape  and  determined  to 
explore  gradient  induced  (flexoelectric)  effects  to  see  if  useful  properties  could  result  which 
would  be  very  different  from  conventional  ceramics. 

Throughout  the  full  period  of  the  grant,  Fousek  and  Cross  have  cooperating  in  assembling 
material  for  a  book  on  Domain  Walls  and  Switching  in  Ferroic  Materials.  Unfortunately,  the 
magnitude  of  the  task  of  writing  a  truly  authoritative  work  on  this  topic  was  originally  grossly 
underestimated  and  the  original  proposed  time  scale  has  now  expanded  from  3  to  5  years.  The 
whole  significant  new  area  of  domain  studies  in  thin  films  proved  a  major  difficulty  and  we 
were  fortunate  to  be  joined  by  Dr.  Alexander  Tagantsev  as  a  co-author.  Now  all  chapters  are 
undergoing  second  revisions  and  the  deadline  for  submission  of  the  manuscript  is  July  2003. 

A  third  area  of  activity  concerned  the  important  task  of  completing  the  new  IEEE  standard 
on  Ferroelectricity  which  involves  a  number  of  MRL/MRI  faculty.  The  proper  definition  of 
polarization  in  a  pyroelectric  but  non-ferroelectric  polar  crystal  is  non  trivial  but  of  significant 
importance  as  ferroelectrics  can  occur  from  polar  prototypic  groups.  Dr.  Fousek  was  of  great 
assistance  in  removing  the  inconsistencies,  which  had  vitiated  the  original  standard. 

2.  APPROACH 

2.1  DOMAIN  RELATED  STUDIES 

The  task  for  this  program  was  not  to  complete  specific  studies,  but  rather  to  catalyze  and 
explore  new  areas  related  to  current  MRL  programs  which  could  have  strong  relevance  to  navy 
needs.  It  was  natural  with  the  major  database  being  generated  on  domains,  domain  walls  and 
switching  for  the  book,  to  be  in  discussion  on  the  fascinating  domain  related  properties  of  the 
important  PZN:PT  and  PMN:PT  single  crystals  being  generated  in  the  laboratory. 

The  publication  “Engineering  Multidomain  Ferroic  Samples”  by  Fousek  and  Cross 
(Appendix  1)  was  probably  the  first  attempt  to  systematize  the  possibilities  of  controlling 
domain  structure  to  provide  regular  spacial  arrays  we  called  domain  geometry  engineered 
samples  where  the  domain  is  tuned  to  correspond  to  the  k  vectors  of  fields  propagating  through 
the  sample  and  domain  average  engineered  samples  where  the  crystal  is  sub  divided  into  a  very 
large  number  of  domains  so  as  to  yield  on  averaged  “texture”  symmetry.  The  concept  is 
further  refined  in  the  letter  “Domain  Geometry  engineering  and  domain  average  engineering  of 
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ferroics.”  (Appendix  2)  by  Fousek,  Litvin  and  Cross.  Here  the  possibilities  for  generating  a 
whole  range  of  averaged  symmetries  for  particular  sub  sets  of  domains  are  discussed. 

The  importance  of  the  domain  average  engineered  structure  in  the  001  field  poled 
rhombohedral  phase  of  the  0.955  lead  zinc  niobate  0.045  lead  titanate  for  high  strain  actuation 
was  underscored  in  the  paper  “Domain  and  Phase  Change  Contributions  to  Response  in  High 
Strain  Piezoelectric  Actuators”  by  Cross.  (Appendix  3)  Here  there  is  no  switching  force  on  the 
domain  walls  in  the  001  poled  condition  so  that  the  induced  strain  is  intrinsic  and  shown  to  be 
associated  with  a  rotation  of  the  1 1 1  polarization  vector  in  the  domain  towards  the  001  field 
direction.  The  calculated  titling  of  Ps  was  later  confirmed  by  optical  studies  by  Z.  Ye  (1). 

A  further  extension  of  the  symmetry  arguments  to  consider  explicitly  the  polar  character 
of  the  domain  is  in  the  paper  “Symmetry  and  antisymmetry  in  Electroceramics”  by  Newnham 
and  Cross  (Appendix  4)  where  the  sample  geometry  is  used  to  control  modes  of  vibration  and 
their  interaction  with  the  polar  structure  is  a  Domain  Geometry  Engineered  context. 

A  fascinating  study,  triggered  by  interaction  between  Fousek  and  Cao  concerned  the 
characterizations  of  the  strange  S  wall  in  Ferroelectries  (Appendix  5).  The  S  walls  refer  to  the 
domain  walls  whose  orientations  do  not  coincide  with  simple  crystallographic  planes  but  are 
planes  of  elastic  strain  compatibility  and  change  with  temperature.  A  simple  continuum  model 

is  constructed  to  describe  S  wall  structures  in  a  proper  ferroelastic  (m3  m—+  mmm  ferroelastic 
phase  transition).  This  model  describes  local  values  of  elastic  strain  components  and  the 
temperature  dependence  of  the  S  wall  orientation.  Variation  of  strain  across  the  S  wall  are  also 
derived. 

2.2  STRAIN  GRADIENT  RELATED  PROPERTIES 

The  major  area  of  activity  on  strain  gradient  effects  in  dielectrics,  flexoelectricity  was 
initiated  through  discussion  with  Dr.  Fousek  where  we  were  looking  for  an  explanation  of  the 
shape  induced  piezoelectricity  in  composites,  discussed  in  Appendix  6.  To  determine  practical 
feasibility  of  such  a  composite  measurement  of  coefficients  puiu  Pi  122,  P1212,  in  lead 
magnesium  niobate:  lead  titanate  ceramic  Appendix  7  were  initiated. 

The  very  high  values  of  the  coefficient  \in  (pi  122)  as  compared  to  earlier  measurements  on 
polymers  lead  us  to  suspect  the  relaxor  character  in  PMN:PT  Appendix  8  and  a  strong  role  of 
micro-polar  regions.  Further  studies  on  excellent  Bao.7Sro,3Ti03  ceramics  however,  revealed 
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constants  even  10  times  higher  than  in  PMN:PT  (Appendix  9).  A  very  interesting  feature  of 
the  (BaSr)Ti03  data  was  a  steep  upturn  in  the  measured  \in  in  the  15°  C  region  above  but  close 
to  Tc  where  C  departs  from  regular  Curie  Weiss  behavior. 

In  further  extension  of  these  studies  it  was  decided  to  explore  an  unpoled  ferroelectric 
ceramic  under  very  large  flexing  stress  to  see  if  it  would  be  possible  to  induce  piezoelectric 
poling  (Appendix  10)  using  PZT  5  H  obvious  ferroelastic  poling  could  be  induced  at  the  higher 
stress  levels,  and  a  clear  softening  of  the  modulus  was  evident  suggesting  the  onset  of 
ferroelastic  wall  motion.  Closely  coupled  to  this  change  was  an  obvious  upturn  in  the 
flexoelectric  coefficient  from  0.5  pc/m  to  2.0  pc/m  strongly  suggesting  an  extrinsic 
contribution  to  pn  triggered  by  the  ferroelastic  wall  motion. 

We  are  excited  that  NUWC  at  Newport  Rhode  Island  have  excellent  facilities  for  very 
high  strain  4  point  bending  tests  and  are  interested  to  cooperate  in  extending  the  flexoelectric 
studies  to  other  ceramic  and  single  crystal  samples. 

A  third  area,  which  was  of  considerable  help,  concerned  the  importance  task  of  updating 
the  IEEE  standard  on  Ferroelectricity,  a  task  taken  on  by  the  MRL  faculty.  Here  the  proper 
definition  of  polarization  in  a  pyroelectric  but  no  ferroelectric  crystal  is  non  trivial,  but  very 
important.  The  original  standard  was  demonstrably  wrong  in  this  context,  but  with  Fousek’s 
help  we  were  able  to  arrive  at  a  physically  consistent  definition. 

3 .  SUMMARY  OF  SCIENCE  AND  TECHNOLOGY  INITIATED 

The  task  for  the  visits  of  Professor  Jan  Fousek  was  not  to  complete  studies  on  this 
program,  but  rather  to  initiate  and  advance  important  new  concepts. 

(i)  Domain  Geometry  and  Domain  Average  Engineering 

Important  MRL/MRI  programs  initiated  in  contact  with  Fousek  are  summarized  in  Appendixes 
1-4.  The  first  initiative  in  the  laboratory  on  the  theory  of  strange  S  walls  in  Appendix  5  opens 
the  more  complex  topic  of  strange  walls  in  ferroelectrics  and  the  rubber  like  elasticity  which 
can  be  generated. 

(ii)  Flexoelectricity:  This  whole  exciting  new  mode  of  charge  separation  under  strain 
gradient  stems  from  Fousek’s  curiosity  regarding  shape  induced  piezoelectricity  in  composites. 
(Appendix  6-9).  Clearly  there  is  a  need  for  more  measurements  in  high  permittivity  high 
symmetry  systems  to  understand  the  origins  of  the  phenomenon  and  to  provide  reliable 
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theoretical  underpinning.  New  measurements  on  unpoled  high  symmetry  systems  suggest  that 
the  DC  flexoelectric  effects  may  give  us  new  modes  for  separating  intrinsic  and  extrinsic 
domain  related  behaviors  (Appendix  10). 

(iii)  The  importance  of  a  proper  Ginsburg:  Landau  formulation  for  antiferroelectrics  is 
highlighted  in  the  works  by  Hatt  and  Gao.  (Appendix  11,12)  Unlike  earlier  Landau 
formulations  this  approach  permits  the  proper  modeling  of  the  group  symmetry  and  of  twin  and 
antiphase  boundaries, 

(iv)  Fousek’s  most  effective  help  is  in  formulating  vital  parts  of  the  new  IEEE  standard 
on  ferroelectrics  took  a  major  load  from  MRL  ONR  related  faculty  and  helped  to  clarify  the 
thinking  of  both  students  and  faculty  in  MRI. 

4.  IMPACT/NAVY  RELEVANCE 

Ferroelectric  ceramic  systems  are  very  widely  used  in  Navy  undersea  sonar  systems.  By 
its  very  nature  the  poled  ferroelectric  piezoelectric  ceramic  has  strong  local  gradients  in 
electric  and  elastic  boundary  conditions.  Flexoelectric  effects  must  be  present  and  in  many  of 
the  best  “soft”  compositions  relaxor  effects  occur.  We  will  be  providing  the  quantitative 
background  to  understand  the  importance  of  these  effects.  In  the  newer  single  crystal 
transducers  domain  averaging  is  exploited  in  the  ultra  high  coupling  rhombohedral 
compositions.  The  complexity  of  possible  averaging  was  ill  appreciated. 

5.  PLANNED  SCIENCE  and  TECHNOLOGY 

It  is  anticipated  that  although  the  formal  program  was  terminated  in  September  2001 
personal  contacts  will  continue  with  exchange  of  ideas  on  both  fundamental  and  practical 
issues  in  ferroelectric  crystals  and  ceramics. 

6.  TECHNOLOGY  TRANSFER 

A  continuous  flux  of  small  and  mid-sized  companies  are  involved  with  material  and 
device  studies  in  MRI/MRL  and  afford  many  avenues  for  commercialization.  For  Navy  needs 
close  cooperation  with  ARL  at  Penn  State  is  well  established  and  joint  programs  are  now 
extending  to  NUSC  through  Ewart  and  Amin. 
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The  existence  of  domains  is  essential  in  many  practical  applications  of  ferroics.  Here  we  dis¬ 
cuss  devices  in  which  &  fixed  spatial  distribution  of  domains  plays  the  significant  role. 
Depending  on  the  prevailing  attributes  of  multidomain  single  crystals,  three  different  possi¬ 
bilities  can  be  distinguished.  In  domain-geometry-engineered  samples  the  spatial  distribution 
of  domains  is  tuned  to  correspond  to  the  ^-vectors  of  fields  propagating  through  the  material. 
In  domain-average-engineered  samples  the  crystal  is  subdivided  into  a  very  large  number  of 
domains,  representing  a  limited  number  of  domain  states.  In  domain-wall-engineered  sam¬ 
ples  the  characteristics  of  static  walls  can  play  an  essential  role  in  the  averaged  macroscopic 
properties.  Examples  illustrating  these  approaches  are  given. 

Keywords:  domain  engineering;  domain-geometry  engineering;  domain-average-engineer¬ 
ing;  domain-wall-engineering;  multidomain  ferroics;  static  domain  pattern 


1.  INTRODUCTION 

Practical  applications  of  ferroics  (undergoing  a  phase  transition  from 
the  point  group  G  to  F)  are  of  two  basically  different  characters:  those 
which  rely  on  properties  of  single  domain  samples  and  those  in  which 
the  presence  of  domains  is  essential.  The  latter  can  be  categorized  into 
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devices  based  on  dynamical  domain  processes  or  those  in  which  a  static 
distribution  of  domains  plays  a  significant  role.  Here  we  discuss  the  last 
mentioned  case,  general  characteristics  of  multidomain  samples  with 
fixed  spatial  distribution  of  domains. 

Depending  on  the  prevailing  attributes  of  multidomain  single  crystal 
samples,  several  possibilities  can  be  distinguished.  Those  discussed  in 
the  following  sections  2  to  4  differ  in  general  features  of  the  geometry 
of  domains  and  it  is  assumed  that  domain  walls  are  of  negligible 
thickness.  In  the  last  considered  case,  the  thickness  of  domain  walls  is 
finite  and  specimens  exhibit  large  wall  density.  Aspects  of  preparation 
and  properties  of  such  samples  are  discussed  in  the  section  5. 


2.  DOMAIN-GEOMETRY-ENGINEERED  FERROIC  SAMPLES 

As  early  as  in  1964,  Miller1'1  showed  that  a  regular  pattern  consisting  of 
180  domains  in  BaTi03  with  a  period  corresponding  to  the  coherence 
length  could  substantially  increase  the  effectiveness  of  optical  second 
harmonic  generation.  The  idea  is  based  on  two  factors:  a)  domain  states 
differ  in  the  sign  of  nonlinear  optical  coefficient,  b)  the  width  of 
domains  is  tuned  to  the  coherence  length.  It  was  Feng  Duan  ct  al.  ^ 
who  succeeded  in  manufacturing  a  periodic  domain  pattern  in  LiNbCb 
311(1  P*°ved  its  efficiency  in  nonlinear  optics.  More  recently,  S.  N.  Zhu 
et  al.  initiated  an  essential  progress  in  this  field  by  producing  a 
domain  pattern  whose  geometry  corresponds  to  a  Fibonacci  superlattice 
and  which  makes  it  possible  to  realize  second  harmonic  generation  for 
multiple  wavelengths.  But  the  significance  of  domain  patterns  with 
engineered  geometry  was  proved  also  in  acoustics.  Meeks  at  el.{4] 
produced  tunable  acoustic  systems  based  on  spatial  modulation  of 
elastic  coefficients  in  a  periodic  domain  pattern  in  NdP50i4  while  Y.  Y. 
Zhu  et  al.  succeeded  in  producing  transducers  up  to  800  MHz  based 
on  spatial  modulation  of  piezoelectric  coefficients  in  multidomain 
LiNb03. 

In  all  these  applications  the  periodicity  of  domain  patterns 
corresponds  to  ^-vectors  of  propagating  waves  and  the  multidomain 
domain  systems  represent  just  two  domain  states.  Up  to  now  less  than 
10  ferroics  have  been  utilized  in  this  area,  and  in  any  of  these  materials 
the  total  number  of  domain  states  v=  |G|  /  |F|  =  2  (here  |G|  and  |F|  stand 
for  the  order  of  the  parent  and  ferroic  phase  point  groups,  resp.).  This  is 
understandable  since  to  produce  a  regular  pattern  in  ferroics  with  v  >  2 
in  which  only  two  domain  states  are  involved  is  not  trivial. 
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To  solve  this  task  in  particular  cases  experimentally  is,  however, 
possible.  In  fact  the  nature  itself  shows  that  patterns  with  a  regular 
geometry  of  a  limited  number  \i  <  v  of  domain  states  can  be 
materialized.  Several  observations  have  been  made  of  domain  patterns 
with  regular  geometry,  fulfilling  this  requirement.  As  an  example,  we 
refer  to  the  Forsbergh’s  square-net  pattern161  shown  in  Fig.l,  which  has 
been  repeatedly  observed  in  BaTiCb  single  crystals  as  well  as  ceramic 
grains.  It  might  be  inspiring  to  examine  its  macroscopic  properties  for 
external  fields  with  both  k  =  0  and  k  ^  0,  based  on  its  symmetry 
characteristics. 


FIGURE  1  Example  of  a  “natural”  multiaxial  domain-geometry 
engineered  system:  the  Forsbergh’s  square-net  pattern  in  tetragonal 
BaTi03.  Left:  microscopic  picture;  right:  arrangement  of 
tetrahedral  building  blocks.  From  Ref.  [6]. 

Several  other  multiaxial  threedimensional  and  reproducible  patterns 
have  been  observed  under  natural  conditions  and  there  is  little  doubt 
that  they  could  be  produced  artificially.  To  mention  just  one  more 
example,  we  refer  to  the  Arlt’sm  patterns  a  and  P;  till  now  they  have 
been  observed  only  in  ceramic  grains  but  very  probably  they  could  be 
created  in  crystals  under  properly  designed  external  forces.  Again,  in 
order  to  deliberate  about  their  properties,  the  symmetry  analysis  would 
be  the  first  step  to  take. 

Indeed,  regular  domain  systems  can  offer  unexpected  symmetry 
properties.  Thus,  for  instance,  the  well  known  Hallbach  array  of 
magnets  with  asymmetric  distribution  of  magnetic  field  can  have  a 
simple  analogy  in  thin  ferroelectric  plates  containing  a  regular  system 
of  90°  domain  pairs.  While  in  the  latter  case,  because  of  the  existence  of 
free  charges,  we  do  not  expect  pronounced  external  depolarizing  field 
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effects,  macroscopic  properties  of  the  array  could  offer  new  applicable 
aspects. 


3.  DOMAIN- AVERAGE-ENGINEERED  FERROIC  SAMPLES 

In  contrast  to  the  previous  cases,  in  domain-average-engineered 
samples  of  ferroic  crystals  the  specimen  is  subdivided  into  a  very  large 
number  of  domains,  representing  just  p  <  v  domain  states.  Such 
situations  can  be  achieved  by  cooling  samples  through  their  phase 
transition  temperatures  under  properly  oriented  stresses,  electric  fields 
or  their  combinations.  The  geometry  of  domains  is  irregular.  The 
sample’s  response  to  external  fields  is  roughly  described  by  tensorial 
properties  averaged  over  all  involved  domain  states.  Thus,  e.g.  for  the 
piezoelectric  response  we  can  write  in  the  zeroth  approximation 

sJt=dSkE„ 

y  a= 1 

Here  V(a)  denotes  the  volume  occupied  by  the  domain  state  a. 

In  recent  years.  Park  et  al.[8,9]  ,  Yin  and  Cao1101  and  other  authors 
considered  a  case  of  this  character  to  discuss  piezoelectric  properties  of 
PZN-PT  single  crystals  poled  along  one  of  the  {001}  directions. 
Assuming  that  the  material  went  through  the  phase  transition  from  G  — 
m  3  m  to  F  =  3m,  the  poling  along  [100]  supports  the  coexistence  of  four 
domain  states  with  spontaneous  polarization  along  the  directions  [111], 
[1 T 1],  [1 1 T]  and  [1 T  T],  with  equal  probability.  In  this  statement,  it  is 
assumed  that  the  domain  wall  orientation  aspects  (e.g.  mechanical 
compatibility)  can  be  neglected.  In  fact,  in  samples  of  ferroelastic 
crystals  strictly  speaking,  only  one  set  of  mechanically  permissible 
parallel  domain  walls  is  allowed^1  ^  while  in  real  samples  walls  of 
various  orientations  coexist  connected  with  additional  elastic  strains, 
paid  by  increased  elastic  energy. 

The  symmetry  aspects  of  domain-average-engineered  samples  can 
be  discussed  in  a  general  way.  The  task  is  to  determine  the  average 
point  symmetry  H,  i.e.  the  symmetry  of  the  subset  of  domain  states 
contained  in  the  multidomain  sample.  This  was  addressed  by  Fousek  et 
alJ12^  and  the  procedure  can  be  facilitated  by  the  use  of  the  computer 
programme  of  Schlesmann  and  LitvinJ13^ 

As  an  example,  in  Fig.2  we  reproduce^12*  the  average  symmetry 
groups  H  of  six  selected  subsets  of  domain  states  (out  of  13  leading  to 
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different  symmetries  F)  which  arise  in  the  phase  transition  specified 
above  (here  pi  =  8).  States  are  represented  by  self-explanatory  numbers 
and  dots  indicating  the  directions  of  Ps  vectors. 


Representative 
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FIGURE  2  Examples  of  subsets  of  domain  states  corresponding  to 
the  transition  m3 m  —  3m  and  their  symmetries. 

Taking  into  account  the  distribution  of  polar  vectors  and 
corresponding  strain  tensors,  one  can  determine  which  external  forces 
should  be  applied  in  order  to  obtain  any  of  the  domain-average- 
engineered  systems.  The  trivial  example  is  the  subset  [1]  produced  by 
the  electric  field  E  along  [111].  The  subset  [1368],  discussed18'1  1 
before,  will  be  produced  by  the  field  along  [001]  while  the  combination 
[15]  requires  the  application  of  a  uniaxial  stress  along  [111].  The  subset 
[13]  calls  for  the  application  of  both  electric  field  along  [001]  and 
uniaxial  stress  along  [1  10]. 

In  addition  to  the  applied  fields,  electrical  or  mechanical,  there  are 
other  approaches  which  can  eliminate  particular  domain  states  in  a 
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given  ferroic  material  The  stimulating  example  is  based  on  chiral 
dopants.  Keve  et  aJ.  showed  that  doping  TGS  with  Z-a-alanine 
prefers  just  one  of  the  two  domain  states  with  antiparallel  Ps.  Zikmund 
and  Fousek  1  generalized  this  approach  and  showed  that  chiral 
substitients  can  reduce  the  number  of  domain  states  by  a  factor  of  two 
in  a  number  of  ferroic  species. 

It  is  understood  that  m  the  mentioned  approaches  we  leave  behind 
the  problems  of  domain  coexistence  connected  with  their  mechanical 
compatibility.  It  seems  obvious  that  domain  average  engineering  can 
succesfully  lead  to  formation  of  crystalline  systems  with  new  desired 
properties  m  particular  in  crystals  where  the  domain  size  is  small. 
Crystalline  systems  exhibiting  tweed  microstructures  similar  to  those 
observed  m  La-modified  lead  titanate^.or  PLZTll8J  might  serve  as 
candidates  for  this  approach. 


4.  PHASE-  AND-  DOMAIN-AVERAGE  ENGINEERED  FERROICS 

The  domain-average-engng  concept  can  be  generalized  to  systems  in 
which  the  prescribed  domain  states  represent  two  or  more  ferroic 
species.  Such  multiphase  situations  occur  in  PZT  ceramics  near  the 
morphotropic  boundary  and  originate1191  in  concentration  gradients  as 
well  as  in  the  independent  nucleations  of  the  ferroic  phases  since  the 
transition  is  of  the  1st  order.  It  appears  that  PZN-PT  single  crystals  with 
a  pronounced  piezoelectric  response  contain  blocks  of  both  tetragonal 
and  rhombohedml  symmetry™.  In  the  basic  approximation,  the 
piezoelectric  coefficient  of  a  properly  poled  sample  is  then  described  by 


d ijk 


I 

a=\  lJ 


L  dU?vW 


.  ViYm1* 

M ijk  (walls)  +  A diJk  ( phase  bndrs ) 

where  pi  and  p2  are  numbers  of  domain  states  represented  in  the  phase 
I  and  2,  resp.,  after  poling.  The  high  piezoelectric  response  of  PZN-PT 
smgle  crystals  originates  in  the  combination  of  intrinsic  coefficients 


A 


ijk  35  as  in  the  extrinsic  contributions  due  to  the  motion 


d%> 

ijk 

of  domain  wails  and  phase  boundaries.  The  question  what  role  is  played 
by  any  of  these  contributions  is  still  to  be  solved. 
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When  addressing  the  problem  of  average  symmetries  of  phase-and- 
domain-average-engineered  samples  one  can  follow1213  a  similar 
approach  as  mentioned  above.  Consider  that  two  species  coexist, 
namely  m  3  m-3m  and  m  3  m-4mm.  It  can  be  shown  that  when  poling 
along  principal  directions,  i.e.  E  ||  [001]  or  E  ((  [Oil]  or  E  ||  [111], 
regions  of  the  two  species  differ  in  the  systems  of  Ps  vectors  but  are  of 
the  same  averaged  symmetries  4 mm,  mm2  or  3m,  resp. 

It  has  to  be  stressed  that  in  both  domain-average  enginnered  and 
phase-and-domain-average  engineered  systems,  electrical  and 
mechanical  compatibility  conditions  play,  in  the  energy  evaluations,  a 
significant  role.  It  is  beyond  the  scope  of  this  presentation  to  discuss 
these  problems  in  detail  and  the  subject  will  be  addressed  in  another 
paper. 


5.  DOMAIN- WALL-ENGINEERING 


In  the  last  considered  case,  the  thickness  of  domain  walls  is  finite  and 
the  specimens  exhibit  large  wall  density.  As  before,  the  representation 
of  specific  walls  can  be  influenced  by  external  forces.  In  such  domain- 
wall-engineered  samples  the  characteristics  of  static  walls  can  play  an 
essential  role  in  the  averaged  macroscopic  properties. 

It  was  predicted  by  Walker  and  Gooding1221  that  Dauphine  domain 
walls  in  nonpolar  quartz  can  carry  a  dipole  moment  and  this  was  later 
demonstrated  experimentally  by  Snoeck  et  al.3233,  in  the 
incommensurate  phase  of  SiC>2.  In  fact  it  is  easy  to  demonstrate  this 
possibility  for  ferroelastic  walls.  If  in  the  free  energy  function  the 
invariant  /%i(Ss}j/3xk)-Ei  is  allowed  by  symmetry,  a  ferroelastic  wall  will 
carry  polarization 


Wyw 


dxh 


Here  the  /Men sor  describes  the  flexoelectric  effect.  As  an  example, 
consider  the  species  m3m  -  4 fmmm  represented  by  crystals  of  SrTiCb, 
or  CsPbCl3  and  the  species  m3m  —  4mm  describing  the  properties  of 
BaTiOj.  In  both  cases,  ferroelastic  walls  (90°  wall  in  the  latter  case)  are, 
by  symmetry,  allowed  to  carry  polarization 

OX 
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this  is  demonstrated  in  Fig.  3.  Symmetries  of  nonferroelastic  domain 
wails  from  which  their  possible  macroscopic  properties  can  be 
envisaged  have  been  discussed  by  PfivratskA  and  Janovec1241. 


x 


FIGURE  3  Due  to  flexoelectricity,  domain  wall  separating  two 
ferroelastic  domains  of  species  m3m  -  4/mmm  or  m3m  -  4mm 
can  cany  polarization  represented  by  the  arrow. 

Should  macroscopic  properties  of  domain  walls  play  a  non- 
negligible  role  in  the  properties  of  a  ferroic  sample,  they  have  to  occupy 
a  sizable  volume.  Realistic  values  of  domain  width  and  domain  wall 
thickness  are  a*  =  1  to  10  pm  and  =  3  nm,  resp.  Suppose  that  the 
required  relative  volume  occupied  by  domain  walls  is  10%  and  that 
there  so  no  way  how  to  increase  the  wall  thickness.  Then  we  look  for 
methods  how  to  reduce  the  average  domain  width  to  about  30  nm  i.e. 
how  to  increase  the  wall  density  30  times  or  more. 

It  appears  that  there  is  a  significant  number  of  experimental  methods 
which  could  address  this  problem.  Here  we  mention  in  passing  only 
some  of  those  which  might  be  considered.  The  density  of  defect- 
induced  order  parameter  gradients  could  be  increased  at  crystal 
growth  .  The  size  of  AFM-written  domains1251  could  be  further 
reduced,  their  density  increased.  Attempts  to  freeze-in  high  density 
tweed  structures  in  some  ferroelastics  close  above  Tc  or  high  density 
pattern  of  discommensuradons  in  modulated  phases  close  above  the 
lock-in  temperatures^71  appear  attractive.  Fixing  photorefiactive 
gratings  by  domains1"'  in  ferroelectrics  with  small  ps  (small  defect- 
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assisted  nucieation  energies  to  produce  domains  in  real  time)  looks  also 
promissing. 


6.  CONCLUSIONS 

Till  now,  the  field  of  domain  geometry  engineering  has  been 
suecesfully  developed  and  applied  to  ferroics  representing  five  species 
of  ferroelectrics  and  one  species  of  ferroelastics,  all  with  two  domain 
states.  Such  specimens  proved  to  be  competitive  in  die  field  of 
nonlinear  optics  and  promissing  in  the  field  of  ultrasonic  generation  and 
detection.  But  much  more  complex  domain  systems  are  obtainable  and 
have  not  yet  been  considered  and  investigated. 

Two  new  areas  of  research  appear  to  be  very  promissing.  Domain- 
and  phase-average  engineered  systems  offer  increased  values  of 
macroscopic  properties,  probably  strengthened  by  domain  wall  and/or 
phase  boundary  induced  motions.  We  still  miss  detailed  date  about  the 
real  structure  of  existing  compounds  like  PZN-PT  as  well  as  theoretical 
analysis  of  multiple  domain  states  compatibility,  the  more  so  for 
multiple  phase  systems.  Domain-wall  engineered  samples  with  high 
density  of  walls  also  promise  a  new  interesting  research  and  application 
area.  At  present,  particular  systems  useful  in  selected  applications  could 
be  specified  and  methods  to  produce  high  density  domain  patterns 
investigated  experimentally. 
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Abstract 

Multidomain  samples  of  ferroies  (ferroelectrics,  ferroelastics,  and  related 
materials)  with  fixed  geometrical  distribution  of  domains  can  offer  new 
macroscopic  properties  required  for  particular  applications.  Two  extreme  cases 
of  such  applications  are  defined.  In  domain-geometry-engineered  samples 
of  ferroic  crystals,  the  spatial  distribution  of  domains  and  thus  the  spatial 
distribution  of  tensorial  properties  is  tuned  to  correspond  to  the  ^-vectors  of 
applied  electric,  optical  or  acoustic  fields.  For  a  given  wavelength,  the  size, 
geometry,  and  distribution  of  domains  give  rise  to  a  qualitatively  new  kind  of 
response  specified  by  the  symmetry  of  the  multidomain  system.  In  domain- 
average-engineered  samples  of  ferroic  crystals,  the  specimen  is  subdivided 
into  a  very  large  number  of  domains,  representing  pi  domain  states  where  pu 
is  smaller  than  the  theoretically  allowed  maximum  number,  and  forming  a 
regular  or  irregular  pattern.  Its  response  to  external  fields  is  roughly  described 
by  tensorial  properties  averaged  over  all  of  the  domain  states  involved.  The 
effective  symmetry  of  the  domain-average-engineered  system  is  given  by  a 
point  group  H  and  we  show  how  it  can  be  determined.  As  an  example,  all 
groups  H  are  specified  for  domain-average-engineered  samples  which  can 
arise  in  a  material  undergoing  the  phase  transition  with  symmetry  change  from 
m3m  to  3m. 


Ferroic  materials  (and  here  we  concentrate  on  non-magnetic  materials,  i.e.  on  ferroelectrics, 
ferroelastics,  and  higher-order  ferroies)  play  an  essential  role  in  a  number  of  technical  appl¬ 
ications,  In  some  of  them,  dynamic  domain  processes  are  essential  (e.g.  thin-film  memories, 
electron  emitters)  while  in  others  the  static  distribution  domains  in  the  sample  play  the  crucial 
role.  In  this  contribution  we  concentrate  on  the  latter  case  and  wish  to  specify  a  clear  distinction 
between  two  kinds  of  such  static  multidomain  system. 

We  have  in  mind  materials  undergoing  a  structural  phase  transition  from  the  parent  phase 
of  point  group  G  into  the  ferroic  phase  of  symmetry  F  C  G  (such  a  material  is  referred  to 
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as  belonging  to  the  species  G  —  F).  This  leads  necessarily  to  the  possibility  of  a  coexistence 
of  v  domain  states  [1];  v  =  |G|/|F|  where  \A\  is  the  order  of  the  group  A.  Assuming  that 
domain  walls  are  of  negligible  thickness  compared  with  the  size  of  the  domains,  there  are  two 
different  ways  of  specifying — and  utilizing — properties  of  a  multidomain  sample  with  &  fixed 
distribution  of  domains;  we  propose  to  refer  to  them  as  domain  geometry  engineering  and 
domain  average  engineering ,  respectively.  It  is  the  purpose  of  this  contribution  to  give  their 
definitions  and  in  particular  to  specify  symmetry  properties  of  domain- average-engineered 
multidomain  samples. 

First,  we  consider  domain  geometry  engineering.  Consider  a  multidomain  sample  for 
which  the  geometry  of  the  spatial  distribution  of  domains  and  therefore  that  of  the  tensorial 
material  coefficients  is  specified.  Macroscopic  responses  of  such  samples  to  external  fields 
(forces)  of  defined  frequency  are  determined  by  this  distribution.  If  the  applied  fields  are 
static,  the  response  of  a  multidomain  sample  is  primarily  determined  by  the  spatial  distribution 
of  tensorial  properties  (domains)  and  of  the  applied  field,  and  codetermined  by  the  boundary 
conditions  along  domain  walls.  (Only  in  cases  of  the  simplest  geometry  of  domains,  namely 
a  single  system  of  parallel  domain  walls  in  the  case  of  ferroelastic  species,  can  the  latter  be 
eliminated.)  As  an  example,  consider  a  multidomain  piezoelectrically  active  sample.  The 
spatial  distribution  of  strain  is 

Sjk(r)=dijk(r)E,(r)  (la) 

diJk(r)  =  (lb) 

where  dtjk  stands  for  the  piezoelectric  tensor,  the  factor  f^a)(r)  =  0  or  1,  and  a  denotes 
the  domain  state:  a  =  1, 2, . . . ,  /z  with  /z  ^  v.  Domain  geometry  engineering  related  to 
dynamic  external  fields  is  of  particular  interest.  The  /:-vector  of  the  applied  fields  defines 
the  wavelength  whose  magnitude  is  chosen  to  be  appropriately  related  to  the  size  of  domains 
and  whose  direction  is  correlated  with  the  domain  geometry.  For  a  given  wavelength,  the 
size,  geometry,  and  distribution  of  domains  give  rise  to  a  qualitatively  new  kind  of  response 
specified  by  the  symmetry  of  the  multidomain  system.  In  the  example  specified  above,  we 
expect  the  presence  of  new  piezoelectric  resonance  frequencies.  It  was  this  case  which  was 
suggested  by  Newnham  et  al  [2],  offering  new  resonance  modes  of  a  two-domain  sample. 
Recently,  more  involved  multidomain  piezoelectric  systems  were  suggested  and  realized  [3] 
in  crystals  of  LiNbC>3  and  LiTa03.  Referred  to  as  acoustic  superlattices,  they  can  be  used  to 
generate  and  detect  ultrasonic  waves  with  frequencies  in  the  range  up  to  several  hundreds  of 
MHz.  Another  example,  which  has  received  unusual  attention,  is  quasi-phase-matched  optical 
multipliers.  When  the  conventional  phase- matching  condition  (n^  =  nw)  cannot  be  realized 
in  a  particular  material  because  of  unsuitable  dispersion  of  refractive  indices,  often  a  quasi- 
phase-matched  system  can  be  constructed  which  offers  a  high  integrated  non-linear  optical 
response  leading  to  frequency  doubling  [4].  This  requires  that  a  periodic  domain  pattern 
be  fabricated  with  a  period  twice  the  coherence  length  lc.  Such  domain- shape-engineered 
systems  are  now  widely  used.  An  even  more  intricate  geometry-engineered  domain  pattern  has 
been  designed  [5]  in  which  two  geometrically  different  building  blocks  A,  B,  each  containing 
two  domains  with  antiparallel  spontaneous  polarizations,  are  arranged  to  form  a  Fibonacci 
sequence.  This  leads  to  the  possibility  of  second-harmonic  light  generation  simultaneously 
for  several  optical  frequencies. 

While  domain-geometry-engineered  systems  have  been  repeatedly  realized  and  theor¬ 
etically  analysed,  the  alternative  approach  to  studying  and  utilizing  multidomain  ferroic 
samples  with  static  chaotic  distribution  of  domains  has  only  recently  become  extremely 
attractive.  By  the  term  domain  average  engineering  we  mean  a  situation  in  which  the  ferroic 
sample  is  subdivided  into  a  very  large  number  of  domains,  representing  /jl  domain  states  where 
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fi  <  v,  and  forming  a  regular  or  irregular  pattern.  Ideally,  the  domain  size  is  expected  to  be 
much  smaller  than  the  wavelength  of  externally  applied  fields.  Here,  in  contrast  to  the  case 
for  domain  shape  engineering,  the  spatial  distribution  of  tensorial  material  coefficients  is  not 
defined  or  is  irrelevant.  The  response  of  the  sample  can  be  to  some  approximation  described 
by  tensorial  properties  averaged  over  all  domain  states  involved.  Considering,  as  an  example, 
again  a  sample  consisting  of  piezoelectrically  active  domains,  we  expect  for  the  multidomain 
sample 


£jk  —  dijkEi 

(2a) 

ij*  =  ^it€via) 

a~\ 

fi  <  V. 

(2b) 

Both  equations  (1),  (2)  can  be  easily  generalized  for  tensors  of  higher  order. 

Recently,  a  case  of  this  character  was  considered  [6-8]  to  allow  discussion  of  piezoelectric 
properties  of  PZN-PT  single  crystals  poled  along  one  of  the  {001}  directions.  Assuming  that 
the  material  went  through  the  phase  transition  from  G  =  m3m  to  F  —  3m,  poling  along 
[100]  supports  the  coexistence  of  four  domain  states  with  spontaneous  polarization  along 
the  directions  [111],  [111],  [111],  and  [ill],  with  equal  probability.  In  this  statement,  it 
is  assumed  that  the  domain  wall  orientation  (i.e.  mechanical  compatibility)  aspects  can  be 
neglected.  In  fact.  In  samples  of  ferroelastic  crystals,  strictly  speaking,  only  one  set  of 
mechanically  permissible  parallel  domain  walls  is  allowed  [9]  while  in  real  samples  walls 
of  various  orientations  coexist  connected  with  additional  elastic  strains,  paid  for  by  increased 
elastic  energy  [10]. 

We  now  discuss  the  symmetry  of  domain-average-engineered  samples  in  which  the 
volumes  of  the  domain  states  represented  are  identical,  i.e.  V(a)  =  V  jfi.  Such  situations  can 
be  achieved  by  cooling  samples  through  their  phase  transition  temperatures  under  properly 
oriented  stresses,  electric  fields  or  combinations  of  these.  We  introduce  a  classification  of 
domain-average-engineered  ferroic  samples  and  determine  their  average  point  symmetries. 
This  average  symmetry  is  taken  to  be  the  symmetry  of  the  subset  of  domain  states  contained 
in  the  multidomain  ferroic  sample. 

Consider  the  phase  transition  from  G  to  F.  The  symmetry  analysis  is  based  on  the  coset 
decomposition  of  the  point  group  G  with  respect  to  its  subgroup  F,  I.e. 

G  —  F  +  g2F  +  g3F  +  *  *  *  +  gyF 

where  the  elements  gt  are  the  coset  representatives  of  the  decomposition  and  g\  =  1.  We 
denote  the  v  domain  states  which  may  arise  at  the  transition  as  Si ,  S2, . . . ,  Sp.  The  symmetry 
groups  Fi  of  the  domain  states  and  the  relative  orientations  of  the  domain  states  and  their 
polarizations  P;  are  all  determined  by  the  coset  representatives,  i.e.  Fx  —  F,  F;  =  giF\g~\ 
Si  =  giSi ,  and  Pt  ~  gt*  P\ ,  i  =  2, 3, ...» v.  The  closure  of  the  group  G  under  multiplication 
implies  a  permutation  of  the  cosets  of  the  coset  decomposition  and  in  turn  a  permutation  of 
the  domain  states  S;  under  elements  g  of  G.  The  action  of  an  element  g  of  G  on  Sj  is  defined 
as  gSt  —  ggiS\  =  gjfSi  —  g j S j  =  Sj,  where  /  is  an  element  of  F,  and  the  domain  state  St* 
is  transformed  by  the  element  g  into  the  domain  state  Sj .  The  action  of  an  element  g  of  G  on 
a  subset  of  domains  is  denoted  by  g[Si ,  S2, . . . ,  S^}  =  {gS\ ,  gS2, . . . ,  gSM}. 

Two  subsets  of  domains  {Si,  52,...,  SM]  and  {SJ,  S^, . . . ,  S^}  are  said  to  belong  to  the 
same  class  of  subsets  of  domains  if  there  exists  an  element  g  of  G  such  that  g  {5! ,  5s, . . . ,  SM}  = 
{g Su  gS2f . , . ,  g5M]  =  {SJ,  S£, . . . ,  Sp.  The  symmetry  group  if  of  a  subset  of  domains 
{Si,  52, ,  SM}  is  defined  as  the  group  of  all  elements  g  of  G  which  leave  the  set  invariant, 
i.e.  g{Si,  S2, . . . ,  SM}  =  {Si,  S2, . . . ,  SM).  The  group  H  represents  the  effective  symmetry  of 
the  domain-average-engineered  system  consisting  of  the  subset  of  domains  {Si ,  S2, . . . ,  SM}, 
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As  an  example,  we  consider  the  phase  transition  from  G  =  m3m  to  F  =  3xyzmiy.  Here 
v  =  8.  The  indexing  of  the  domain  states,  the  corresponding  coset  representatives  of  the 
coset  decomposition  of  G  with  respect  to  F ,  the  symmetry  groups,  and  the  corresponding 
polarizations  in  each  domain  state  are  given  in  table  1. 

Table  1.  Domain  state  index,  coset  representative,  symmetry  group,  and  polarization. 


Index  i 

Coset  representative  g, 

Pi  —  8iP\  g~ 1 

Pi=giPl 

1 

1 

3xyzmxy 

(A,  Ay  A) 

2 

2X 

3xyzmyz 

1 

1 

3 

2j 

3xyzmxy 

(-A,  —Ay  A) 

4 

2, 

2xyzmxz 

{-Ay  Ay -A) 

5 

I 

3xyzmxy 

{-Ay  -Ay  -A) 

6 

mx 

3xyzmyz 

(—A,  Ay  A) 

7 

mz 

3 xyzmxy 

(A,  A,  —A) 

8 

rtly 

^xyzmxz 

(A,  -Ay  A) 

All  subsets  of  these  domain  states  have  been  classified  into  classes  as  defined  above.  In 
table  2  we  list  one  subset  of  domain  states  from  each  class.  Each  subset  is  denoted  by  listing, 
between  square  brackets,  the  indices  of  the  domain  states  contained  in  that  subset,  the  indices 
having  been  given  in  table  1,  e.g.  the  subset  {5j,  53,  S5 )  is  denoted  by  [135],  In  the  right-hand 
column  is  the  subgroup  H  of  elements  of  G  which  leave  the  corresponding  subset  invariant. 
This  table,  in  fact,  represents  the  list  of  domain-average-engineered  systems  which  can  arise 
in  a  material  undergoing  a  phase  transition  from  m3m  to  3m. 

Table  2.  Representative  subsets  of  domain  states  for  the  species  m3m  —  3 m  and  the  subgroups  of 
m3m  which  leave  them  invariant. 


Representative  subset 

Symmetry  H  of  the  subset 

[1]  or  [2345678] 

3jc yzmxy 

[13]  or  [245678] 

IKxytoxy  2r 

[15]  or  [234678] 

3 xyzmxy 

[16]  or  [234578] 

mxmyZ2yZ 

[123]  or  [45678] 

3  xyzmxz 

[135]  or  [24678] 

rtlxy 

[136]  or  [24578] 

mxy 

[1234] 

43  m 

[1235] 

Wxz 

[1238] 

3xyzmxz 

[1356] 

2 Mz 

[1357] 

TTl  x  y  Ffl  jj  y  fit  ▼ 

[1368] 

4  zmxmxy 

In  figure  1 ,  for  each  subset  listed  in  table  2,  we  schematically  represent  the  array  of  domain 
states  and  their  polarizations  associated  with  the  domain  states  of  each  subset.  Each  domain 
state  is  denoted  by  a  heavy  dot  at  a  comer  of  the  cube.  This  represents  a  polarization  from 
the  centre  of  the  cube  to  that  comer — that  polarization  given  in  table  1  associated  with  the 
corresponding  domain  state.  Subfigure  [1]  denotes  the  single-domain  state  with  polarization 
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Figure  1.  Graphical  representation  of  the  subsets  of  domain  states  whose  symmetries  are  specified 
in  table  1.  Points  at  the  cube  vertices  represent  spontaneous  polarization  vectors;  the  origin  is  in 
the  centre  of  the  cube.  The  numbering  is  that  of  the  indices  of  the  domains  and  polarizations  given 
in  table  1. 


in  the  [111]  direction.  The  figure  denoted  by  [1368],  e.g,,  denotes  a  multidomain  sample  in 
which  the  following  polarization  vector  directions  are  equally  represented:  [111],  [Hi],  [111], 
and  [111],  The  corresponding  symmetry  groups  of  all  these  multidomain  systems  are  listed, 
as  already  pointed  out,  in  the  right-hand-side  column  of  table  2. 

Taking  into  account  the  distribution  of  polarization  vectors  and  corresponding  strain 
tensors,  one  can  determine  which  external  forces  should  be  applied  in  order  to  obtain  any 
of  the  domain-average-engineered  systems  listed  in  table  1.  The  trivial  example  is  the  system 
[1]  produced  by  the  electric  field  E  along  [111],  The  system  [1368],  discussed  above  (see 
references  [6-8])  will  be  produced  by  die  field  along  [001]  while  the  system  [16]  requires  the 
application  of  the  field  along  [Oil],  Hie  combination  [15]  requires  the  application  of  a  uniaxial 
stress  along  [1 1 1]  while  the  system  [13]  calls  for  the  application  of  both  an  electric  field  along 
[001]  and  a  uniaxial  stress  along  [1 10], 

It  is  understood  that  in  this  symmetry  approach  we  leave  behind  problems  of  coercive 
fields  and  stresses  as  well  as,  as  already  mentioned,  problems  of  domain  coexistence  connected 
with  their  mechanical  compatibility.  It  seems  obvious  that  domain  average  engineering  can 
successfully  lead  to  the  formation  of  crystalline  systems  with  new  desired  properties,  in 
particular  in  crystals  where  the  domain  size  is  small. 

Each  of  the  methods  of  domain  engineering  specified  above  can  open  a  new  vista  of 
materials  research  possibilities  in  die  area  of  ferroic  materials  and  lead  to  multidomain 
assemblies  with  new  desired  properties. 


L38 


Letter  to  the  Editor 


JF  appreciates  the  support  of  the  Ministry  of  Education  of  the  Czech  Republic  (Project 
No  VS96006)  and  the  assistance  of  the  Materials  Research  Laboratory,  The  Pennsylvania  State 
University;  DBL  acknowledges  the  hospitality  of  Dr  V  Kopsky  and  the  Institute  of  Physics  of 
the  Czech  Academy  of  Sciences,  and  the  support  of  the  Czech  Ministry  of  Education  under 
Grant  ME336(1999)  and  of  the  National  Science  Foundation  under  Grants  No  DMR-9722799 
and  No  DMR  0074550. 

References 

[  1]  Aizu  K  1970  Phys.  Rev.  B  2  754-72 

[2]  Newnham  R  E,  Miller  C  S,  Cross  L  E  and  Cline  T  W  1975  Phys.  Status  Solidi  a  32  69-78 

[3]  Chen  Y  F,  Zhu  A  N,  Zhu  Y  Y,  Ming  N  B,  Jin  B  B  and  Wu  R  F  1997  Appl  Phys .  Lett.  70  592-4 

[4]  Duan  F,  Ming  N  B,  Hong  J  F,  Yang  Y  S,  Zhu  J  S,  Yang  Z  and  Wang  Y  N  1980  Appl  Phys.  Lett.  37  607-9 

[5]  Zhu  S  N,  Zhu  Y  Y,  Qin  Y  Q,  Wang  H  F,  Ge  C  Z  and  Ming  N  B  1997  Phys.  Rev.  Lett.  78  2752-5 

[6]  Park  S  E  and  Shrout  T  R  1997  J.  Appl.  Phys.  82  1804-1 1 

[7]  Wada  S,  Park  S  E,  Cross  L  E  and  Shrout  T  R  1999  Ferroelectrics  221  147-55 

[8]  Yin  J  and  Cao  W  2000  J.  Appl  Phys.  87  7438^1 

[9]  Fousek  J  and  Janovec  V  1969  J.  Appl  Phys.  40  135-42 

[10]  Salje  E  K  H  and  Ishibashi  Y  1996  J.  Phys.:  Condens.  Matter  8  8477-95 


APPENDIX  3 


Domain  and  Phase  Change  Contributions  to 
Response  in  High  Strain  Piezoelectric  Actuators 

L.  Eric  Cross 


Evan  Pugh  Professor  Emeritus  of  Electrical  Engineering 
Materials  Research  Laboratory 
The  Pennsylvania  State  University 
University  Park.  PA  16S02  USA 


Abstract.  Current  solid  state  actuators  are  briefly  compared  to  traditional  actuator 
technologies  to  highlight  the  major  need  for  enhanced  strain  capability.  For  the 
ferroelectric  piezoelectric  polycrystal  ceramics,  the  balance  of  evidence  suggests  a 
large  entrinstc  contribution  to  the  field  induced  strain  from  ferroelectric- 
ferroelastic  domain  wall  motion.  Here-to-fore  the  intrinsic  single  domain 
contribution  has  been  derived  indirectly  from  phenomenological  analysis.  Now, 
new  evidence  of  a  stable  monoclinic  phase  at  compositions  very  close  to  the  MPB 
suggest  that  the  previous  assessment  will  need  to  be  revised. 

Actuator  behavior  in  the  new  iead  zinc  niobate-lead  titanate  (PZN:PT) 
single  crystal  show's  most  unusual  anisotropic  behavior.  For  1 1 1  oriented  field 
poled  crystals  in  the  rhombohedral  phase  normal  low  induced  strain  is  observed. 
For  001  field  poled  crystals  however  massive  (0.6%)  quasi-linear  anhysteritic 
strain  can  be  induced.  Since  the  001  oriented  field  in  the  rhombohedral  phase  can 
not  drive  ferroelastic  domain  walls  it  is  suggested  that  the  strain  must  be  intrinsic. 
The  suggestion  is  that  it  is  due  to  an  induced  monoclinic  phase  in  which  the  Ps 
vector  tilts  under  increasing  field  up  to  more  than  20*  from  111,  before  the  vector 
switches  to  the  tetragonal  001  direction.  Such  a  polarization  rotation  mechanism 
has  also  been  suggested  by  Fu  and  Cohen  [1].  Calculations  of  induced  single 
domain  strain  using  measured  eiectrostriction  constants  agree  well  writh  observed 
behavior. 

Recent  measurements  by  Park  et  al  [2]  and  Wada  et  al  [3]  on  single  crystal 
BaTiOs  show  strongly  enhanced  piezoelectricity  at  temperatures  near  the 
ferroelectric  phase  transitions.  Of  particular  relevance  is  the  inverse  experiment 
forcing  the  tetragonal  over  to  the  rhombohedral  phase  with  high  111  oriented 
field.  From  this  result  it  is  suggested  that  both  cubic  and  dodecahedral  mirrors 
participate  in  the  reorientation  through  orthorhombic  to  the  rhombohedral  state 
giving  rise  to  different  value  of  the  induced  dgj  at  different  field  levels. 
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INTRODUCTION 


Characteristics  of  stress  levels,  strain  capability  and  efficiency  for  a  number  of 
actuator  systems  are  summarized  in  Table  1.  Traditional  electromagnetic,  pneumatic  and 
hydraulic  systems  have  high  efficiency  good  stroke  (strain  levels)  and  adequate  force. 
They  are  however  bulky,  rather  slow,  and  require  extensive  backup  generator  equipment. 
Solid  state  electrical  piezoelectric,  electrostrictive  and  phase  switching  actuators  clearly 
have  adequate  force,  tolerable  efficiency,  however  the  stroke  permitted  by  induced  strain 
levels  is  miniscule.  These  materials  are  still  important  as  the  actuation  can  be  very  fast, 
inertia  is  low'  and  compact  systems  can  be  engineered,  however  there  is  still  a  major  need 
for  enhanced  strain  capability. 

In  this  paper,  the  mechanisms  responsible  for  strain  in  the  polycrystal  ferroelectric 
ceramics  will  be  discussed.  Currently  the  balance  of  evidence  points  to  a  major 
contribution  from  ferroelastic  domain  wall  motion  particularly  in  the  soft  PZTs,  but  the 
discovery  of  a  new  monoclinic  phase,  just  at  the  critical  morphotropic  phase  boundary 
composition  requires  a  re-assessment  of  the  intrinsic  contribution  to  response.  Evidence 
from  the  new  lead  zinc  niobatedead  titnate  (PZN:PT)  single  crystals  suggests  that  the 
ultra  high  strain  in  001  field  poled  crystals  may  be  dominantly  intrinsic  involving  the 
field  induction  of  a  monoclinic  phase  in  which  Ps  is  strongly  tilted  in  the  dodecahedral 
mirror  plane  before  finally  switching  to  the  001  polar  tetragonal  form.  Recent 
experiments  on  BaTiC>3  where  the  inverse  experiment,  forcing  the  tetragonal  over  to 
rhombohedral  symmetry  by  high  (111)  oriented  E-field  at  room  temperature  lends  further 
support  for  a  monoclinic  mechanism  for  the  switching,  but  now  involving  both  cubic  and 
dodecahedral  mirrors. 


Polarization  and  Strain  Mechanisms 

For  the  simple  linear  nonferroelectric  insulator  the  induced  deformation  under 
electric  field  may  be  described  by 

xij  =  sijkl^ij  +  dmijEm  +  MmnijEmEn  ( 1 ) 

where  xjj  are  components  of  induced  strain,  Xjj  components  of  applied  stress  Em,  En 
components  of  electric  field  Sjjjj  the  elastic  compliance,  d^j  the  piezoelectric  tensor  and 

Mmnij  the  electrostriction  coefficients.  Obviously  if  the  dmij  are  non  zero,  in  the  absence 
of  stress 


TABLE  1.  Salient  Features  of  Traditional  and  Solid  State  Actuators 


TRADITIONAL  TECHNOLOGIES 

Stress  fMPa) 

Strain 

Effieiencv 

Electromagnetic 

0.02 

0.5 

90% 

Hydraulic 

20 

0.5 

80% 

Pneumatic 

0.7 

0.5 

90% 

SOLID  STATE  ACTUATORS 

Stress  (MPa) 

Strain 

Effieiencv 

Shape  Memory  Alloys 

200 

01 

3% 

Piezoelectric 

35 

0.002 

50% 

Electrostrictive 

50 

0.002 

50% 

Phase  Switching 

50 

0.004 

7 

Eleetroceramie 

Electrostrictive 

0.1 

0.10 

20% 

Polyurethane 

Contractile  Polymer 

03 

0.5 

30% 

Muscle 

035 

0.2 

30% 

which  is  the  piezoelectric  actuation  function  and  if  the  dmjj  are  zero  by  symmetry 

xy  =  MmnijEmEn  (3) 

and  the  system  will  be  electrostrictive. 

In  a  polycrystal  ensemble,  as  for  the  randomly  axed  ceramic,  the  system  is  much 
more  complex.  With  any  symmetry  in  the  crystallites  (fig.  la)  the  random  axial 
arrangement  will  introduce  two  infinite  fold  rotation  axes  («*»)  in  the  macro-ensemble 
symmetry  wiping  out  all  piezoelectricity.  One  of  the  few  way  out  of  this  dilema  is  to 
consider  a  ferroelectric  crystallite  where  in  the  ferroelectric  state  spontaneous  electric 
polarization  is  distributed  in  domains  with  different  equivalent  orientation  states  that  can 
be  reoriented  by  an  external  poling  field  (fig.  lb).  The  remanent  polar  vector  drops  the 
macro-symmetry  to  conical  («m)  which  is  piezoelectric  with  non  zero  d3  j  =  d32,  d33,  d]5 
=  d24- 

All  practical  poiycrystal  ceramic  actuators  use  compositions  in  the  lead 
zircon  ate:  lead  titanate  (PZT)  solid  solution  system  (fig.  1c).  The  reason  for  the 
infatuation  with  PZT  is  the  near  vertical  (morphotropic)  phase  boundary  close  to  the 
52:48  Zr/Ti  cation  ratio,  separating  tetragonal  and  rhombohedral  ferroelectric  phases. 
Qualitatively  two  clear  advantages  stem  from  working  with  compositions  near  this 
boundary,  and  this  is  where  most  practical  systems  are  pitched. 


Unnoied  Poling 


(a)  the  symmetry  is  °°<»m  which  is  centric  and 
forbids  piezoelectricity 

Macro  Symmetry  and  the  importance 


(b)  the  symmetry  is  oom  which  is  non  centric 
(polar)  and  permits  piezoelectricity. 
of  poling  of  a  ferroelectric  species. 


20  30  40  50  60  70  80  90  PbTi03 

MOLE  %  PbTiQ3 

(c)  Lead  zirconate  titanate  phase  diagram. 


POSSIBLE  ORIENTATION  STATES  tN  PER0VSK1TES 


TETRAGONAL  4mm  LBaTi03:  CaTi033 

POLARIZATION  ALONG  — —  K33  -  0.48 
6  EQUIVALENT  <100>  PR  -  8/U.C/CM2 


PbZr  03:  PbTi03 
k33  ~  O.  75 
PR-4qtxC/CM2 


ORTHORHOMBIC  2mm  t  KNb03  :  NoNb03 3 

POLARIZATION  ALONG  ^ —  K33  ~  0.65 
12  EQUIVALENT  <IIO>  pR  -  20/xC/CM2 


RHOMBOHEDRAL  3m 
POLARIZATION  ALONG 
8  EQUIVALENT  <lll> 


(d)  Importance  of  many  oriemation  states  for  polycrvstal  poiability  and  performance. 


FIGURE  1.  Important  features  in  piezoelectric  polycrvstal  ceramic  materials. 
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Firstly,  the  MPB  is  a  first  order  phase  change,  so  that  in  the  co-existence  region  for  the 
two  phases  it  is  possible  in  the  poling  process  to  make  use  of  14  domain  orientation  states 
(fig.  Id)  leading  to  more  complete  poling.  Secondly,  for  properly  chosen  compositions 
one  can  stay  close  to  the  phase  boundary  over  a  wide  temperature  range  and  use  its 
destabilizing  influence  to  enhance  the  intrinsic  polarizability  of  the  single  domain  states. 

From  these  basic  advantages  stem  two  fundamental  problems  in  understanding  the 
behavior.  Firstly,  what  is  the  balance  between  intrinsic  single  domain  and  extrinsic 
contributors  to  response;  and  secondly,  what  exactly  is  the  nature  of  the  dominant 
extrinsic  contributor. 

In  the  randomly  axed  polycrystal  ceramic,  ferroelectric  domain  wall  motion  is 
essential  to  permit  the  induction  of  strong  remanent  polarization  essential  tor 
piezoelectricity.  For  the  quasi-reversible  behavior  working  about  this  induced  remanent 
polar  state  the  question  is  to  what  extent  are  ferroelastic:ferroelectric  wall  motions  the 
primary  extrinsic  contributor  to  both  dielectric  and  piezoelectric  response. 

No  practical  PZTs  are  in  fact  pure  lead  zirconate  titanate  solid  solutions,  all  are 
modified  by  aleovalent  cation  additions  (fig.  2).  We  speak  of  excess  charge  cation  added 
or  donor  doped  soft  PZTs,  and  of  charge  deficient  cation  or  Acceptor  doped  hard  PZTs. 
The  terminology  “hard"  versus  “soft”  is  derived  from  magnetism  and  clearly  implying 
difficult  or  easy  domain  wall  motion  as  in  the  magnetic  counterparts.  For  the  hard  PZTs 
however,  there  is  compelling  evidence  [4,5],  that  stabilization  by  weakly  mobile  oxygen 
vacancy:defect  dipole  pairs  is  of  the  whole  domain  not  just  the  wall.  A  mechanism  that 
has  no  counterpart  in  magnetism  but  does  however  indirectly  make  the  walls  more 
difficult  to  move. 

In  soft  donor  doped  PZTs,  there  is  very  strong  evidence  for  enhanced  extrinsic 
response  in  both  dielectric  and  piezoelectric  behavior,  however  evidence  that  it  is 
associated  with  domain  wall  motion  is  indirect.  It  is  clear  that  at  fields  well  below  the 
coercivity,  the  response  is  hysteritic  [6]  and  can  be  well  described  by  the  Rayleigh  law 
[7]  suggesting  a  domain  wall  origin.  There  is  however  also  evidence  from  transmission 
electron  microscopy  that  the  donor  dopants  tend  to  break-up  the  normal  domain  structure 
leading  to  tweed  like  and  micro-polar  relaxor  like  structures  [8]. 

Earlier  a  rather  complete  phenomenology  has  been  presented  for  the  whole  PZT 
family  [9-13]  from  which  single  domain  intrinsic  values  can  be  denied  for  both 
tetragonal  and  rhombohedral  phases  adjacent  to  the  MPB  and  averaged  values  for 
permittivity  and  piezoelectric  response  derived,  it  was  on  this  basis  that  some  80%  of 
dielectric  response  in  hard  PZTs,  but  only  of  order  30%  of  response  in  soft  PZTs  was 
predicted  to  be  intrinsic  (fig.  3). 

Very  recently  the  whole  basis  for  evaluating  both  contribution  to  response  has  been 
changed  by  the  discovery  using  very  precise  x-ray  analysis  [14,15]  of  a  new  monoclinic 
phase  just  on  the  rhombohedral  side  of  the  MPB  (fig.  4).  In  this  phase,  the  axial 
components  of  spontaneous  polarization  Ps  are 

Pj  =  P22^  O'  pj  &  0  but  pj  &  pj  (4) 
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DONOR  ADDITIVES 
Nb205  PbNb206 

Ta205  PbTa206 

W03 
Bi203 
Sb203 
La203 


Wall  Motion  Promoters 

d33 

593  pC/N 

d3i 

273  pC/N 

Tc 

195‘C 

e3 

3.450 

Kp 

0.69 

Tan  8 

0.02 

Ec 

6  Kv/cm 

(b) 

ACCEPTOR  ADDITIVES 
Fe203 
A1203 
Cr203 
Mn02 
MgO 
NiO 

(c) 

Wall  Motion  Inhibitors 
1 50  pC/N 
60  pC/N 
350'C 
1,000 
0.52 
0.004 
15  Kv/cm 


FIGURE  2.  (a)  Typical  polarization  and  strain  curves,  (b)  Donor  and  Acceptor  Dopants, 
(c)  Dielectric  and  Piezoelectric  Parameters 
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f  KHz 


The  polar  vector  lies  within  the  dodecahedral  minror  plane  of  the  prototype  m3m  and 
there  are  24  ^orientation  states  for  this  species.  From  symmetry,  the  transitions 
rhombohedral  ->m  <-*  monoclinic  m  and  tetragonal  mm  2  <-»  monoclinic  m  may  be  either 
first  or  second  order. 

The  presence  of  this  newly  discovered  phase  does  require  a  revision  of  the 
phenomenology,  as  in  the  earlier  treatment  all  phases  other  than  tetragonal  and 
rhombohedral  were  chosen  to  be  metastable  in  the  composition  ranae  of  the  MPB. 
Similarly  the  possible  extrinsic  contributions  to  response  now  need  further  consideration. 
The  possibility  for  the  E-field  to  drive  phase  boundary  motion  of  a  weakly  hvsteritic  first 
order  3m  h  m  or  4mm  m  phase  change  could  contribute  here-to-fore  unimagined 
extrinsic  response  that  could  mimic  domain  wall  motion  but  contribute  very  stronaly  to 
elastic  strain  and  piezoelectricity.  Either  intrinsic  or  extrinsic  contributions  from  this  very 
narrow  insert  of  monoclinic  phase  could  explain  the  very  sharp  peaking  of  k33  near  the 
MPB  that  is  evident  in  pure  PZTs  and  has  been  most  difficult  to  explain  on  earlier 
theories. 


Alternative  Single  Crystal  Systems 

The  complexity  of  the  elasto-dielectric  behavior  of  the  polycrystal  PZTs  highlight 
the  need  for  single  crystal  and  single  domain  studies  to  properly  explore  the  phenomena 
in  lead  based  perovskite  structures.  PZT  has  proven  markedly  intractable  for  crystal 
growth,  but  more  recent  studies  of  the  lead  zinc  niobaterlead  titanate  PZNrPT,  and  lead 
magnesium  niobaterlead  titanate  (PMN:PT)  have  yielded  “respectable”  crystals  for 
compositions  close  to  the  pseudo-morphotropic  boundaries  between  rhombohedral  and 
tetragonal  phases  in  both  systems  [16,17,18], 

Dielectric  and  piezoelectric  properties  have  been  well  characterized  [19,20]  and  the 
massive  anisotropy  in  response  between  111  and  001  E-field  poled  crystals  in  the 
rhombohedral  phase  underscored.  The  piezoelectric  d33  >  2000  pC/N,  k33  >  92%, 
e33  >  5,000  for  001  E-field  rhombohedral  phase  crystal  at  composition  close  to  the  MPB 
are  of  strong  practical  interest.  Of  major  concern  here  is  the  actuator  performance  where 
it  has  been  shown  that  in  001  field  poled  PZNrPT  strains  up  to  1.7%  can  be  induced  at 
high  fields  (fig.  5a). 

There  is  good  evidence  that  at  the  highest  field  levels  the  crystal  is  switched  over 
into  the  tetragonal  phase  through  a  field  induced  first  order  phase  change  that  is 
necessarily  hysteritic.  Of  special  interest  is  the  lojjig  quasi-Iinear  anhysteritic  strain,  up  to 
0.6 /o  (fig.  5b)  and  polarization  change  10  Uc/cm~  (fig.  5c)  that  can  be  induced  by  fields 
up  to  60  Kv/cm. 

It  is  suggested  that  the  sequence  of  phase  changes  induced  by  the  001  high  E-field 
is  as  in  fig.  6.  In  the  virgin  state,  the  crystal  is  a  relaxor  with  presumable  a  largely  frozen 
rhombohedral  nano/micro  domain  structure  (fig.  6a).  On  first  poling  with  001  oriented 
field_  the  _Pr  which  is  almost  exactly  1/Vi  of  the  rhombohedral  states,  i.e. 
111.  Ill  111,  and  111  oriented  polarizations  have  been  induced  (fig.  6b).  Note  that  this 
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(a)  E',  field  induced  longitudinal 
strain  in  001  Poled  PZN~92:PTS. 


(b)  Electric  Field  induced 
polarization  under  001  field  in 
poled  rhombohedrai  PZN  95.5  PT 
4.5. 


(c)  Electric  Field  Induced  Strain 
in  001  field  Poled  PZN  95.5  PT 
4.5  anhysteritic  for  repeated  field 
application. 


FIGURE  5.  Actuator  behavior  of  00 1  electric  field  poled  rhombohedrai  PZN :PT  single  crystal. 

structure  requires  a  high  concentration  of  charged  walls,  and  gives  rise  to  a  domain 
averaged  tetragonal  4mm  macro  symmetry.  Since  the  high  001  field  will  change  the 
energy  of  each  of  the  4  symmetry  equivalent  domains  in  the  same  way,  there  is  no 
driving  force  to  move  the  ferroelastic  walls  between  the  remaining  domain  states. 
Obviously  the  high  field  will  induce  a  tilting  of  the  Ps  vector  of  the  domain  into  the 
appropriate  dodecahedral  mirror  plane  givng  rise  to  monoclinic  domain  states  (fig.  6c)  in 
which 

P]  =  P22*0  Pj^O  P\>P2i  (5) 


Again  the  monoclinic  domains  will  average  to  macro  4mm  symmetry.  Eventually 
one  expects  the  field  to  induce  a  true  tetragonal  monodomain  state  in  which  macro  and 
micro  symmetries  coincide  (fig.  6d). 

An  interesting  question  is  whether  the  proposed  sequence  of  changes  can  account 
for  the  observed  large  strain  changes  without  extrinsic  domain  wall  contribution. 

Since  intrinsic  shape  changes  are  electrostrictive  it  is  only  necessary  to  know  the 
total  polarization  and  the  electrostrictive  coefficients  Qi  \,  Qi2>  2nd  Q44.  For  m3m,  it  is 
simple  to  show  that  the  volume  strain  sv  in  ail  possible  ferroelectric  phases  is  given  by 

s*-(Qu+*Qn)rb»*  ® 
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(a)  Virgin  Crystal  Cubic 
m3m  macro-symmetry 
rhombohedral  micro  domains 
all  orientation  states 


c)  Region  of  Linear  Strain 

Tetragonal  macro-symmetry 
4mm  monoclinic  domains 
m  monoclinic  polydomain 


b)  Poled  Crystal  Tetragonal 
4mm  macro-symmetry 
rhombohedral  domains 


d)  Switched  Tetragonal 
monodomain  at  very 
high  field  levels 


FIGURE  6. 

Sequence  of  states  in  high  001  field  poled  rhombohedral  PZN  95.5  PT  4.5  under  increasing  001  E-field. 


Measurements  by  Park  [20]  show  that  over  the  linear  strain  range  up  to  35  Kv/cm 
the  volume  change  in  PZN  0.95.5  PT  0.045  is  less  than  0.05%.  Taking  as  a  first 
aproximation  sv  =  0,  the  polarization  vector  is  just  rotating  under  field  and  from  fig.  5c 
the  tilt  angle  in  the  induced  monoclinic  phase  can  be  deduced  and  is  presented  in  Table  2. 

Taking  Qj  \  from  the  total  strain  and  total  polarization  in  the  induced  tetragonal 
monodomain  state 

Q]1  =  0.0535 

taking  Qh  =  0  then  yields 

^Q12  =  0.0267. 

In  the  field  induced  monoclinic  phase 

s!=(Q,1^Q;2)p]  +  oi2p] 

s  2  =  (Q  u +  Q  n)P2i  +  Q 12P3  (7) 

S3  =  (Q;i  +  p])  +  2Q12P] 
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TABLE  2.  Tilt  angle  and  polarization  components  under  increasing  001  E3  field  in  rhombohedral 
PZN;95.5  PT  4.5  single  crystal.  Tilt  is  measured  from  the  001  axial  direction. 


Field 

kV/cm 

Polarization  P3 
p.c/cm" 

AngleS 

Polarization  P\ 

ue/etrf 

5 

27.S 

50.1 

23.47 

10 

29.8 

46.5 

22.21 

15 

31.4 

43.4 

21.04 

20 

33 

40.34 

19.82 

25 

34.3 

37.5 

18.64 

30 

35.6 

34.6 

17.42 

35 

36.75 

31.9 

16. 19 

40 

37.85 

29.1 

14.87 

From  equation  7  and  the  known  Pj  and  P3,  components  of  the  strains  s^,  S2?  S3  can 
be  computed  for  the  induced  monoelinie  phase  in  the  anhysteritic  range  up  to  35  Kv/cm 
E3  field  (fig.  7).  Clearly  the  intrinsic  strains  in  the  monoclinic  form  are  adequate  to 
describe  the  induced  total  deformation. 


Longitudinal 


Transverse 
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40  60 

Electric  Field  (kV/cm) 


FIGURE  7.  Calculated  intrinsic  strain  in  the  field  induced  monoclinic  phase  (a)  S3  (b)  $\ 


If  this  proposed  explanation  for  high  strain  behavior  is  valid  it  should  be  expected 
that  other  peroskite  structure  ferroelectrics  would  show  high  intrinsic  strain  effects  at 
temperatures  close  to  phase  boundaries.  Strong  enhancement  of  <133  at  temperatures 
below  but  close  to  both  rhombohedral/orthorhombic  and  orthorhomic/tetragonal  phase 
changes  in  BaTi'03  single  crystals  have  been  demonstrated  [2],  Of  particular  relevance  is 
the  study  by  Wada  et.  al.  [3]  of  the  inverse  transition,  i.e.  tetragonal  to  rhombohedral  in 
111  cut  BaTi03  crystal  under  increasing  1 1 1  oriented  electric  field.  The  sequence  of 
changes  taken  from  their  paper  is  depicted  schematically  in  fig.  8.  The  1 1 1  poling  field 
first  induces  100,  010,  and  001  orientation  of  Ps  in  the  tetragonal  state.  Under  this 
condition  for  fields  below  5  Kv/cm  at  0.1  Hz  a  monoclinic  phase  is  induced  with  intrinsic 
d33  ~  203  pC/N  (fig.  8a). 

At  higher  fields,  the  system  becomes  hysteritic  and  the  family  of  1  JO,  101,  and  01 1 
orthorombic  domain  orientations  is  finally  induced  (fig.  8b). 

A  new  linear  regimen  of  strain  is  evident  for  fields  up  to  45  Kv/cm,  again  exhibiting 
monoclinic  symmetry,  until  the  system  moves  by  a  hysteritic  first  order  change  to 
rhombohedral  symmetry  (fig.  8c).  It  is  interesting  to  note  that  the  second  linear  strain 
range  has  a  different  slope  corresponding  to  a  d33  -  295  pC/N. 

We  suggest  (fig.  8)  that  because  of  the  closeness  to  room  temperature  of  the 
4mm/mm2  change  in  BaTiC>3,  the  polarization  vector  under  1 1 1  field  first  moves  out  into 
the  cube  mirror  (mj)  catalyzing  the  change  to  mm2  symmetry,  but  then  at  higher  field 
must  move  in  the  dodecahedral  mirror  (m2)  to  achieve  the  rhombohedral  orientation. 
After  achieving  the  rhombohedral  state,  it  would  appear  the  Ps  chooses  to  stay  in  the 
dodecahedral  mirror  which  connects  directly  to  the  tetragonal  state.  Since  the  two 
monoclinic  phases  accessible  from  m3m  are  quite  distinct  this  would  account  for  the 
different  piezoelectric  slopes  measured  in  the  experiment. 

Since  there  is  a  very  complete  Gibbs  Free  Energy  Function  for  BaTiC>3  [21]  it 
should  be  a  simple  matter  to  do  the  full  three  dimensional  energy  plots,  to  calculate  the 
lowest  energy  trajectories  for  the  polarization  change  and  thus  to  verify  the  major  role  of 
the  E*P  term  in  Energy  for  these  “soft”  ferroelectric  perovskites. 
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f>  XO  20  30  40  ~~SQ 

E  (kV/cxn) 

(c)  Strain  vs  electric-field  curve  for  [111]  oriented  BaTi03 
crystal  under  unipolar  electric  field  below  45  kV/cm 
withO.1  Hz  at  25*C. 

FIGURE  8.  Sequence  of  domain  and  strain  changes  in  1 1 1  poled  BaTiOs  crystal 
at  room  temperature  after  Wada  et.  al.  (2). 


Single  domain. 
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P32  i-  0  Pi2  =  P22  =  0  Tetragonal 

p32  jt  0  P22  ^  0  P32  t  P22  Pi2  s  0  monoclinic 

(cube  m) 


P32  =  P22  5s  0  Pi2  =  0 
P32  “  P22  T  0  Pi2  #  0 


P32  =  P22  =  Pi2  1 0 


Orthorhombic 

p32  monoclinic 

(Dodecahedron  m) 

Rhom  bohedral 


111 


FIGURE  9.  Suggested  model  for  the  inverse  field  induced  transition  from 
4mm  =»  3m  in  single  crystal  BaTiO;  at  25'C. 


Conclusions 

In  polycrystal  ceramic  PZT  piezoelectrics  the  current  balance  of  evidence  still 
suggests  a  strong  contribution  to  dielectric  and  piezoelectric  responses  from 
ferroelastic:ferroelectric  domain  wall  motion.  New  evidence  of  a  stable  monoclinic  phase 
just  on  the  rhombohedral  side  of  the  MPB  does  however  require  revision  of  the  earlier 
thermodynamic  phenomenology  to  describe  the  intrinsic  single  domain  responses  and 
also  contributes  the  possibility  of  previously  unimagined  new  extrinsic  contributions  to 
response. 

For  the  single  crystal  lead  zinc  niobate  lead  titanate  at  the  4.5%  lead  titanate 
composition  the  very  high  linear  anhysteritic  strain  induced  by  001  oriented  E-field  is 
susaested  to  be  intrinsic  and  to  result  from  the  induction  of  a  monoclinic  phase  in  which 
PsTs  substantially  tilted  from  111.  In  BaTi03  for  the  “inverse  experiment”  forcing 
tetragonal  to  rhombohedral  symmetry  at  room  temperature  it  is  suggested  that  with 
increasing  E-field  Ps  moves  first  in  the  cube  mirror  to  orthorhombic,  then  in  the 
doecahedral  mirror  towards  rhombohedral  symmetry,  but  at  very  high  field  may  chose  to 
accomplish  complete  switching  in  the  alternative  dodecahedral  mirror. 
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Many  types  of  symmetry  are  utilized  in  ceramic  art  and  in  ceramic  science.  In 
addition  to  the  common  mirror,  rotation,  inversion  and  translation  symmetry 
elements  observed  in  crystals  and  textured  ceramics,  artists  and  engineers  often 
make  use  of  the  antisymmetry  elements  involved  in  color  groups.  We  illustrate 
the  application  of  color  symmetry  in  controlling  the  vibration  modes  in 
ferroelectric  and  ferromagnetic  ceramics  prepared  as  spheres,  hemispheres, 
and  rings  conforming  to  black-and-white  Curie  groups. 

SYMMETRY  IN  ART  AND  SCIENCE 

While  visiting  the  city  of  Kutahya,  one  of  the  great  centers  for  traditional 
ceramic  tiles  and  porcelain  products,  we  had  the  opportunity  to  view  contem¬ 
porary  Turkish  artwork.  Islamic  art  has  long  been  famous  for  its  geometric 
patterns,  but  we  were  especially  fascinated  with  the  porcelains  made  by  Kervan 
-  Chini,  who  incorporates  many  interesting  symmetries  into  his  pottery.  The 
colorful  dish  shown  in  Fig.  1  contains  a  number  of  the  symmetries  and  broken 
symmetries  observed  in  electroceramic  materials  and  devices:  polar  symmetry, 
chiral  symmetry,  symmetry  of  scale,  quasisymmetry,  the  pseudosymmetry  of 
incommensurate  phases,  and  the  antisymmetric  elements  of  color  symmetry 
groups. 

Piezoelectric  ceramics,  pyroelectric  glass  ceramics,  and  certain  function¬ 
ally-graded  ceramics  have  polar  symmetry  like  that  of  a  concave-shaped  disk 
(Fig.  2a).  Chiral  symmetry  is  present  in  the  porcelain  dish  as  well  (Fig.  2b).  The 
leaf  pattern  of  the  flowers  and  the  sleeping  sheep  near  the  edge  are  both 
arranged  in  a  counterclockwise  fashion  leading  to  handedness.  Chiral  phenom¬ 
ena  such  as  optical  activity,  acoustic  activity,  and  the  Faraday  Effect  are  well 
known  in  crystal  physics.  Symmetry  of  scale  is  important  in  all  classes  of 
materials.  The  pattern  in  Fig.  1  violates  the  self-similarity  principle  introduced 
by  Benoit  Mandelbrot  to  describe  fractal  geometry.  Beginning  at  the  center, 
the  porcelain  dish  shows  5-fold  then  9-fold,  10-fold,  and  19-foid  rotational 
symmetry.  Five-fold  symmetry  is  characteristic  of  quasicrystals  and  the  circle 
possesses  a  oc-fold  axis,  one  of  Curie  group  symmetries  found  in  textured  poly¬ 
crystalline  ceramics.  The  change  in  symmetry  with  scale  is  illustrated  in  Fig.  2c. 
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The  symmetry  of  scale  is  an  important  one  in  many  electronic  systems 
as  component  sizes  become  smaller  and  smaller.  Current  trends  in  multilayer 
capacitors  are  shown  in  Fig.  3  and  4.  The  lateral  dimensions  in  MLC  chips  are 
now  about  1  mm,  approaching  the  limit  for  many  of  the  pick-and-place 
machines  used  in  assembling  circuits.  Layer  thicknesses  and  grain  sizes  are 
also  a  concern.  Ceramists  are  now  testing  BaTi03  MLCs  with  layers  of  1  pm 
thickness  which  is  pressing  the  limits  of  tape-casting  technology.  Particle  sizes 
are  in  the  0.1  pm  range  which  raises  some  fundamental  questions  regarding 
size  effects  in  ferroelectrics.  In  large  grain  size,  say  1  pm  and  larger,  each  grain 
contains  many  domains.  The  number  of  domains  and  the  type  of  domain 
walls  change  as  grain  size  drops  below  a  micron.  Eventually  each  grain  becomes 
a  single  domain  with  profound  changes  in  dielectric  constant  and  switching 
behavior.  In  the  nm  range  the  ferroelectric  phase  transition  becomes  diffuse 
and  the  symmetry  of  BaTi03  appears  pseudocubic.  Size  effects  in  primary 
ferroics  are  illustrated  schematically  in  Fig.  5.  Dielectric  constants  drop 
substantially  for  grain  sizes  less  than  0.1  pm. 

Returning  to  the  porcelain  plate  in  Fig.  1,  the  outer  portions  of  plate 
provide  an  excellent  illustration  of  pseudosymmetry  or  "almost"  symmetry.  In 
the  central  leaf  pattern  there  are  10  leaves  and  between  neighboring  leaves 
are  1 0  flowers.  Near  the  outside  rim  of  the  plate  are  red  hearts  which,  at  first 
glance,  appear  to  be  in  register  with  the  1 0  leaves  and  with  the  1 0  flowers,  but 
such  is  not  in  the  case.  The  artist  has  made  a  subtle  change  in  symmetry  with 
only  19  hearts  instead  of  20  (Fig.  2d).  Many  ferroelectric  and  ferrimagnetic 
oxides  possess  incommensurate  structures  in  which  the  local  polarization  or 
magnetization  vectors  are  out  of  register  with  the  lattice  periodicity  (Fig.  6). 
These  so-called  incommensurate  phases  often  exhibit  unusual  physical 
properties  because  of  their  abnormal  symmetry. 

FERROIC  CRYSTALS  AND  ANTISYMMETRY  ELEMENTS 

The  symmetry  elements  discussed  thus  far  are  purely  spatial  trans¬ 
formations  such  as  mirror  planes,  rotation  axes,  and  inversion  centers.  Spatial 
symmetry  elements  are  all  that  are  required  for  the  usual  crystallographic  and 
limiting  point  groups,  but  additional  symmetries  occur  in  ferroic  electro¬ 
ceramics.  These  additional  symmetries  can  be  described  by  color  symmetry. 
The  ceramic  dish  from  Turkey  possesses  decagonal  color  symmetry  in  which 
the  flower  rotates  by  36°  and  changes  color  (Fig.  2e).  In  this  paper  we  discuss 
color  groups  and  the  antisymmetry  elements  found  in  ferroic  ceramics  with 
complex  domain  patterns.  The  10'  antisymmetry  element  in  Fig.  2e  belongs 
to  this  type  of  color  group. 

More  than  a  century  ago,  the  basic  relationships  between  symmetry 
and  physical  properties  were  established  by  Neumann  and  the  brothers  Gu- 
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rie,  and  systematized  in  Woldemar  Voigt's  monumental  "Lehrbuch  der 
Kristallphysik."  (Voigt,1928)  Using  Neumann's  Principle,  the  symmetry 
restrictions  for  piezoelectricity,  magnetostriction,  and  other  linear  and  nonlinear 
tensor  properties  have  been  enumerated,  and  have  led  to  the  development  of 
many  types  of  sensors,  actuators,  and  transducers.  For  single  crystals,  the 
properties  are  governed  by  the  32  crystallographic  point  groups,  and  for 
textured  polycrystalline  materials,  we  use  the  seven  limiting  groups  first 
described  by  Pierre  Curie.  In  1 974  we  extended  these  symmetry  arguments  to 
ferroic  crystals,  laying  out  the  symmetry  changes  involved  in  the  phase 
transformations  leading  to  the  hysteretic  phenomena  created  by  domain  wall 
motion  (Newnham,  1974).  There  are,  for  example,  15  symmetry  changes 
consistent  with  pure  ferroelastoelectric  behavior  in  which  adjacent  domains 
differ  only  in  the  orientation  of  third  rank  polar  tensor  properties.  A  "pure" 
ferroelastoelectric  does  not  exhibit  any  other  primary  or  secondary  ferroic 
behavior  (Newnham  and  Cross,  1974}. 

In  crystals  with  long  range  magnetic  order,  the  addition  of  a  time  reversal 
operator  leads  to  the  generation  of  the  90  magnetic  point  groups  used  to 
describe  ferromagnetic,  ferrimagnetic,  and  antiferromagnetic  substances  (Birss, 
1 964).  Introduction  of  the  antisymmetric  spin  reversal  operator  makes  it  possible 
to  apply  Neumann's  Principle  to  magnetic  crystals  and  classify  the  symmetry 
restrictions  for  magnetic  properties  such  as  magnetoelectricity,  pyromagrietism 
and  piezomagnetism.  All  three  of  these  cross-coupled  phenomena  are  axial 
tensor  properties  and  since  spin  reversal  corresponds  to  the  reversal  of  an 
axial  vector,  all  three  properties  are  strongly  influenced  by  the  antisymmetry 
operator. 

In  this  paper  we  introduce  antisymmetry  elements  for  other  types  of 
ferroic  materials  with  movable  domain  walls.  Ferroelectricity,  ferroelasticity, 
and  the  six  types  of  secondary  ferroics  are  also  controlled  by  tensorial 
antisymmetry  operators  (Table  1).  For  a  ferroelectric,  polarization  -  a  polar 
vector  -  is  the  key  operator,  since  ferroelectricity  is  defined  by  polarization 
reorientation  between  symmetry  related  states.  For  a  ferroelastic,  the  key 
operator  is  strain  -  a  second  rank  polar  tensor  -  since  ferroelasticity  is  defined 
by  strain  reorientation  between  symmetry  related  states.  For  a  ferrobielastic 
secondary  ferroic,  the  key  operator  is  a  fourth  rank  tensor,  etc. 

The  procedure  will  adopt  for  introducing  these  antisymmetry  elements 
is  as  follows.  Each  operator  is  combined  with  the  spatial  symmetry  elements 
of  the  32  crystallographic  point  groups  or  to  the  seven  Curie  groups  to  gene¬ 
rate  a  family  of  antisymmetric  groups  analogous  to  the  90  magnetic  point 
groups.  Having  determined  the  relevant  symmetries,  Neumann's  Principle  is 
then  applied  to  determine  the  polar  and  axial  property  matrices,  enabling  one 
to  predict  the  useful  symmetries  that  will  optimize  the  material  for  a  given 
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FIGURE  7  -  A  ceramic  dish  made  by  Kervan  Chini  which  illustrates 
many  types  of  unusual  symmetry. 
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FIGURE  2a  -  Polarity. 
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FIGURE  2b  -  Handedness. 
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FIGURE  2c  -  Size  effects. 


FIGURE  2d  -  Aperiodicity. 
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FIGURE  2e  -  Color  patterns. 
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1980  1985  1990  1995  2000 

YEARS 

FIGURE  3  -  The  sizes  of  ceramic  capacitors  have  decreased  steadily  with  time  while 
the  market  continues  to  grow  exponentially  following  Moore's  Law. 
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FIGURE  4  -  Layer  thicknesses  in  multilayer  capacitors  are  now  nearing  one  micron  with 
grain  sizes  around  0.2  j.i m. 
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FIGURE  5  -  Profound  changes  in  properties  take 
place  in  ferroic  ceramics  when  the  grain  sizes 
■drop  below  the  micron  range. 


FIGURE  6  -  Formation  of  incommensurate  phase 
as  a  result  of  the  "freezing"  of  a  displacement 
wave  with  a  length  incommensurate  with  the 
unit-cell  parameter. 
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TABLE  7  -  Nine  types  of  ferroics  in  which  domain  walls  can  be  moved  by  magnetic 
fields ,  electric  fields ,  or  mechanical  stress.  Each  is  controlled  by  a  different  term  in  the 
free  energy  function  with  an  antisymmetry  operator  relating  domain  state  orientations. 


Type 

Domain  Difference 

Antisymmetry  Tensor 

Primarv  Ferroics 

Ferromagnetic 

Magnetization 

Axial  First  Rank 

Ferroelectric 

Polarization 

Polar  First  Rank 

Ferroelastic 

Strain 

Polar  Second  Rank 

Secondary  Ferroics 

Ferrobielectric 

Permittivity 

Polar  Second  Rank 

Ferrobi  magnetic 

Permeability 

Polar  Second  Rank 

Ferrobielastic 

Elasticity 

Polar  Fourth  Rank 

Ferroelastoelectric 

Piezoelectricity 

Polar  Third  Rank 

Ferromagnetoelectric 

Magnetoelectricity 

Axial  Second  Rank 

Ferromagnetoelastic 

Piezomagnetism 

Axial  Third  Rank 

application. 

Knowing  the  optimum  symmetry  groups,  how  does  one  go  about  creating 
these  symmetries  in  a  real  material?  Here  one  applies  Curie's  Principle  of 
Symmetry  Superposition  and  examines  various  composite  ferroics  (Newnham 
and  Ciniewicz,  2000).  Symmetry  patterns  can  be  generated  at  many  different 
size  scales  ranging  from  electron  spin  arrangements,  crystal  structures,  defect 
structures,  textured  grains,  tailored  domain  patterns,  electrode  geometries, 
connectivity  patterns,  and  various  morphologies  using  processing  methods  such 
as  extrusion,  tape-casting,  and  photolithography.  Many  of  the  more  unusual 
point  groups  can  be  accessed  when  ferroic  materials  are  prepared  as  fibers, 
films,  tubes,  disks,  multilayers  and  hollow  spheres.  Using  Curie's  Principle  of 
Symmetry  Superposition,  many  unusual  cross-coupled  properties  can  be 
optimized  by  symmetry  control.  The  new  feature  we  are  developing  is  a  family 
of  antisymmetry  operators  (polar  vector  reversal,  polar  second  rank  tensor 
reversal,  axial  second  rank  tensor  reversal,  etc.)  which  Will  augment  our 
understanding  of  the  symmetry  of  solids,  and  enable  us  to  predict  new  families 
of  functional  materials  for  biomedical  and  industrial  applications.  Our  goal  is 
to  control  the  acoustic  and  electromagnetic  resonance  phenomena  in  a  wide 
range  of  electroceramic  products. 
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ANTISYMMETRIC  CURIE  GROUPS 

In  this  section  we  describe  a  family  of  small  electroceramic  transducers, 
sensors,  and  actuators  based  on  black-and-white  Curie  group  symmetries  (Fig. 
7).  One  of  the  goals  of  this  investigation  is  to  control  the  acoustic  vibration 
modes  and  radiation  patterns  by  controlling  the  symmetries  of  the  transducer's 
material,  its  external  shape  and  poling  pattern,  together  with  the  driving 
electrode  geometry.  Later  we  intend  to  extend  these  symmetry  arguments  to 
electromagnetic  wave  modes  as  in  millimeter  wave  resonators  and  dielectric 
antennae. 

We  begin  with  an  example  to  illustrate  the  basic  idea.  Consider  a 
spherical  shell  made  of  a  piezoelectric  ceramic  such  as  PZT  with  randomly 
oriented  grains.  When  statistically  averaged,  randomly  oriented  grains  in  an 
unpoled  ceramic  conform  to  spherical  symmetry  ooooml '  where  oooo  represents 
an  infinite  number  of  infinite-fold  rotation  axes,  m  is  an  infinite  number  of 
mirror  planes,  and  1 '  is  the  color  reversal  operator  changing  white  to  black, 
and  vice  versa. 

In  ferroelectric  substances  such  as  PZT  we  can  use  color  to  represent 
polar  vector  reversal.  When  electrically  poled  in  a  strong  DC  field  near  T  the 
altered  domain  structure  changes  the  symmetry  to  oo/m'm,  the  symmetry  of 
an  electric  field.  This  is  also  the  overall  symmetry  of  an  electric  dipole  with  its 
positive  and  negatively-charged  ends.  The  oo-fold  axis  is  parallel  to  the  dipole. 


FIGURE  7  -  The  seven  Curie  groups  (solid  boxes)  and  fourteen  black  and  white  Curie 
groups  (dashed  boxes).  The  interconnecting  lines  denote  group-subgroup  relationships. 
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m'  is  a  mirror  plane  accompanied  by  charge  (color)  reversal  perpendicular  to 
the  oo-fold  axis.  In  the  symbol  oo/m'm,  m  refers  to  the  infinite  number  of  mirror 
planes  parallel  to  the  oo-fold  dipole  axis.  This  is  the  symmetry  of  the  material, 
and  that  in  turn  determines  the  physical  properties  through  Neumann's  Law 
from  crystal  physics.  That  of  course  is  why  poled  PZT  is  piezoelectric  and 
pyroelectric. 

But  that  is  not  the  only  symmetry  that  is  important  in  transducer  design. 
The  vibration  modes  and  radiation  patterns  depend  on  the  symmetry  of  the 
transducer's  shape,  poling  pattern,  and  driving  electrodes  as  well.  If  we  consider 
a  spherical  shell  as  an  example,  the  shape  conforms  to  spherical  symmetry 
coooml ',  and  if  it  is  poled  radially,  the  symmetry  group  changes  to  oooom.  This  is 
a  black-  and-  white  symmetry  group  which  exists  in  two  physically-distinct 
states  corresponding  to  a  hollow  sphere  poled  inside-out  or  outside-in.  The 
principal  piezoelectric  and  pyroelectric  coefficients  differ  in  sign  for  the  two 
states  which  are  related  to  one  another  by  the  symmetry  element  m<t,  reflection 
accompanied  by  dipole  reversal. 

The  symmetry  of  the  radially-poled  sphere  (oooom)  together  with  the 
symmetry  of  the  material  (co/m'm)  determines  the  vibration  modes  of  the 
transducer.  The  material  symmetry  determines  the  non-zero  piezoelectric 
coefficients:  d33,  d31=d32,  d24=d15.  For  a  radially-poled  sphere,  the  poling 
direction  X3  and  the  driving  field  E3  are  in  the  radial  direction  r,  while  Xv  and 
X2  are  orthogonal  tangential  directions  designated  by  the  symbol  t.  Therefore 
d33  =  dif  and  d31  =  d32  =  dft.  When  driven  electrically,  the  radially-poled  sphere 
vibrates  in  two  fundamental  modes:  a  high-frequency  wall  thickness  mode 
governed  by  drr  and  a  low-frequency  breathing  mode  controlled  by  drt.  For  a 
hollow  sphere  2  mm  in  diameter  with  a  wall  thickness  of  about  0.1  mm,  these 
resonances  occur  near  20  MHz  and  600  kHz,  respectively  (Alkoy,  1999).  As 
discussed  later,  a  much  more  complex  mode  spectrum  with  ellipsoidal  modes 
(symmetry  group  oo/mm)  is  observed  when  the  sphere  is  electroded  at  the  N 
and  S  poles  and  poled  tangentially. 

Three  additional  types  of  symmetry  control  can  also  be  illustrated  with 
this  example.  The  first  involves  re-electroding.  The  driving  electrode  pattern 
need  not  be  the  same  as  the  poling  electrode  pattern.  Poling  is  generally  carried 
out  at  elevated  temperatures  under  high  voltage.  After  cooling  and  aging  the 
poling  pattern  is  relatively  stable.  The  poling  electrodes  can  then  be  removed 
and  a  new  set  of  driving  electrodes  installed  which  create  different  vibration 
modes.  This  is  indicated  schematically  in  Fig.  8.  In  this  way  a  radially  poled 
sphere  might  be  driven  tangentially. 

The  second  modification  involves  the  use  of  a  second  transducer  attached 
to  the  first  transducer  and  driven  in  a  different  mode.  It  is  relatively  easy  to 
build  transducers  in  multilayer  form  with  a  common  ground  electrode  and 
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FIGURE  8  -  Outline  of  experimental  approach. 


independently-driven  external  electrodes.  This  enables  us  to  superimpose  two 
vibration  modes  driven  with  different  phases  and  amplitudes. 

The  practical  aim  of  our  program  is  to  control  the  external  radiition 
pattern  by  controlling  the  symmetry.  If  one  transducer  is  driven  in  an  omni¬ 
directional  fundamental,  and  the  second  transducer  in  a  dipole  mode,  the  net 

result  is  a  cardioid  mode.  In  this  way  one  can  direct  the  radiation  in  the  desired 
direction. 

All. of  these  patterns  and  their  symmetries  can  be  controlled  indepen¬ 
dently  by  choosing  the  material  and  its  processing.  Our  current  objective  in 
Ibis  research  is  to  access  the  black-and-white  symmetries  in  Fig.7  to  generate 
interesting  acoustic  vibration  modes  using  piezoelectric  and  magnetostrictive 
transducers. 


FERROIC  SPHERES,  HEMISPHERES,  AND  RINGS 

Our  initial  aim  is  to  study  the  vibration  modes  and  radiation  patterns 
from  the  fourteen  black-and-white  Curie  groups  oooom,  oooom',  oooo,  oo/m'm,  00/ 
mm',  oo/m'm',  00/mm,  co/m,  co/m',  oo2,  oo2',  00m',  oom,  and  00.  All  fourteen 
symmetries  can  be  accessed  with  ferroic  hollow  spheres  or  modified  hollow 
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spheres  shaped  as  rings  or  hemispheres.  When  polished  down  on  one  side, 
hollow  spheres  (symmetry  group  ooooml')  become  hollow  hemispheres 
(symmetry  group  ooml  ')•  When  polished  symmetrically  on  two  opposite  sides, 
hollow  spheres  become  rings  (symmetry  group  oo/mml  ')• 

Following  the  processing  steps  laid  out  in  the  experimental  procedure 
section,  polycrystalline  hollow  spheres  are  prepared  from  strongly  piezoelectric 
and  magnetostrictive  substances  such  as  PbZr1xTix03  (PZT)  and  Tb1x  DyxFe2 
(Terfenol-D).  After  poling,  magnetizing,  and  polishing  all  14  of  the  black-and- 
white  Curie  groups  can  be  generated  according  to  the  game  plan  shown  in 
Table  2. 


TABLE  2  -  Fourteen  black-and-white  symmetry  groups  to  be  obtained  from  processed 
hollow  spheres.  Initially  we  are  concentrating  on  the  piezoelectric  elements ,  and  then 
the  magnetostrictive  transducers.  Combined  elements  are  to  be  examined  in  the  third 
stage  using  multilayer  processing. 


Antisymmetric  Curie  Group 

Ferroic  Morphology 

cooom 

Radially-poled  piezoelectric  sphere 

oooorn' 

Radially-magnetized  magnetostrictive  sphere 

COOO 

Concentric  radially-poled  piezoelectric  sphere  and 
radially-magnetized  magnetostrictive  sphere 

oo/mm 

Radially-poled  piezoelectric  ring 

oo/mm' 

Longitudinally-magnetized  magnetostrictive  ring 

oo/m'm 

Longitudinally-poled  piezoelectric  ring 

oo/rn'm' 

Radially-magnetized  magnetostrictive  ring 

oo/m 

Concentric  radially-poled  piezoelectric  and 
longitudinally-magnetized  magnetostrictive  rings 

.  oo/m' 

Concentric  longitudinally-poled  piezoelectric  and 
radially-magnetized  magnetostrictive  rings 

oo2 

Concentric  radially-poled  piezoelectric  and 
radially-magnetized  magnetostrictive  rings 

oo2' 

Concentric  longitudinally-poled  piezoelectric  ring 
and  longitudinally-magnetized  magnetostrictive 

ring 

com 

Radially-poled  piezoelectric  hemisphere 

oom' 

Radially-magnetized  magnetostrictive  hemisphere 

OO 

Concentric  radially-poled  piezoelectric  and 
radially-magnetized  magnetostrictive  hemisphere 
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The  next  step  is  to  determine  the  resonant  vibration  modes  experi¬ 
mentally  using  impedance  spectroscopy  and  to  verify  the  measurements 
theoretically  using  the  ATILA  finite  element  code.  Typical  calculated  spectra 
fora  radially-poled  piezoelectric  hollow  sphere  and  a  tangentially  poled  hollow 
sphere  are  illustrated  in  Figures  9  and  10.  For  the  radially-poled  sphere,  the 
low-frequency  mode  near  600  kHz  is  a  breathing  mode  driven  by  piezoelectric 
coefficient  d31,  and  the  high-frequency  20  MHz  wall-thickness  mode  is 
governed  by  d33.  Miniature  hydrophones  and  minipumps  make  use  of  the 
large  hydrostatic  piezoelectric  coefficient  dh  of  these  transducers.  Hollow  sphere 
sensors  and  actuators  have  a  stress  amplification  factor  proportional  to  the 
sphere  radius  divided  by  the  wall  thickness,  which  is  typically  more  than  an 
order  of  magnitude  in  our  transducers. 

Similar  amplification  factors  apply  to  the  tangentially  poled  sphere  in 
Fig.  10.  This  nicely  illustrates  how  to  control  the  vibration  modes  and  their 
symmetry  with  the  electrode  pattern.  In  this  case  different  size  circular 
electrodes  were  applied  to  the  top  and  bottom  of  the  sphere,  conforming  tq 
antisymmetry  group  oo/m'm.  When  driven  electrically,  a  series  of  ellipsoidal 
and  egg-shaped  vibrations  were  excited,  as  illustrated  in  Figure  10.  The  ATILA 
code  is  very  helpful  in  identifying  and  visualizing  the  various  modes  of  vibration 
of  complex  patterns  and  shapes. 


FIGURE  9  -  The  (a)  admittance  spectrum  of  a  hollow  sphere  electroded  inside  and 
outside  and  poled  radially.  The  principal  vibration  modes  are  the  (b)  low  frequency 
breathing  mode  and  the  (c)  high  frequency  wall  thickness  mode.  Both  modes  are 
spherically  symmetric  conforming  to  the  shape  and  electrode  pattern.  The  dashed  lines 
correspond  to  the  rest  position  and  the  solid  line  to  the  vibration  motion.  Scales  are 
arbitrary.  (Alkoy,  1999) 
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FIGURE  70  -  Comparison  of  the  (a)  calculated  and  measured  admittance  response  of 
(b)  a  tangentially-poled  hollow  sphere  with  asymmetric  electrodes.  The  displacement 
fields  of  four  vibration  modes ,  two  of  which  are  ellipsoidal  in  shape  (bo /mm'),  one  is 
spherical  (cocom'l  and  one  egg-shaped  (oom'J.  This  illustrates  Neumann's  Principle  which 
says  that  the  symmetry  of  a  property  must  include  the  symmetry  of  the  material  and 
the  way  in  which  it  is  driven  (Alkoy,  7  999). 


TRANSDUCER  FABRICATION 

The  fabrication  process  is  designed  to  produce  spherical  electroceramic 
devices  with  cofired  inner  electrodes.  For  commercial  applications  we  need  a 
process  capable  of  producing  large  quantities  of  hollow  sphere  transducers  in 
sizes  ranging  from  a  fraction  of  a  millimeter  to  several  millimeters.  The  process 
is  illustrated  schematically  in  Fig.  11. 
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a  slurry  containing  metal  powders.  The  electrode  is  then  dried  and  further 
serves  as  a  mold  for  the  ceramic  slurry.  A  ceramic  shell  is  then  deposited  on 
top  of  the  dried  electrode  by  dipping  into  a  well  dispersed  ceramic  slurry 
containing  binder  and  dispersant.  The  ceramic  slurry  can  be  aqueous  or  non- 
aqueous.  The  ceramic  layer  thickness  is  controlled  by  immersion  time,  the 
solid  loading  of  the  ceramic  slurry  and  its  viscosity.  When  a  large  wall  thickness 
is  desired,  multiple  coatings  can  be  performed.  When  the  sphere  is  removed 
from  the  slurry,  excess  slurry  on  the  sphere  can  be  removed  by  spinning  the 
sphere.  Spinning  also  improves  the  uniformity  of  the  wall  thickness.  The  porous 
polystyrene  spheres  and  the  electrode  layer  absorb  water  from  the  ceramic 
layer.  The  capillary  suction  allows  a  dense  green  body  of  ceramic  shell  to  be 
formed,  which  is  crucial  in  sintering  a  fully  densified  ceramic  body. 

After  forming,  the  spheres  are  dried  and  taken  through  a  polymer  burnout 
process  in  the  temperature  range  between  300  and  700  °C.  A  ceramic  shell 
with  an  inner  electrode  is  obtained  after  this  process.  Then  the  ceramic  shell  is 
cofired  with  its  electrode.  The  sintering  conditions  vary  according  to  the  ceramic 
composition.  In  most  cases,  the  ceramic  and  electrode  materials  are  sintered 
in  the  temperature  range  between  1200  and  1400  °C.  The  immersion  time 
and  optimum  sintering  temperature  must  be  adjusted  experimentally  to  achieve 
the  maximum  sintering  density  of  ceramic.  For  lead  containing  compositions, 
the  ceramic  shell  needs  to  be  fired  in  a  closed  container  in  a  lead  rich 
atmosphere  to  reduce  lead  loss  from  the  ceramic  body. 

After  sintering,  the  outer  surface  of  the  shell  is  coated  with  electrode 
metal  through  sputtering  or  dipping.  A  tiny  lead  wire  is  then  attached  to  the 
outer  electrode.  The  ceramic  shell  is  poled  in  the  case  of  piezoelectric  materials 
by  subjecting  the  ceramic  shell  to  high  electric  field  at  elevated  temperatures 
in  the  range  between  1 00  and  200  °C.  This  completes  the  process  of  fabricating 
a  hollow  spherical  transducer. 

Multilayer  hollow  sphere  transducers  (Fig.  12)  can  be  fabricated  by  a 
multiple  dipping  process  to  give  a  structure  with  alternating  electrode  and 
ferroic  ceramic  layers.  After  the  first  ceramic  layer  is  deposited,  it  is  dried  and 
then  coated  with  another  layer  of  electrode  material,  followed  by  a  coating  of 
another  layer  of  ceramic.  Care  must  be  taken  to  make  certain  the  second 
electrode  layer  does  not  touch  the  platinum  wire  and  a  hole  must  be  left  in 
the  second  ceramic  layer  for  electrical  connection.  The  goal  is  to  make 
multilayer  hollow  spheres  with  many  thin  layers  similar  to  multilayer  capacitors 
and  multilayer  actuators. 

After  polymer  burnout  and  sintering,  the  electrodes  are  connected  to 
electrical  lead  wires.  Different  connections  are  available  depending  on 
applications.  For  receive  applications,  the  layers  can  be  connected  in  series  to 
improve  the  receive  response.  For  projector  applications,  the  layers  can  be 
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FIGURE  12-  A  multilayer  hollow  spherical  transducer 
prepared  using  a  multiple  coating  process. 


connected  in  parallel  to  increase  the  pressure  output  to  transmit  directional 
beams  with  the  electrodes  driven  with  different  phases  and  amplitudes. 

To  access  certain  of  the  antisymmetric  Curie  group  symmetries  (Table  2) 
it  is  necessary  to  convert  the  hollow  spheres  into  hemispheres  or  rings.  This 
step  is  carried  out  by  mounting  the  hollow  sphere  in  a  polymer  matrix  and 
polishing  the  sphere  down  on  one  side  (for  the  hemisphere)  or  two  sides  (for 
die  ring).  The  polymer  matrix  provides  the  mechanical  strength  required  for 
the  lapping  operation. 

Limited  re-electroding  can  be  accomplished  after  removing  the  outer 
electrodes  used  in  poling  the  transducer.  This  enables  us  to  apply  driving  fields 
in  directions  other  than  the  poling  field  directions,  and  thereby  generate 
different  symmetries  and  different  vibration  modes. 

ACOUSTIC  APPLICATIONS 

Applications  are  not  the  main  focus  of  this  discussion  but  it  is  important 
to  have  practical  goals  in  mind  to  justify  the  research.  The  increasing  use  of 
ultrasonic  techniques  in  biomedical  and  underwater  applications  and  the  soph¬ 
istication  of  these  techniques  has  increased  the  need  to  characterize  the  acoustic 
fields  using  ultrasonic  transducers.  In  the  last  forty  years  several  different  tech¬ 
niques  and  probe  designs  have  been  employed  for  this  purpose,  including 
rod-guided  waves,  thermoacoustic  sensors,  fiberoptic  sensors  utilizing  acousto- 
optic  interaction,  and  miniature  acoustic  probes  made  from  ceramic  and 
polymer-based  piezoelectrics.  Miniature  ultrasonic  probes  have  been  used 
both  for  mapping  the  field  of  a  hydrophone  as  well  as  the  nonacoustic  field  of 
turbulent  flow.  There  are,  however,  several  important  requirements  for  micro¬ 
probe  sensors  in  these  applications.  In  detecting  underwater  signals  or  biological 
pressure  waves,  omnidirectionality  is  highly  advantageous.  Accurate  mapping 
of  an  acoustic  field  requires  that:  (i)  the  physical  dimensions  of  the  probe 
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should  be  smaller  than  the  acoustical  wavelength  of  Interest,  (ii)  the  resonance 
frequencies  of  the  probe  should  be  well  above  the  frequency  range  of  interest, 
(iii)  adequate  sensitivity  with  an  acceptable  signal-to-noise  ratio  and  (iv)  wide 
bandwidth. 

Although  volume  expanders  with  spherical  shape  are  thought  to  be  the 
best  way  to  achieve  omnidirectionality,  there  are  problems  associated  with 
fabricating  hollow  spheres  with  sizes  in  the  millimeter  range,  as  pointed  out 
by  previous  researchers.  There  are  a  number  of  examples  of  transducers 
prepared  from  piezoelectric  plates,  but  these  transducers  have  a  pronounced 
directivity  even  when  the  probe  dimensions  are  smaller  than  a  wavelength. 
Solutions  proposed  previously  for  the  directivity  problem  include  solid-core 
spherical  probes  and  hollow  cyclindrical  probes. 

In  addition  to  the  underwater  applications,  spherically  shaped  transducers 
are  also  desirable  in  medical  ultrasound  applications.  Examples  include  a  quasi- 
omnidirectional  polymer-based  transducer  developed  for  ultrasonic  guidance 
of  intravascular  catheters  and  a  focused  transducer  for  biomedical  ultrasonic 
imaging,  prepared  from  a  spherical  ceramic  shell.  Every  year  millions  of  medical 
procedures  are  performed  in  the  United  States  using  catheter  systems  a  few 
millimeters  in  diameter.  Typically,  they  are  inserted  through  the  skin  into  blood 
vessels  or  into  various  body  cavities  to  deliver  drugs  or  therapeutic  devices. 

At  present  catheters  are  guided  with  real-time  X-ray  imaging  (fluoroscopy) 
despite  the  advantages  of  ultrasound  imaging  over  X-ray  imaging  in  terms  of 
cost,  safety  and  availability.  Because  the  wavelengths  used  in  medical  ultrasonic 
imaging  (0.2-0.5  mm)  are  several  times  smaller  than  the  catheter,  it  acts  as  a 
reflector  and  the  ultrasonic  visibility  of  the  catheter  is  highly  direction  dependent 
with  regard  to  external  ultrasonic  beams.  One  solution  to  this  problem  is  to 
mount  an  omnidirectional  ultrasonic  transducer  on  the  catheter  that  can  ser¬ 
ve  as  an  ultrasonic  marker  to  help  locate  the  catheter  without  the  use  of 
hazardous  X-rays.  The  applications  of  these  ultrasonic  imaging  catheters  inclu¬ 
de  guiding  balloon  angioplasties  of  the  leg,  guiding  catheters  inside  the  heart 
to  ablate  incorrectly  functioning  cardiac  tissue,  and  guiding  catheters  within 
the  uterus  to  inject  fluid  into  the  fallopian  tubes  to  test  for  tubal  blockage. 
Ceramic  hollow  sphere  transducers  have  the  required  omnidirectionality  and 
high  sensitivity  and  can  be  electrically  matched  to  the  electronic  systems. 

It  is  our  hope  that  the  miniature  hollow  sphere  transducers  described  in 
this  paper  will  help  satisfy  the  size,  directionality  and  acoustic  impedance 
matching  requirements  for  biomedical  and  underwater  applications.  Later  we 
hope  to  manufacture  the  transducers  in  large  numbers  by  a  simple  inexpensive 
process,  which  makes  it  possible  to  mass-produce  throw-away  transducers. 
The  basic  understanding  of  vibration  modes  in  small  spherical  transducers 
provided  by  this  symmetry-guided  approach  will  help  immensely. 
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The  S-walls  (also  called  W'-walls)  refer  to  those  domain  walls  whose  orientations  do  not 
coincide  with  crystallographic  planes  but  rather  change  with  temperature.  A  simple  contin¬ 
uum  model  has  been  constructed  to  describe  the  S-wall  structures  that  occur  at  a 
m3m  mmm  ferroelastic  phase  transition.  The  model  describes  local  values  of  the  elas¬ 
tic  strain  components  and  the  temperature  dependence  of  the  S-wall  orientation.  Distortions 
of  unit  cells  across  an  S-wall  was  obtained  numerically. 

Keywords:  S-walls;  W-walls;  Landau-Ginzburg  model;  crystal  symmetry;  ferroelastics; 
domain  walls 

INTRODUCTION 

The  S-walls  (also  called  W-walls)  occur  in  certain  ferroic  species 
and  generally  do  not  coincide  with  crystallographic  planes.  The  orien¬ 
tation  of  S-walls  (SWs)  in  ferroelastics  depends  on  the  values  of 
spontaneous  strain,  hence,  will  change  with  temperature. 

Because  the  "strange"  physical  nature  and  the  high  application  po¬ 
tential  for  these  SW  structures,  there  have  been  intense  experimental 
studies  on  them  and  many  interesting  results  have  been  obtained.  For 
example,  intersections  between  S-walls  and  W-walls  in  ferroelectric 
KNb03  w,  temperature  induced  SW  switching  in  ferroelastic  phases  of 


[349J/137 


138/P50] 


JIRlERHART  etal. 


NaNb03 151  and  AgNb03  {6J.  S-walls  have  also  been  observed  in  many 
other  orthorhombic  materials,  such  as  PbZr03  [7],  PbHf03 
PKZri.xSnJOj 1,1  and  Pb(Ybo.sNbo.s)03  "01. 


FIGURE  1.  Distortion  near  the  S-wa!l  for  4/m->2/m  phase  transition1"1. 

One  earlier  work  on  the  SW  structure  was  done  for  a  4/m->2/m  fer- 
roelastic  phase  transition  by  using  a  two-dimensional  continuum 
model  lul.  More  recently,  high  resolution  transmission  electron  mi¬ 
croscopy  has  been  used  to  characterize  the  SWs  in  monoclinic  LaNb04 
(4/m—>2/m  ferroelastics)  ^12*. 

In  this  paper,  we  present  a  simple  continuum  model  for  the  S-walls 
resulting  from  an  /w3/w— >  nttrim  phase  transition,  and  also  provide 
numerical  solutions  for  the  S-wall  structure  and  strain  variations  near 
the  SWs. 


THEORETICAL  MODEL 

The  m3ro-»  mmm  ferroelastic  phase  transition  has  a  shear  strain  com¬ 
ponent  as  the  order  parameter;  the  other  components  of  the  elastic 
strain  tensor  serve  as  secondary  order  parameters. 
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We  consider  twin  structure  characterized  by  the  following  two 
spontaneous  strains P1, 


fs, 

^6 

(P 

fs3 

0 

0  N 

sll)  = 

0 

,  s(2!  = 

0 

V 

-S' 

1° 

0 

u 

-S< 

St  J 

In  order  to  match  the  elastic  strain  at  the  interface,  the  orientation  of 
the  domain  wall  plane  will  depend  on  the  values  of  spontaneous  strain 
according  to  the  following  relation 

(S,  - S3)(x,  -x3)  +  2 S6x2  =  0  (2) 

All  strain  components  used  here  are  expressed  in  the  same  cubic  coor¬ 
dinate  system  of  parent  paraelastic  phase. 

Simplified  free-energy  density  (invariant  in  mint)  can  be  written 
as 

F(nij,i?ij.k)  =  ccl(Til  +rj\  +ijI  +  tj2)  + 

<z7{tj}{ti25+7]1  +T)t  -2 -2 vl)}  (3) 

a30?4  +J?i  +Tll)  +  g{(c^-f  +&L)2  +($1)2}- 

CDCj 

The  linear  elastic  strains,  %,  are  defined  by  the  derivative  of  die  elas¬ 
tic  displacement  u 

1,  ■ 

ty®  j(uu+uj.i)  (4) 

and  expressed  using  the  Voigt  notation  tji=  r\u,  tj2=  t]22,  r}3=  t]33,  rj4= 
2tj 23 ,  tj5=  2rfi3,  rj6=  2r}t2.  All  the  expansion  coefficients  in  Equation  (3) 
are  assumed  to  be  temperature  independent  except  a;  =  cto  (T-To), 
with  clq>0.  In  order  to  reduce  the  number  of  independent  coefficients 
in  the  free  energy  expansion,  we  have  normalized  the  free  energy  den¬ 
sity  by  the  coefficient  of  the  second  term. 

SW  problem  is  translationally  invariant  along  the  wall  plane  and 
hence  all  quantities  characterizing  the  SW  should  depend  on  the  dis- 
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tance  from  SW  only.  We  choose  a  new  coordinate  system  with  one  of 
the  axes  (x'/)  perpendicular  to  the  SW,  and  the  remaining  two  axes  (x;, 
x ’3)  parallel  to  the  SW  plane  (Figure  2).  Due  to  this  choice  the  orienta¬ 
tion  of  the  coordinate  system  depends  on  the  values  of  spontaneous 
strains,  i.e.  on  temperature. 


S-wall 

FIGURE  2.  Relationship  between  the  new  and  old  coordinate  systems. 

The  corresponding  coordinate  transformation  can  be  chosen  in  a 
special  form  given  by 


/-«  \ 
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(6a,b) 


The  components  %  of  the  elastic  strain  in  the  free-energy  must  be 
transformed  to  the  new  coordinate  system  (denoted  by  a  prime, 

rf\  >7,2»J7,3*etc  )- 


S-WALL  SOLUTIONS 

The  continuum  SW  solution  can  be  derived  from  the  Euler-Lagrange 
equations  taking  into  account  proper  boundary  conditions 
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-|L  =  0,  (/,;  =  1,2,3)  (7) 

In  materials  without  defects  the  following  compatibility  relations  must 
also  be  satisfied 

M*=°.  (i,jXl.m,n  =  1,2,3)  (8) 

where  is  the  Levi-Civita’s  tensor  (permutation  symbol).  Equations 
(7)  and  (8)  are  given  in  invariant  tensor  form,  independent  of  the 
choice  of  the  coordinate  system.  In  the  rotated  coordinate  system  the 
compatibility  relations  are  reduced  to 

*7*33.11  =  0,  *7*23.1 1 =  ®  ^ 

because  all  strains  depend  on  x’  i  only  due  to  the  translational  symme¬ 
try. 

In  non-homogeneous  system,  five  components  of  the  spontaneous 
strain  in  Equation  (1)  (transformed  to  the  rotated  coordinates)  are 
equal  in  both  domains,  but  the  shear  strain  component  S’ 6  (order  pa¬ 
rameter)  has  opposite  signs  in  the  two  domains.  The  boundary  condi¬ 
tions  for  the  elastic  strains  should  match  spontaneous  strains  of  domain 
1  asx'r-4+eoand  domain  2  asx'i-+-co 


_8_  JF_ 
dxk  drjljk 


lim  *7 ’«  (**, )  =  Sf ,  lim  *7*6  (**, )  =  ±Sf  Act  =  1,2,. ..,5)  (10a,b) 


The  compatibility  relations  Equation  (9)  can  be  integrated  and  the 
strain  tensor  components  parallel  to  the  SW  (rf 2 ,if 3 ,tj\)  must  be 
constant  due  to  the  required  boundary  conditions  Equation  (10a).  For 
convenience,  we  introduce  the  dimensionless  strain  components 
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Dimensionless  Space  Variable  x"j 

FIGURE  3.  Normalized  strain  components  across  an  S-wall  at  a  specific  tem¬ 
perature  below  T0. 


Finally,  the  SW  problem  reduces  to  solving  the  following  three  coupled  dif¬ 
ferential  equations  (dimensionless  coordinate  x\  is  defined  by  x,"  =  ^j= ) 

yfg 

d2y 

=  Ra(y\>ys>y(,\  («  = 1,5,6)  (ii) 

The  functions  Rj,  R$,  Rg  are  third-order  polynomial  functions  in  terms 
of  the  relative  elastic  strains  y',,  y's,  and  y'g.  The  coefficient  of  each 
term  is  a  complicated  function  of  the  coefficients  cr/,  cr?,  and  as  and 
will  change  with  temperature.  Typical  solution  of  a  twin  with  an  SW  is 
given  in  Figure  3. 
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Dimensionless  Space  Variable  x,fi 


FIGURE  4.  Displacement  profile  across  an  S-wall  at  a  specific 
temperature  below  To. 

The  elastic  displacements  can  be  found  by  numerically  integrating 
the  corresponding  strain  components  in  rotated  coordinates,  i.e.  Equa¬ 
tion  (4).  Typical  solutions  of  the  elastic  displacements  are  shown  in 
Figure  4. 


SUMMARY 

A  simple  three-dimensional  continuum  model  has  been  developed  for 
the  ferroelastic  S-wall  structure  produced  by  a  m3m->  mmm  ferro- 
elastic  phase  transition.  The  elastic  strains  as  well  as  displacements 
near  an  S-wall  have  been  calculated  using  this  model.  In  general,  latti¬ 
ce  deformation  near  the  S-wall  not  only  contains  longitudinal  strain 
but  also  includes  nonzero  shear  strain  components.  Although  the  SW 
solution  is  only  numerical,  our  three-dimensional  model  provides  a 
unique  theoretical  tool  to  study  SW  structure  in  the  systems  produced 

through  a  m3m->  mmm  ferroelastic  phase  transition. 
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Abstract 

Current  piezoelectric  composite  materials  contain  two  or  more  phases  out  of  which  at  least  one  reveals  piezoelectric 
properties  in  itself.  We  show  that  this  is  in  fact  not  a  necessary  condition.  The  mechanism  of  the  linear  stress-polarization 
response  averaged  over  a  composite  sample  can  be  also  based  on  flexoelectric  properties  of  one  or  mom  constituents.  Proper 
shaping  of  the  composite  constituents  is  required,  such  that  the  system  as  a  whole  acquires  a  symmetry  allowing  for  nonzero 
piezoelectric  coefficients  even  if  none  of  the  components  is  piezoelectric.  Externally  applied  stress  is  transformed,  due  to 
proper  geometry  of  the  constituents  with  different  elastic  properties,  into  a  strongly  nonhomogeneous  distribution  of  induced 
strain.  Flexoelectric  properties  which  are,  by  symmetry,  allowed  in  all  materials,  transform  the  strain  gradient  into 
polarization.  The  proposed  piezoelectric  composite  falls  into  the  category  of  composites  with  product  properties  since  it 
involves  different  assets  of  the  phases  (elastic,  flexoelectric  and  dielectric)  and  the  interaction  between  the  phases, 
determining  the  inhomogeneous  distribution  of  stress,  is  essential  €  1999  Elsevier  Science  B.V.  All  rights  reserved. 


PACS:  77.84  -s:  77.90  -k 
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1.  Introduction 

Composites  are  multiple-phase  solids  which  com¬ 
bine  materials  of  different  chemical  composition  and 
macroscopic  properties  with  the  aim  to  produce  sam¬ 
ples  with  the  desired  average  response.  Figures  of 
merit  of  die  final  composite  can  be  tuned  by  choos¬ 
ing  component  phases  with  the  right  properties  and 
coupling  them  in  an  optimum  manner.  Newnham  et 
al.  [1]  offered  a  classification  of  composites  based  on 
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three  criteria.  The  most  important  aspects  are  the 
macroscopic  properties  of  the  constituents,  e.g,.  their 
response  to  electric,  magnetic  and  elastic  fields.  This 
determines  the  final  assets  of  the  composite.  The 
second,  connectivity,  indicates  the  way  in  which 
each  phase  connects  to  itself.  It  is  essential  for  the 
magnitude  and  symmetry  of  the  composite's  re¬ 
sponse.  The  third  is  scale,  which  determines  the 
response  of  the  composite  when  wavelengths  of 
propagating  waves  become  comparable  with  the 
characteristic  dimensions  of  any  of  the  constituents. 

Many  composites  have  been  considered  in  con¬ 
nection  with  their  piezoelectric  properties  [1*2).  To 
discuss  or  model  the  piezoelectric  response  of  a 
composite,  it  was  generally  assumed  that  at  least  one 

All  rights  reserved. 
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of  the  components  was  piezoelectric.  In  this  paper 

we  reconsider  this  assumption. 

*• 

2.  Discussion 

For  a  system  to  be  piezoelectric,  it  has  to  fulfill 
certain  symmetry  criteria.  If  it  has  a  crystalline  struc¬ 
ture.  the  material  must,  by  symmetry,  belong  to  one 
of  20  crystal  classes.  The  remaining  12  classes  do 
not  show  piezoelectric  properties;  these  are  the  11 
centrosymmetric  classes  and  the  class  432  in  which 
piezoelectricity  is  forbidden  by  the  combined  sym¬ 
metry  elements. 

In  Nye‘s  widely  used  overview  of  the  equilibrium 
properties  of  the  32  crystal  classes  [3],  the  properties 
of  an  isotropic  medium  are  also  included.  The  sym¬ 
metry  of  an  isotropic  medium  is  primarily  character¬ 
ized  by  the  presence  of  arbitrarily  oriented  symmetry 
axes  of  infinite  order.  Depending  on  whether  its 
symmetry  elements  do  or  do  not  include  arbitrarily 
oriented  mirror  planes,  the  isotropic  medium  repre¬ 
sents  one  of  the  Curie  symmetry  groups  (limiting 
groups),  namely,  =c,  oc /m  or  respectively.  In 


Table  1 

Matrices  of  d,jk  in  Curie  groups  which  are  piezoelectric 

5C 

0  0  0  du  0 

0  0  0  d j*  0 

d\-  d't,  j  d^i  0  0  0 

0  0  0  rf,4  0  0 

0  0  0  0  -  c/,4  0 

0  0  0  0  0  0 

0  0  0  0  </„  0 

■0  0  0  dl5  0  0 

d,_  dyx  dy,  0  0  0 


In  all  groups  the  *  axis  is  taken  as  .v; :  the  axes  .t,,  .t2  are 
perpendicular  to  xy  and  to  each  other,  otherwise  their  orientation 
is  arbitrary. 


co  m 


Fig.  1.  Characteristic  forms  representing  symmetry  of  the  Curie 
groups  which  allow  for  piezoelectricity  [4]. 


both  these  groups  no  nonzero  piezoelectric  coeffi¬ 
cients  are  possible. 

In  addition  to  these  two  groups,  however,  we  may 
consider  systems  representing  the  remaining  Curie 
groups,  namely  »,  *>/m,  »2,  *m  and  */mm.  Out  of 
these,  sc /m  and  »/ mm  do  not  allow  for  the  exis¬ 
tence  of  nonzero  components  of  a  third-rank  polar 
tensor  dri  of  symmetry  V  [V:].  i.e..  of  the  piezo¬ 
electric  tensor.  In  the  remaining  groups  nonzero 
components  are  possible,  as  shown  in  Table  1.  It  is 
useful  to  illustrate  symmetry  properties  of  these  point 
groups  by  characteristic  forms  [4]:  these  are  shown 
in  Fig.  1.  We  realize  that  systems  revealing  the 
symmetries  *  or  *>2  can  exist  in  two  forms,  left-  and 
right-handed. 

Thus,  for  instance,  a  composite  with  connectivity 
0-3.  in  which  the  phase  *0’  is  represented  by  cone- 
shaped  particles  whose  =c-axes  are  parallel  to  each 
other  but  which  are  randomly  distributed  in  the 
phase  *3\  has  the  symmetry  *m.  Next  we  can 
imagine  that  the  cones  are  subject  to  helical  deforma¬ 
tions  so  that  spiral-shaped  panicles  result.  This  low¬ 
ers  the  symmetry  to  ».  The  system  can  exist  in  two 
forms,  right-  or  left-handed.  The  third  piezoelectric 
Curie  group  can  be  visualized  stoning  again  with  a 
composite  of  connectivity  0-3  in  which  the  phase 
‘O’  is  represented  by  cylinders:  its  symmetry  is 
nonpiezoelectric  x>/rrun.  If  now  all  cylinders  are 
subject  to  a  helical  deformation,  the  symmetry  is 
reduced  to  ^2,  which  again  can  exist  in  two  forms 
differing  in  handedness.  Fig.  2  shows  such  compos¬ 
ites  schematically.  These  particular  models  are  based 
on  the  0-3  connectivity  but  similar  considerations 
can  be  made  for  other  connectivities  as  well. 

It  thus  appears  easily  possible  to  manufacture 
composites  whose  symmetry  properties  allow  for  the 
existence  of  piezoelectric  tensor  although  they  con¬ 
sist  of  components  which  by  themselves  need  not  be 
made  of  piezoelectric  materials. 
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Fig.  2.  Simple  models  of  0-3  composites  allowing  for  piezoelec¬ 
tricity,  There  is  an  infinite  number  of  shapes  of  the  0-constituent 
that  could  be  tested  for  a  maximum  response  of  flexoelectric 
polarization. 


A  long  time  ago,  Shubnikov  et  al.  introduced  the 
concept  of  piezoelectric  textures ;  see  Ref.  [5]  and 
the  first  chapter  of  Ref.  [6].  This  notion  denotes 
systems  composed  of  crystallites  which  show  piezo¬ 
electric  properties,  some  crystal  axes  of  which  are 
chaotically  oriented  in  a  given  way,  leading  to  spe¬ 
cific  averaging  of  properties  characterizing  the 
piezoelectric  effect.  In  the  subsequent  parts  of  Ref. 
[6],  two  specific  kinds  of  materials  are  discussed  as 
examples.  Plate-like  samples  containing  crystallites 
of  Rochelle  salt  separated  by  amorphous  layers 
(chapter  2  by  Konstantinova  and  SiTvestrova)  could 
indeed  be  considered  a  0-3  composite.  Polarized 
ceramic  samples  of  barium  titanate  (chapter  3  by 
Zheludev),  on  the  other  hand,  whose  symmetry  is 
s=m,  represent  a  piezoelectric  texture  by  Shubnikov’ s 
definition  but  could  only  be  included  into  the  family 
of  composites  if  the  grain  boundaries  had  an  appre¬ 
ciable  volume. 

In  both  these  cases  basic  components  are  piezo¬ 
electric  by  themselves  (point  symmetries  2  for 
Rochelle  salt  and  4mm  for  BaTi03).  On  the  other 
hand,  in  the  general  symmetry  approach  this  is  not  a 
specific  requirement  As  shown  above,  the  possibil¬ 
ity  of  a  piezoelectric  response  in  a  0-3  composite  is 
assured,  from  the  point  of  view  of  symmetry,  already 
by  shaping  the  particles  of  zero  connectivity.  We 
thus  have  to  look  for  alternative  mechanisms  which 
would  lead  to  formation  of  an  average  polarization 
proportional  to  an  applied  stress  for  a  0-3  composite 
made  of  nonpiezoelectric  materials,  with  properly 
shaped  particles. 

For  several  decades,  the  effect  of  inducing  polar¬ 
ization  by  imposing  spatially  nonuniform  strain  was 
repeatedly  discussed  in  the  literature.  Originally  dis¬ 
covered  experimentally  in  centrosymmetric  (and 
therefore  nonpiezoelectric)  liquid  crystals,  it  was 


termed  the  flexoelectric  effect  and  described  by  the 
equation 

dSij 

p,~ftisk!irk  (1) 

from  which  it  is  obvious  that  the  tensor  p  has  the 
general  symmetry  [V2]  F2,  Le.,  pijkl  =  pjikl  is  the 
only  requirement  imposed  by  symmetry.  A  tensor  of 
this  symmetry  has  nonzero  components  in  all  crystal 
classes.  The  first  attempt  to  observe  the  flexoelectric 
effect  in  a  solid  crystal  of  point  symmetry  4/m , 
namely  CaW04  was  made  by  Zheludev  et  al.  [7].  As 
shown  by  Tagantsev  [8],  the  static  effect  includes  a 
bulk  tod  a  surface  contribution.  The  bulk  part  is  due 
to  the  .fact  that  the  crystal  lattice  which  has  been 
nonhomogeneously  deformed  in  accordance  with  the 
laws  of  the  theory  of  elasticity  is  not  in  equilibrium 
from  the  point  of  view  of  displacements  in  the  unit 
cell.  The  displacements  that  are  necessary  to  reach 
true  equilibrium  give  rise  to  a  dipole  moment  of  the 
cell,  Le.,  to  polarization.  In  addition,  the  deformation 
of  the  surface  of  a  finite  sample,  whose  electrical 
neutrality  in  the  original  state  was  achieved  by  com¬ 
pensating  free  charges,  leads  to  a  surface  contribu¬ 
tion  which  can  be  expected  to  be  of  the  same  order 
of  magnitude  as  the  bulk  part  of  the  effect  The 
simplest  estimates  [8]  for  a  common  insulator  indi¬ 
cate  that  both  contributions  to  the  value  of  p  axe  of 
the  order  of  the  ratio  of  the  electron  charge  to  the 
lattice  constant 

We  now  have  in  mind  a  0-3  composite  made  of 
nonpiezoelectric  constituents,  in  which  the  0-ele- 
ments  are  shaped  and  oriented  in  such  a  way  that  the 
overall  symmetry  is  one  of  the  Curie  groups 
and  so2.  As  an  example,  consider  a  plate-like  sample 
of  composite  of  symmetry  in  which  the  orienta¬ 
tion  of  the  0-constituents  is  such  that  the  oo  axis  is 
perpendicular  to  the  major  plane.  Since  the  tensor 
pijki  has  nonzero  components  even  for  continuous 
groups  and  soc© m  [4],  one  can  imagine  that  both 

the  0-component  and  3-components  are  made  of 
isotropic  materials.  Their  shaping  is  such  that  when  a 
load  cr33appl  is  applied,  the  spatial  distribution  of 
stress  will  be  nonfaomogeneous,  leading  to  gradients 
of  strain  in  both  constituents.  To  be  concrete,  we  can 
imagine  that  a  plate-like  sample  of  thickness  d  is 
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divided  into  a  regular  system  of  cubes  of  linear  size 
d0>  each  cube  containing  one  particle  of  the  compo¬ 
nent  with  zero  connectivity  and  only  this  component 
is  assumed  to  have  nonzero  flexoelectric  properties. 
The  dipole  moment  of  each  0-particle  within  one 
cube  will  be  given  by 

r  deij(r ) 

*--4>'‘“*-£rdr  <2> 


where  the  strain  gradient  will  be  determined  by  a 
factor  aiJk  which  reflects  the  shape  of  the  compo¬ 
nents  0  of  the  composite  and  depends  on  elastic 
compliances  of  both  constituents: 


8g.y(r)  ,  v 

‘~aijk\r)  ^33appl  - 


(3) 


The  induced  charge  density  on  the  electrode  of  the 
plate  will  be 


so  that  the  effective  polarization  of  the  sample  P3  — 
Q  will  be  given  by 

^3  “  ^33flex  ^33appl  *  (5) 

where 


dents  of  its  constituent  phases)  and  product  proper¬ 
ties.  In  the  latter  case  the  composite  property  coeffi- 
tient  involves  different  properties  of  the  constituent 
phases  with  interactions  between  them.  It  appears 
that  piezoelectric  composites  based  on  flexoelectric- 
ity  fall  into  this  last  category  and  the  effect  might  be 
referred  to  as  a  ‘shape-controlled  product  property’; 
indeed  the  combined  effect  involves  different  proper¬ 
ties  of  die  constituent  phases  (elastic,  flexoelectric 
and  dielectric)  and  the  interaction  between  the  phases 
is  essential;  here  it  is  die  nonhomogeneous  distribu¬ 
tion  of  stress  which  depends  primarily  on  the  shapes 
of  constituents  and  on  their  elastic  tensors.  The 
following  sequence  of  phenomena  describes  the 
combined  effect:  homogeneous  applied  stress 
inhomogeneous  stress  in  the  O-constituents  — ► 
polarization  in  the  O-constituents  due  to  flexoelectric 
effect  — >  nonhomogeneous  distribution  of  polariza¬ 
tion  in  the  sample  depending  also  on  spatial  distribu¬ 
tion  of-  pennittivity  — >  nonhomogeneous  surface 
bound  charge  — ►  averaged  surface  bound  charge  den¬ 
sity  defining  effective  polarization.  Considering  a 
stress  <rzz  applied  perpendicularly  to  a  plate-like 
composite  sample,  we  have  the  sequence 

^cc.appi  grud pje ( r)  — *  P( r)  — >  ^(  x,y)surf 


^33flex  d3JidiaiJk(r)^iJ3dr.  (6) 

A  piezoelectric  composite  based  on  flexoelectricity 
will  be  useful  if  a  reasonably  high  value  of  d33flcx 
could  be  reached,  e.g.,  100  pC/N.  It  follows  from 
the  preceding  formulae  that  the  latter  can  be  influ¬ 
enced  by  a  proper  tuning  of  several  independent 
factors:  selecting  materials  with  high  values  of  those 
components  of  which  are  involved  in  a  particu¬ 
lar  geometry  of  the  constituents,  choosing  a  high 
density  of  the  constituents  0,  but  also  by  achieving 
large  factors  aijk  which  depend  on  the  shape  of  the 
0-components  and  on  the  elastic  tensors  of  both  0 
and  3  constituents. 

Newnham  [2]  classified  properties  of  composite 
materials  into  three  groups:  sum  properties  (the  com¬ 
posite  property  coefficient  depends  on  the  corre¬ 
sponding  coefficients  of  its  constituent  phases),  com¬ 
bination  properties  (the  composite  property  coeffi¬ 
cient  depends  on  two  or  more  corresponding  coeffi- 


*  4?surf  ®  ^.surf  ^  ^cr.appl  * 

At  this  stage  very  few  data  on  the  tensor  fiiJkl  in 
solids  seem  to  be  available.  A  fairly  strong  flexoelec¬ 
tric  response  was  reported  for  crystals  of  Cd2W04 
[7].  Marvan  and  Havrdnek  [9]  studied  the  flexoelec¬ 
tric  effect  in  elastomers  of  isotropic  symmetry.  Sam¬ 
ples  in  the  form  of  truncated  pyramids  were  de¬ 
formed  by  axial  pressure  along  the  axis  3.  Then  for 
constant  volume  of  the  sample  the  only  active  coeffi¬ 
cient  is  ^3333  which  was  estimated  to  be  of  the  order 
10”n  to  10~10  C/m.  Experiments  with  0-3  com¬ 
posites  in  which  the  0-constituent  or  both  compo¬ 
nents  of  the  composite  would  be  a  polymer  might  be 
worthwhile. 

We  may  also  note  that  such  composite  samples 
might  be  interesting  to  investigate  in  which  one  of 
the  constituents  is  piezoelectric;  due  to  flexoelectric¬ 
ity,  its  induced  dipole  moment  could  be  considerably 
enhanced  by  proper  shaping  to  optimize  the  nonho¬ 
mogeneous  distribution  of  strain. 
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The  relationship  between  elastic  strain  gradient  and  electric  polarization  was  investigated  in  the 
relaxor  ferroelectric  lead  magnesium  niobate  ceramic.  Experimental  studies  indicated  that 
flexoelectric  polarization  is  linearly  proportional  to  the  applied  strain  gradient,  and  the  flexoelectric 
coefficient  fil2  is  calculated  to  be  4X1CT6  C/m,  which  is  much  higher  than  the  early 
phenomenological  estimations.  ©  2001  American  Institute  of  Physics.  [DOI:  10.1063/L1356444] 


The  piezoelectric  effect  is  the  linear  relationship  be¬ 
tween  electric  polarization  and  applied  elastic  stress  (direct 
effect)  or  between  elastic  strain  and  applied  electric  field 
(converse  effect).  Due  to  the  symmetry ‘requirements,  the 
piezoelectric  effect  may  only  exist  in  20  crystal  point  groups 
where  center  of  symmetry  is  absent  More  than  three  decades 
ago,  phenomenological  analysis1  stated  that  electric  polariza¬ 
tion  might  be  induced  by  inhomogeneous  strain  in  the  mate¬ 
rials.  Later  the  term  flexoelectric  was  coined  to  describe  such 
an  effect  in  analogy  to  an  effect  discovered  in  the  nonpiezo- 
electric  liquid  crystals,2  and  tensor  properly3  and  micro¬ 
scopic  theory4  were  investigated  theoretically  for  single  crys¬ 
tals.  The  flexoelectric  effect  is  defined  by  the  following 
relation: 

d€n 

PnHjtrT**  (i) 

dxk 

where  Pi  is  the  /th  component  of  flexoelectric  polarization, 
Pijki  is  the  flexoelectric  coefficient,  €fj  is  the  elastic  strain, 
and  xk  is  the  position  coordinate.  Flexoelectricity  is  con¬ 
trolled  by  a  fourth-rank  tensor,  and  thus  is  a  property  of  all 
insulating  solids,  but  so  far  no  flexoelectric  coefficients  have 
in  fact  been  measured.  In  the  late  1960s,  attempts  were  made 
to  investigate  the  electrical  polarization  in  nonpiezoelectric 
CaW04  crystals  during  torsional  deformation,  which  might 
be  closely  associated  with  the  flexoelectric  effect,5  but  the 
polarization  had  been  attributed  to  the  disappearance  of  the 
center  of  symmetry  after  torsion,  rather  than  the  effect  of  the 
strain  gradient.  In  the  late  1980s,  a  tiny  flexoelectric  effect 
was  observed  in  an  isotropic  elastomer  with  a  static  flexo¬ 
electric  coefficient  of  the  order  of  magnitude  of 
10“ : li- 10" 1 10  C/m.6 

Recently,  an  idea  was  developed  for  designing  0-3  pi¬ 
ezoelectric  composites  based  on  die  flexoelectric  effect 
where  none  of  the  components  are  piezoelectric.7  If  the 
flexoelectric  coefficients  are  large  enough  for  practical  appli¬ 
cations  a  range  of  properly  engineered  flexoelectric  compos¬ 
ite  structures  could  provide  completely  new  piezoelectric  ca¬ 
pability.  Recent  studies8,9  have  indicated  that  the  relaxor 
ferroelectric  perovskites  like  PbCMgj^Nb^Os  (PMN)  are 
exceedingly  sensitive  to  elastic  stress  in  the  region  of  tem¬ 
perature  between  the  dielectric  permittivity  maximum  Tm 
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and  the  freezing  temperature  Tf  from  the  Vogei-Fulcher  re¬ 
lationship.  In  this  letter,  we  report  the  investigation  of  flexo¬ 
electric  effect  in  relaxor  PMN  ceramics. 

The  flexoelectric  measurements  were  earned  out  at  room 
temperature  by  using  an  experimental  setup  shown  schemati¬ 
cally  in  Fig.  1.  The  PMN  ceramic  bars  (76.3  mm  length,  12.7 
mm  width,  and  2.5  ram  thickness)  were  provided  by  die  TR3 
company  of  State  College,  Pennsylvania.  A  very  thin  layer  of 
sputtered  gold  was  used  as  electrode,  the  bottom  surface  of 
the  sample  is  fully  covered  with  gold  while  on  the  top  a 
series  of  3  mm  diameter  electrodes  were  prepared.  Strain 
gradient  was  generated  in  the  samples  along  the  thickness 
direction  by  a  loudspeaker.  A  Micro  Strain  DVRT  transducer 
was  used  to  measure  the  displacement  at  several  positions 
along  the  sample  bar.  The  loudspeaker  was  driven  by  an 
HP467A  power  amplifier  with  a  1  Hz  sinusoidal  signal  from 
an  SR830  DSP  lock-in  amplifier.  The  generated  current  was 
monitored  by  die  lock-in  amplifier.  The  generated  polariza¬ 
tion  may  be  calculated  using  the  following  equation: 

i 


where  i  is  the  alternating  current,  A  is  the  electrode  area,  and 
/is  the  driving  frequency  of  the  loudspeaker. 

The  macroscopic  symmetry  of  PMN  between  Tm  and  Tf 
is  cubic  m3m,  while  die  unpoled  ceramics  are  regarded  as 
isotropic.  Thus,  the  flexoelectric  tensor  matrices  may  be 
written  as 


Loudspeaker 


FIG.  1.  Experimental  setup  for  the  measurement  of  die  flexoelectric  effect. 
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FIG.  2.  Schematic  diagram  of  a  cantilevered  beam  and  the  axial  arrange¬ 
ment. 
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For  the  present  situation  in  Fig.  1,  only  the  transverse 
effect  exists,  so  the  flexoelectric  equation  may  be  simplified 
as 


A*i2“  •  (3) 

(9*3 

The  strain  in  the  PMN  ceramic  bar  was  calculated  by 
assuming  the  natural  vibration  of  a  cantilevered  beam  shown 
in  Fig.  2.  The  mode  shape10  can  be  written  as 

■ Ar[(sin  /Srl-sinh  fir  I)(sin  fir  x{ 

—  sinh  fir  *1)  +  (cos  /3rZ  +  cosh  fir  L) 


X(C0S  fir  Xj  —  cosh  fir  *,)],  (4) 


where  Ar~Cx  /(sin  fir  L -sinh  firL),  r=  1,2, - 

Here  we  only  consider  die  fundamental  mode  (r=l), 
therefore  we  have  fix  L=  1.875.  Cx  can  be  determined  from 
the  boundary  condition,  i.e.,  the  measured  vertical  displace¬ 
ment  of  the  PMN  bar. 

The  gradient  of  transverse  strain  along  the  thickness  di¬ 
rection  can  be  written  as11 


deu  __  d2W(xx ) 
<9*3  dx\ 


-0.004  0.000  0.004  0.008  0.012  0.016  0.020 

Strain  gradient  (m'1) 


FIG.  3.  The  relationship  between  flexoelectric  polarization  and  strain  gradi¬ 
ent  in  the  relaxor  PbCMg^Nb^^  ceramics. 

to  1  MHz  have  been  investigated  and  the  dielectric  peak 
temperatures  are  well  below  room  temperature.  At  room 
temperature,  the  materials  should  be  cubic  and  the  unpoled 
ceramics  are  isotropic  and  should  not  exhibit  any  piezoelec¬ 
tric  effect  in  principle.  Furthermore,  for  the  current  measure¬ 
ment  configuration,  the  upper  and  lower  parts  of  the  sample 
endure  strain  of  opposite  sign.  So  even  if  there  is  any  rem¬ 
nant  piezoelectric  effect,  the  net  effect  in  the  whole  sample 
should  be  zero  because  the  piezoelectric  polarization  gener¬ 
ated  in  the  upper  and  lower  halves  would  cancel. 

Former  theoretical  analysis  estimated  the  flexoelectric 
coefficient  to  be  of  the  order  of  magnitude  of  e/a , 1,3,4  where 
e  is  the  absolute  value  of  electron  charge  and  a  is  the  lattice 
parameter.  Generally,  e/a  is  around  10” 10  C/m  for  almost  all 
materials,  which  is  very  small.  The  current  experimental  re¬ 
sults  on  the  relaxor  PMN  ceramics  are  much  higher  than  the 
phenomenological  estimations  and  are  very  encouraging. 
Such  high  flexoelectric  effect  observed  in  the  relaxor  PMN 
ceramics  may  be  associated  with  die  high  dielectric  constant 
(around  13  000  at  1  kHz  and  room  temperature),  so  the 
charge  separation  becomes  much  easier  than  those  in  the  low 
permittivity  materials. 

In  summary,  a  flexoelectric  effect  has  been  investigated 
in  the  relaxor  ferroelectric  PMN  ceramic,  the  measured 
flexoelectric  coefficient  is  about  four  orders  of  magnitude 
higher  than  the  earlier  theoretical  estimation. 


Figure  3  showed  the  experimental  results  of  the  strain 
gradient  induced  polarization  at  two  positions  of  the  PMN 
bar,  one  is  near  die  clamped  end  (x{/L  =  0.1 8)  where  the 
strain  gradient  is  supposed  to  be  higher,  while  the  other  is 
near  the  free  end  ( xx  /L= 0.67).  It  is  seen  that  the  flexoelec¬ 
tric  polarization  is  proportional  to  the  strain  gradient.  The 
slopes  of  the  polarization  versus  strain  gradient  curve  ob¬ 
tained  at  two  positions  are  very  close.  At  the  same  driving 
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voltage,  the  clamped  end  generated  higher  polarization  than 
the  free  end.  Based  on  the  experimental  results,  jjl12  was 
calculated  to  be  about  4X  10“6  C/m. 
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Observation  of  the  flexoelectric  effect  in  relaxor  Pb(Mg1/3Nb2,3)03  ceramics 
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The  relationship  between  elastic  strain  gradient  and  electric  polarization  was  investigated  in  the 
relaxor  ferroelectric  lead  magnesium  niobate  ceramic.  Experimental  studies  indicated  that 
flexoelectric  polarization  is  linearly  proportional  to  the  applied  strain  gradient,  and  the  flexoelectric 
coefficient  fil2  is  calculated  to  be  4X10"  6  C/m,  which  is  much  higher  than  the  early 
phenomenological  estimations.  ©  2001  American  Institute  of  Physics,  [DOI:  10,1063/1.1356444] 


The  piezoelectric  effect  is  the  linear  relationship  be¬ 
tween  electric  polarization  and  applied  elastic  stress  (direct 
effect)  or  between  elastic  strain  and  applied  electric  field 
(converse  effect).  Due  to  the  symmetry  requirements,  the 
piezoelectric  effect  may  only  exist  in  20  crystal  point  groups 
where  center  of  symmetry  is  absent  More  than  three  decades 
ago,  phenomenological  analysis1  stated  that  electric  polariza¬ 
tion  might  be  induced  by  inhomogeneous  strain  in  the  mate¬ 
rials.  Later  the  term  flexoelectric  was  coined  to  describe  such 
an  effect  in  analogy  to  an  effect  discovered  in  the  nonpiezo¬ 
electric  liquid  crystals,2  and  tensor  property3  and  micro¬ 
scopic  theory4  were  investigated  theoretically  for  single  crys¬ 
tals.  The  flexoelectric  effect  is  defined  by  fixe  following 
relation: 

d€t  / 

P'ijkl  ?  (1) 

9xk 

where  Pt  is  the  /th  component  of  flexoelectric  polarization, 
fiijki  *s  &e  flexoelectric  coefficient,  is  fixe  elastic  strain, 
and  xk  is  the  position  coordinate.  Flexoelectricity  is  con¬ 
trolled  by  a  fourth-rank  tensor,  and  thus  is  a  property  of  all 
insulating  solids,  but  so  far  no  flexoelectric  coefficients  have 
in  fact  been  measured.  In  the  late  1960s,  attempts  were  made 
to  investigate  the  electrical  polarization  in  nonpiezoelectric 
CaW04  crystals  during  torsional  deformation,  which  might 
be  closely  associated  with  the  flexoelectric  effect,5  but  fire 
polarization  had  been  attributed  to  the  disappearance  of  the 
center  of  symmetry  after  torsion,  rather  than  the  effect  of  the 
strain  gradient.  In  the  late  1980s,  a  tiny  flexoelectric  effect 
was  observed  in  an  isotropic  elastomer  with  a  static  flexo¬ 
electric  coefficient  of  the  order  of  magnitude  of 
10"  n- ID" 10  C/m.6 

Recently,  an  idea  was  developed  for  designing  0-3  pi¬ 
ezoelectric  composites  based  on  the  flexoelectric  effect 
where  none  of  the  components  are  piezoelectric.7  If  the 
flexoelectric  coefficients  are  large  enough  for  practical  appli¬ 
cations  a  range  of  properly  engineered  flexoelectric  compos¬ 
ite  structures  could  provide  completely  new  piezoelectric  ca¬ 
pability.  Recent  studies8*9  have  indicated  that  the  relaxor 
ferroelectric  perovskites  like  Pb(Mgl/3Nb2^)03  (PMN)  are 
exceedingly  sensitive  to  elastic  stress  in  the  region  of  tem¬ 
perature  between  the  dielectric  permittivity  maximum  Tm 
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and  the  freezing  temperature  Tf  from  the  Vogel-Fulcher  re¬ 
lationship.  In  this  letter,  we  report  the  investigation  of  flexo¬ 
electric  effect  in  relaxor  PMN  ceramics. 

The  flexoelectric  measurements  were  carried  out  at  room 
temperature  by  using  an  experimental  setup  shown  schemati¬ 
cally  in  Fig.  L  The  PMN  ceramic  bars  (76.3  mm  length,  12.7 
mm  width,  and  2.5  mm  thickness)  were  provided  by  die  TRS 
company  of  State  College,  Pennsylvania.  A  very  thin  layer  of 
sputtered  gold  was  used  as  electrode,  the  bottom  surface  of 
the  sample  is  fully  covered  with  gold  while  on  the  top  a 
series  of  3  mm  diameter  electrodes  were  prepared.  Strain 
gradient  was  generated  in  the  samples  along  the  thickness 
direction  by  a  loudspeaker,  A  MieroStrain  DVRT  transducer 
was  used  to  measure  the  displacement  at  several  positions 
along  die  sample  bar.  The  loudspeaker  was  driven  by  an 
HP467A  power  amplifier  with  a  1  Hz  sinusoidal  signal  from 
an  SR830  DSP  lock-in  amplifier.  The  generated  currant  was 
monitored  by  the  lock-in  amplifier.  The  generated  polariza¬ 
tion  may  be  calculated  using  die  following  equation: 

i 


where  i  is  the  alternating  current,  A  is  die  electrode  area,  and 
/is  the  driving  frequency  of  the  loudspeaker. 

The  macroscopic  symmetry  of  PMN  between  Tm  and  Tf 
is  cubic  m3m,  while  the  unpoled  ceramics  are  regarded  as 
isotropic.  Ihus,  the  flexoelectric  tensor  matrices  may  be 
written  as 


Loudspeaker 


FIG,  I.  Experimental  setup  for  the  measurement  of  the  flexoelectric  effect. 
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FIG.  2.  Schematic  diagram  of  a  cantilevered  beam  and  the  axial  arrange- 
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For  the  present  situation  in  Fig.  1,  only  the  transverse 
effect  exists,  so  the  flexoelectric  equation  may  be  simplified 
as 

de 

^3  =  ^12^-  0) 

dx3 

The  strain  in  the  PMN  ceramic  bar  was  calculated  by 
assuming  the  natural  vibration  of  a  cantilevered  beam  shown 
in  Fig.  2.  The  mode  shape10  can  be  written  as 

JFr(xi)=j4r[(sin  ^X-sinh  I)(sin  prx{ 

—  sinh  f}r  +  (cos  f3r  L  +  cosh  pr  L) 

X(C0S  Xi” COSh  Xj)],  (4) 

where  Ar-C{  /(sin  /?rZ-sinh  L),  r=  1,2, _ 

Here  we  only  consider  the  fundamental  mode  (r=l), 
therefore  we  have  fii  L-  1.875.  Cx  can  be  determined  from 
the  boundaty  condition,  i.e.,  the  measured  vertical  displace¬ 
ment  of  the  PMN  bar. 

The  gradient  of  transverse  strain  along  the  thickness  di¬ 
rection  can  be  written  as11 


d€\i  d2W{x  j) 
dx3  dx\ 


(5) 


FIG.  3.  The  relationship  between  flexoelectric  polarization  and  strain  gradi¬ 
ent  in  the  relaxor  PMMg^Nb^)^  ceramics. 


to  1  MHz  have  been  investigated  and  the  dielectric  peak 
temperatures  are  well  below  room  temperature.  At  room 
temperature,  the  materials  should  be  cubic  and  the  unpoled 
ceramics  are  isotropic  and  should  not  exhibit  any  piezoelec¬ 
tric  effect  in  principle.  Furthermore,  for  die  current  measure¬ 
ment  configuration,  the  upper  and  lower  parts  of  the  sample 
endure  strain  of  opposite  sign.  So  even  if  there  is  any  rem¬ 
nant  piezoelectric  effect,  the  net  effect  in  the  whole  sample 
should  be  zero  because  the  piezoelectric  polarization  gener¬ 
ated  in  the  upper  and  lower  halves  would  cancel. 

Former  theoretical  analysis  estimated  the  flexoelectric 
coefficient  to  be  of  the  order  of  magnitude  of  el  a, 1,3,4  where 
e  is  the  absolute  value  of  electron  charge  and  a  is  the  lattice 
parameter.  Generally,  el  a  is  around  10“ 10  C/m  for  almost  all 
materials,  which  is  very  small.  The  current  experimental  re¬ 
sults  on  the  relaxor  PMN  ceramics  are  much  higher  than  the 
phenomenological  estimations  and  are  very  encouraging. 
Such  high  flexoelectric  effect  observed  in  the  relaxor  PMN 
ceramics  may  be  associated  with  the  high  dielectric  constant 
(around  13  000  at  1  kHz  and  room  temperature),  so  the 
charge  separation  becomes  much  easier  than  those  in  the  low 
permittivity  materials. 

In  summary,  a  flexoelectric  effect  has  been  investigated 
in  the  relaxor  ferroelectric  PMN  ceramic,  the  measured 
flexoelectric  coefficient  is  about  four  orders  of  magnitude 
higher  than  the  earlier  theoretical  estimation. 


Figure  3  showed  the  experimental  results  of  the  strain 
gradient  induced  polarization  at  two  positions  of  the  PMN 
bar,  one  is  near  the  clamped  end  (xj /L  =  0.18)  where  the 
strain  gradient  is  supposed  to  be  higher,  while  the  other  is 
near  the  free  end  (xj  IL  =  0.67).  It  is  seen  that  the  flexoelec¬ 
tric  polarization  is  proportional  to  the  strain  gradient.  The 
slopes  of  the  polarization  versus  strain  gradient  curve  ob¬ 
tained  at  two  positions  are  very  close.  At  the  same  driving 
voltage,  the  clamped  end  generated  higher  polarization  than 
the  free  end.  Based  on  the  experimental  results,  yu12  was 
calculated  to  be  about  4X 10" 6  C/m. 

The  temperature  dependence  of  dielectric  permittivity 
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Large  flexoelectric  polarization  in  ceramic  lead  magnesium  niobate 

Wenhui  Ma  and  L.  Eric  Cross3^ 

-  Materials  Research  Laboratory,  Pennsylvania  State  University,  University  Park,  Pennsylvania  16802 
(Received  6  August  2001;  accepted  for  publication  4  October  2001) 

Flexoelectric  coefficient  fin  is  greatly  enhanced  in  the  relaxor  ferroelectric  lead  magnesium  niobate 
ceramic.  Phenomenological  analysis  suggests  the  high  dielectric  permittivity  is  inadequate  to 
explain  the  great  enhancement.  Temperature  dependent  measurement  reveals  a  close  relation 
between  the  flexoelectric  polarization  and  the  preexisting  polar  microregions  in  this  relaxor 
ferroelectries.  It  is  proposed  that  the  strain  gradient  might  change  the  Gibbs  free  energy  of  the 
relaxor  system  and  easily  reorient  the  already  existing  polar  microregions,  leading  to  the  greatly 
enhanced  flexoelectric  effect,  ©  2001  American  Institute  of  Physics.  [DOI:  10. 1063/1.1 426690] 


Flexoelectricity  describes  the  generation  of  electric  po¬ 
larization  In  an  insulating  solid  by  an  elastic  strain  gradient. 
Phenomenological  treatment  of  the  flexoelectric  effect  was 
first  made  by  Kogan,1  who  proposed  that  piezoelectric  polar¬ 
ization  in  an  inhomogeneously  deformed  medium  should  be 
a  function  not  only  of  deformation  but  also  of  the  spatial 
derivatives  of  the  deformation  at  the  same  point.  For  a  cen- 
trosymmetrie  crystal  subjected  to  inhomogeneous  deforma¬ 
tion,  only  the  strain  gradient  contributes  to  die  polarization, 
so  the  flexoelectric  effect  can  be  investigated  independently. 
The  flexoelectric  coefficient  is  a  fourth-rank  polar  ten¬ 
sor,  so  that  for  a  cubic  crystal  the  nonzero  components  are 
Pim*  ^U22>  Pun,  or  in  matrix  notation  fiUj  fin, 
and  fiM .  For  homogeneous  solids,  all  earlier  investigators 
agree  that  die  effect  is  very  small  in  simple  low  permittivity 
dielectrics,  and  the  flexoelectric  coefficients  are  of  the 
order  of  10“ 10- 10“ 11  C/m.2“4  In  a  previous  article,  we  an¬ 
nounced  our  measurement  of  the  flexoelectric  effect  in  lead 
magnesium  niobate  (PMN)  ceramic,  where  the  fil2  coeffi¬ 
cient  is  greatly  enhanced,  about  4X10~6C/m  at  room 
temperature,5  In  this  letter,  we  report  the  temperature  depen¬ 
dent  measurement  of  the  flexoelectric  polarization  in  PMN, 
Temperature  dependent  flexoelectric  measurement  was 
performed  by  using  an  experimental  setup  schematically 
shown  In  Fig.  1.  The  measurement  configuration,  sample  di¬ 
mension,  and  electrode  geometry  are  the  same  as  those  used 
previously  5  As  shown  In  Fig.  1,  one  end  of  the  sample  was 
fixed  to  form  a  cantilevered  beam.  Temperature  is  uniform 
all  over  the  sample  bar.  A  loudspeaker  was  used  to  drive  the 
beam  at  a  frequency  of  1  Hz,  so  a  natural  vibration  was 
generated  at  the  fundamental  mode.  By  measuring  the  mode 
shape,  we  can  calculate  the  strain  all  over  the  sample  bar. 
Because  the  top  and  bottom  surfaces  are  subjected  to  strain 
of  opposite  sign,  there  is  a  strain  gradient  along  the  thickness 
direction.  Our  former  measurements5  show  that  the  flexo¬ 
electric  polarization  is  proportional  to  the  applied  strain  gra¬ 
dient  and  the  fxn  coefficients  obtained  from  measurements  at 
different  positions  along  the  beam  are  close. 

Figure  2  shows  the  experimental  result  of  flexoelectric 
polarization  as  a  function  of  temperature  obtained  in  the 
PMN  ceramic.  There  is  a  maximum  around  the  dielectric 
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maximum  temperature,  Tm .  Above  Tm ,  the  flexoelectric  po- 
larization  decreases  rapidly  with  increasing  temperature.  As 
shown  in  Fig.  2,  an  exponential  decay  function  fits  the  ex¬ 
perimental  data  well.  By  extrapolation,  it  is  found  that  the 
flexoelectric  polarization  becomes  close  to  zero  around  Td , 
the  so-called  Bums  temperature.6  At  this  temperature,  the 
polar  microregions  disappear  and  the  relaxor  reverts  to  clas¬ 
sical  Curie-Weiss  behavior.  Our  result  indicates  that  there  is 
a  close  relation  between  die  flexoelectric  polarization  and  the 
already  existing  polar  microregions  in  PMN,  suggesting  that 
the  preexisting  polar  microregions  might  play  a  dominant 
role  in  the  polarization  mechanism. 

hi  order  to  understand  the  origin  of  the  greatly  enhanced 
Pn  coefficients  in  relaxor  PMN,  let  us  make  a  simple  analy¬ 
sis  of  the  contributors  to  the  flexoelectric  effect.  In  the  mid 
1980!s,  Tagantsev  developed  the  theory  for  flexoelectric  ef¬ 
fect  in  single  crystals.23  In  a  crystal  whose  bulk  symmetry 
does  not  allow  the  piezoelectric  interaction  the  appearance  of 
a  polarization  proportional  to  a  strain  gradient  can  be  com¬ 
prised  of  four  components:  a  bulk  static  flexoelectric  effect;  a 
bulk  dynamic  flexoelectric  effect;  a  surface  flexoelectric  ef¬ 
fect,  and  a  surface  piezoelectric  effect.  As  shown  in  the  dis¬ 
cussion  by  Tagantsev,  for  a  simple  homogeneous  ionic  solid 
of  cubic  centric  symmetry  based  on  the  rigid  ion  model,  all 
four  components  of  the  flexoelectric  response  are  of  the  or- 


HG.  1.  Experimental  setup  for  the  temperature  dependent  flexoelectric  mea¬ 
surement:  (1)  PMN  bar,  (2)  loudspeaker,  (3)  driving  arm,  (4)  thermocouple, 
(5)  microstrain  transducer  core,  (6)  high  purity  nitrogen,  (7)  gas  flow  meter, 
(8)  copper  coils  immersed  in  liquid  nitrogen,  (9)  heating  elements. 
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FIG.  2.  Flexoelectric  polarization  as  a  function  of  temperature  in  the  PMN 
ceramic. 


der  of  10" 10  C/m,  although  the  magnitude  of  surface  piezo¬ 
electric  effect  depends  markedly  on  the  nature  of  the  surface. 
For  a  soft-mode  ferroelectric  dielectric  above  Tc  in  the  cubic 
phase  the  flexoelectric  behavior  is  much  simplified,  both 
static  and  dynamic  bulk  coefficients  become  larger  in  propor¬ 
tion  to  the  enhanced  dielectric  susceptibility,  and  the  surface 
flexoelectric  coefficient  becomes  unimportant  as  it  does  not 
scale  with  the  susceptibility.  The  surface  piezoelectric  effect 
might  also  be  enlarged  in  proportion  to  the  dielectric  suscep¬ 
tibility,  but  again  depends  on  details  of  the  surface  character¬ 
istics.  For  the  current  measurement  configuration,  the  piezo¬ 
electric  polarization  from  the  upper  and  lower  halves 
neutralize,  so  there  should  be  no  net  effect  even  in  the  low 
temperature  ferroelectric  phase.  We  have  also  performed  the 
flexoelectric  measurements  with  the  sample  bar  turned  over, 
i.e.,  with  the  upper  and  lower  surfaces  interchanged,  the  re¬ 
sults  obtained  for  the  two  opposite  orientations  are  found  to 
be  in-phase,  evidencing  the  absence  of  a  piezoelectric  effect 
due  to  unbalanced  inhomogeneous  surface  layers. 

Figure  3  displays  the  dielectric  response  of  the  PMN 
ceramic  sample.  As  shown  in  the  dielectric  spectrum,  the 
dielectric  maximum  temperature  increases  with  frequency, 
which  is  typical  for  a  relaxor  ferroelectric  dielectric.  Permit¬ 
tivity  of  PMN  at  room  temperature  is  of  the  order  of  15  000, 
while  the  dielectric  permittivity  for  a  simple  ionic  solid  on 
the  point  charge  model  is  about  5,  so  an  enhancement  of 
3  X  103  can  be  expected  from  the  high  dielectric  permittivity. 
However,  the  flexoelectric  coefficient  jjin  of  PMN  is  of  the 
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FIG.  3.  Dielectric  permittivity  and  loss  tangent  vs  temperature  for  the  PMN 
ceramic. 
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FIG.  4.  Flexoelectric  coefficient  fil2  vs  dielectric  susceptibility  \  in  the 
PMN  ceramic,  where  the  solid  circles  are  the  experimental  data,  the  dashed 
curve  is  the  polynomial  fit,  and  the  solid  line  is  the  estimated  fin  values  by 
phenomenological  analysis. 

order  of  5  X 10"  6  C/m,  some  5  X 104  times  larger  than  that  of 
the  simple  ionic  solid.  It  is  clear  that  the  high  dielectric  per¬ 
mittivity  is  inadequate  to  explain  the  observed  larger  /x12 
value.  Figure  4  shows  the  curve  of  measured  fin  coefficient 
versus  dielectric  susceptibility  and  the  solid  line  is  from 
the  phenomenological  analysis  — ^(e/zz) ,  where  e  is 

the  elementary  charge  and  a  is  the  dimension  of  lattice  cell]. 
It  is  found  that  there  is  discrepancy  between  the  experimen¬ 
tally  observed  fxn  values  and  the  linear  estimations,  which 
increases  with  the  increase  of  dielectric  susceptibility  or  de¬ 
crease  of  temperature  below  Td .  It  may  be  noted  that  similar 
deviations  from  linear  extrapolation  of  high  temperature 


(c)  IP-0  (d)  I  p  *0 

FIG.  5.  Schematic  illustrations  of  (a)  eight  equivalent  <111)  polarization 
orientations  in  a  cubic  lattice  cell,  (b)  Gibbs  free  energy  A G  vs  polarization 
for  four  of  the  eight  equivalent  <11 1}  orientations  in  a  rhombohedral  ferro¬ 
electric  phase,  and  cross  sections  of  the  PMN  bar  showing  the  arrangement 
of  orientation  of  the  nanopolar  regions  in  a  [001]-oriented  grain  (defined  by 
the  dotted  ellipse)  (c)  before  and  (d)  after  being  subject  to  an  elastic  strain 
gradient. 
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properties  have  been  found  during  the  temperature- 
dependent  measurements  of  optic  index  of  refraction6  and 
thermal  expansion7  in  PMN,  which  were  unanimously  attrib¬ 
uted  to  the  appearance  of  localized  polarization. 

Relaxor  ferroelectrics  are  inhomogeneous  media  con¬ 
taining  mieropolar  regions  even  at  temperatures  well  above 
Tm ,  and  the  existence  of  nanoscale  B-site  order  has  been 
directly  confirmed  by  the  TEM  investigations.8  The  dielectric 
behavior  in  relaxor  ferroelectrics  can  be  understood  by  the 
combination  of  die  compositional  heterogeneity  model9  and 
superparaelectrie  model.7  Above  Td,  PMN  is  cubic  both  glo¬ 
bally  and  locally  without  any  polar  regions.  Figure  5(a) 
shows  a  lattice  cell  of  cubic  paraelectrie  phase  of  perovskite 
structure  displaying  eight  equivalent  (111)  orientations.  With 
temperature  going  down  below  Tdi  localized  polarization 
appears.  According  to  the  superparaelectrie  model,  it  is  ex¬ 
pected  that  the  nanopolar  regions  are  dynamically  disordered 
in  the  thermal  field  and  the  chance  of  appearance  of  a  polar 
vector  along  any  of  the  eight  (111)  directions  is  equivalent. 
Figure  5(b)  shows  the  Gibbs  free  energy  for  four  of  the  eight 
(111)  polarizations  for  a  low  temperature  ferroelectric  phase, 
it  is  clear  that  all  energy  minima  are  of  the  same  depth  in  the 
infinite  perfect  crystal.  However,  in  the  relaxor  ferroelectrics, 
due  to  composition  fluctuation  it  is  expected  that  the  local 
symmetry  will  be  lower  than  the  global  symmetry  so  that  the 
eight  (111)  polarization  orientations  are  now  inequivalent  in 
the  local  polar  microregions.  Under  zero  strain  condition,  a 
possible  arrangement  for  the  orientations  of  the  nanosized 
domains  in  a  [001]-oriented  grain  in  the  PMN  bar  is  sche¬ 
matically  shown  in  Fig.  5(c),  the  sum  of  all  those  polar  vec¬ 
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tors  is  zero.  In  Fig.  5(d),  the  sample  bar  is  subject  to  an 
elastic  strain  gradient  whose  direction  is  defined  as  pointing 
from  contraction  to  expansion  and  is  parallel  to  the  normal 
direction  of  sample  surface.  When  the  PMN  bar  is  bent, 
some  polarization  orientations  become  free  energy  unfavor¬ 
able  and  it  is  suggested  that  the  inhomogeneously  deformed 
nanosized  polar  regions  reorient  themselves  towards  the  di¬ 
rection  of  strain  gradient  by  polarization  rotation  between 
adjoining  (111)  directions,  e.g.,  from  [111]  to  [ill],  so  as 
to  reach  the  Gibbs  free  energy  minima  and  there  will  be  net 
polarization  in  the  strain  gradient  direction. 

In  summary,  flexoelectric  coefficient  fil2  has  been  mea¬ 
sured  as  a  function  of  temperature  and  the  results  suggest 
that  the  observed  large  flexoelectric  polarization  might  be 
caused  by  reorientation  of  the  preexisting  polar  clusters  in 
the  relaxor  ferroelectrics  by  the  elastic  strain  gradient 
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search,  Contract  No.  N00014-99-1-1011.  We  thank  M.  Mar- 
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the  PMN  ceramics. 
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Flexoelectric  polarization  of  barium  strontium  titanate 
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The  strain-gradient-induced  polarization  (flexoelectric  effect)  was  investigated  in  Bao  67Sr0  33Ti03 
(BST)  ceramic  at  temperatures  above  the  21  °C  Curie  point.  At  23  °C  the  flexoelectric  coefficient 
fin  was  more  than  one  order  of  magnitude  greater  than  the  highest  value  measured  in  lead 
magnesium  niobate  ceramic.  Over  the  temperature  range  of  linear  Curie- Weiss  behavior,  the 
coefficient  fiu  was  roughly  proportional  to  the  dielectric  permittivity;  however,  the  constant  of 
proportionality  was  higher  than  predicted  for  simple  ionic  solids.  The  unexpected  behavior  in  the 
BST  ceramic  suggests  the  need  for  a  broader  database  of  flexoelectric  coefficients.  ©  2002 
American  Institute  of  Physics.  [DOI:  10.1063/1.1518559] 


Recently  we  investigated  the  flexoelectric  polarization  in 
lead  magnesium  niobate  (PMN)  ceramic,1,2  a  material  which 
is  often  regarded  as  the  prototype  for  a  broad  family  of  re- 
laxor  ferroelectrics.3  The  high  values  measured,  and  particu¬ 
larly  the  observed  temperature  dependence  of  the  flexoelec¬ 
tric  coefficient  fin  suggested  that  the  polarization  was 
associated  with  a  reorientation  of  micropolar  regions  which 
are  known  to  occur  in  the  temperature  region  above  the  weak 
field  dielectric  maximum,  but  below  a  so-called  Bums  tem¬ 
perature  Td  which  is  just  the  temperature  range  where  the 
high  n  12  values  are  observed.  Alternatively  the  high  response 
may  be  a  consequence  of  the  very  high  dielectric  permittivity 
in  the  relaxor  where  the  relative  permittivity  sr  reaches  up  to 
13  000  (measured  at  1  kHz)  at  room  temperature.  Tagantsev4 
suggested  that  the  fitj  may  in  fact  be  proportional  to  the 
dielectric  susceptibility  Xij  following  a  relation  of  the  form 

torrxtq*  (l) 

where  y  is  a  constant  of  value  close  to  unity,  e  is  the  electron 
charge,  and  a  the  lattice  parameters. 

To  test  the  possible  role  of  micropolar  regions  it  was 
desirable  to  measure  the  response  of  a  “normal”  paraelectric 
cubic  perovskite  with  comparable  or  higher  dielectric  permit¬ 
tivity  but  no  lower  frequency  dielectric  dispersion  in  the  tem¬ 
perature  region  above  Tc.  Barium  strontium  titanate  solid 
solutions  have  perovskite  structure,  very  high  non-dispersive 
dielectric  permittivity,  and  Curie  temperature  Tc  which  can 
be  positioned  just  below  room  temperature  by  appropriate 
choice  of  composition.  It  must  be  cautioned  however  that 
this  is  an  exceedingly  difficult  system  to  synthesize  as  a  ho¬ 
mogeneous  solid  solution.  We  were  fortunate  to  be  able  to 
use  a  composition  fabricated  by  Texas  Instruments  (IT),  us¬ 
ing  a  complex  wet  chemical  route  to  the  very  demanding 
specification  for  their  long  wavelength  infrared  dielectric  bo- 
lometric  imager.  In  that  system  local  homogeneity  is  essen¬ 
tial  to  maintain  the  weak  first  order  transition  at  Tc  which 
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keeps  er  at  high  value  under  dc  bias  fields,  and  global  uni¬ 
formity  is  necessary  to  avoid  fixed  pattern  noise  in  the  de¬ 
tector  pixel  array. 

The  BST  ceramic  sample  used  in  this  study  was  cut  from 
a  boule  supplied  by  H  and  was  in  the  form  of  a  bar  68  mm 
long,  12.7  mm  wide,  and  1.8  mm  thick.  As  for  the  earlier 
measurements  on  PMN  ceramics2  3  mm  diameter  thin  sput¬ 
tered  gold  electrodes  were  spaced  along  the  bar.  Dielectric 
measurements  were  carried  out  in  a  delta  design  temperature 
chamber  to  check  local  and  global  uniformity.  Typical  data 
for  permittivity  and  loss  as  a  function  of  temperature  and 
frequency  are  given  in  Fig.  1  and  the  corresponding  recipro¬ 
cal  permittivity  (dielectric  stiffness)  in  Fig.  2.  The  high  qual¬ 
ity  of  the  sample  which  has  a  Curie  point  at  21  °C  is  attested 
by  the  very  high  sharp  permittivity  peak  at  Tc ,  the  low  di¬ 
electric  loss  and  the  complete  absence  of  dielectric  disper¬ 
sion  over  the  frequency  range  from  100  Hz  to  100  kHz  in  the 
paraelectric  range.  Polarization  measurements  confirm  that 
the  Curie  transition  is  weakly  first  order.  It  is  important  to 
note  the  small  upturn  in  dielectric  loss  below  35  °C  (Fig.  1) 
which  correlates  closely  with  the  departure  from  Curie- 
Weiss  behavior.  We  suspect  that  even  in  these  TI  prepared 


FIG.  1.  Dielectric  permittivity  and  loss  in  Bao,67Sr033H03  ceramic  as  a 
function  of  temperature,  measured  at  frequencies  of  100  Hz,  1  kHz,  10  kHz, 
and  100  kHz  at  a  field  level  of  5  V/cm. 
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FIG.  2.  Curie-Weiss  plot  of  reciprocal  dielectric  permittivity  (stiffness)  as  a 
function  of  temperature  in  Bao.^Sro^TiC^ . 


samples  there  is  a  weak  persistance  of  macroferroelectric  re¬ 
gions  above  Tc  reducing  permittivity  from  the  expected 
Curie-Weiss  behavior  and  enhancing  dielectric  loss. 

The  system  used  for  flexoelectric  measurements  is  iden¬ 
tical  to  that  used  in  earlier  studies  of  PMN  and  is  reproduced 
here  in  schematic  form  in  Fig.  3.  The  ceramic  bar  sample  is 
rigidly  clamped  at  one  end  and  driven  into  transverse,  vibra¬ 
tion  at  the  “free”  end  by  a  small  moving  coil  loudspeaker. 
Current  generated  from  the  small  gold  electrode  is  measured 
by  a  lock-in  amplifier,  phase  locked  to  the  loudspeaker  driv¬ 
ing  frequency  of  1  Hz.  A  Microstrain  DVRT  (differential 
variable  reluctance  transducer)  is  used  to  measure  the  ac 
strain  as  a  function  of  position  along  the  bar.  For  the  outer 
1/3  of  the  bar  length  the  measured  vertical  displacement 
curve  correlates  closely  with  the  curve  calculated  for  the  first 
fundamental  mode  of  a  freely  vibrating  cantilevered  bar.1 
Close  to  the  clamped  end,  the  measured  amplitude  is  slightly 
larger,  probably  indicating  that  our  best  clamping  is  not  suf¬ 
ficiently  rigid  for  the  stiff  BST  ceramic.  The  free  bar  model, 
with  the  measured  amplitudes  was  used  to  calculate  the 
strain  gradient  at  the  positions  of  the  electrodes. 


FIG.  3.  Experimental  setup  for  the  temperature  dependent  flexoelectric  mea¬ 
surement:  (1)  BST  bar,  (2)  loudspeaker,  (3)  driving  arm,  (4)  thermocouple, 
(5)  microstrain  transducer  core,  (6)  high  purity  nitrogen,  (7)  gas  flow  meter, 
(8)  copper  coils  immersed  in  liquid  nitrogen,  (9)  heating  elements. 


FIG.  4.  Flexoelectric  polarization  as  a  function  of  applied  strain  gradient 
measured  at  a  flexure  frequency  of  1  Hz. 


Experimental  results  for  the  flexoelectric  polarization  as 
a  function  of  the  transverse  strain  gradient  are  presented  in 
Fig.  4.  The  measurements  were  made  near  the  clamped  end 
of  the  bar  (jc/L=0.25)  at  a  temperature  of  23  °C  where  this 
BST  solid  solution  is  in  the  paraelectric  phase.  As  expected 
the  polarization  is  linearly  proportional  to  the  elastic  strain 
gradient,  but  the  slope  of  the  line  gives  a  magnitude  of  /tI2  of 
100  yuC/m.  This  magnitude  of  /x,12  was  a  major  surprise  as 
the  value  is  not  less,  but  in  fact  some  20  times  larger  than  the 
highest  value  (At12~4  A*C/m)  measured  in  the  PMN  ce¬ 
ramic. 

For  the  BST  ceramic  measurements  of  fil2  as  a  function 
of  temperature  above  Tc  are  shown  in  Fig.  5.  The  break  in 
the  curve  near  35  °C  correlates  well  with  the  departure  from 
Curie-Weiss  behavior.  The  plot  of  fin  versus  relative  dielec¬ 
tric  permittivity  (closely  equivalent  to  susceptibility  x  in 
these  high  K  ceramics)  is  shown  in  Fig.  6  compared  with  the 
expectation  from  Eq.  (1)  with  y—  1.  There  is  a  quasilinear 
range  between  er=2500  and  er—  1 1 000  but  the  measured 
slope  gives  y—  9.3.  The  change  above  sr~  1 1 000  is  not  un¬ 
expected  as  the  departure  from  Curie-Weiss  law  and  the 
enhanced  loss  signal  the  onset  of  an  additional  polarization 


FIG.  5.  Flexoelectric  coefficient  derived  from  plots  such  as  those  in  Fig.  4 
as  a  function  of  temperature  above  the  Curie  point  Tc . 
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Relative  Permittivity 

FIG.  6.  Flexoelectric  coefficient  vs  dielectric  permittivity.  Solid  dots: 
measured  values  on  BST;  dotted  straight  line:  expectation  from  Eq.  (1)  with 
y—  1 ;  and  solid  line  is  plotted  to  highlight  the  quasilinear  range  between 
ar~  2500  and  er—  1 1  000. 

mechanism,  perhaps  due  to  the  survival  of  some  ferroelectric 
domains  above  the  general  Tc .  The  dropoff  of  response  be¬ 
low  sr=2500  is  unexpected  and  not  associated  with  the  in¬ 
adequate  sensitivity  since  PMN  response  levels  below  4 
jitC/m  were  measured  on  the  same  system. 

It  may  be  argued  that  working  close  to  Tc  in  the 
paraelectric  phase  stresses  engendered  by  the  transverse 
drive  may  be  sufficient  to  raise  Tc  and  induce  ferroelectricity. 
In  these  studies  however  the  drive  levels  are  very  low  so  that 
the  maximum  strain  at  the  surface  is  of  order  of  LI 
X  1CT6.  Assuming  values  of  elastic  compliance  and  electros- 
triction  constants  similar  to  Bali03 ,5  the  associated  stress 
would  move  the  Curie  point  —0.03  °C.  With  equal  but  oppo¬ 
site  stresses  at  the  two  surfaces  this  would  produce  an  unde¬ 
tectable  broadening  of  the  response.  That  SrH03  based  ce¬ 
ramics  may  potentially  have  polar  grain  boundaries  is  likely 
as  demonstrated  by  Petzelt6  for  pure  3rTi03  ceramic.  In  the 
volume  of  the  ceramic  one  may  expect  the  Curie  group 
symmetry  of  the  randomly  axed  structure  to  average  the  pi¬ 
ezoeffect  to  zero.  For  the  surface  however,  the  symmetry  is 
lower  and  residual  piezoelectricity  is  possible.  In  the  flexure 
experiment,  the  two  major  surfaces  are  subjected  to  equal 
and  opposite  tractions,  so  the  effect  should  largely  balance 
out.  However  the  two  surfaces  may  not  be  exactly  equiva¬ 
lent,  leaving  a  small  response.  To  test  for  this  possibility  we 
measured  the  phase  of  the  flexoelectric  signal  with  the 
sample  inverted.  Thus  if  the  signal  was  piezoelectric,  the 
phase  should  change  by  180°.  No  change  of  phase  was  mea¬ 
sured,  thus  precluding  the  surface  effect. 

Although  the  stable  narrow  polar  grain  boundary  phase 
will  not  by  symmetry  affect  the  elasto-dielectric  responses, 
as  Petzelt6  suggests,  close  to  Tc  the  surface  polarity  may 
induce  maeropolar  domains  in  some  grains.  Such  domains 
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would  now  contribute  a  new  dielectric  contribution  to  re¬ 
sponse,  strongly  coupled  to  the  elastic  behavior  and  could  be 
the  origin  of  the  departure  from  Curie- Weiss  law  and  the 
enhanced  flexoelectric  response  observed  in  BST  samples 
close  to  Tc . 

There  is  extensive  literature  on  flexoelectricity7"11 
mostly  focused  upon  low  permittivity  simple  ionic  solids  and 
soft  elastomer  systems.  Tagantsev4  appears  to  have  been  one 
of  the  first  to  highlight  the  importance  of  the  dielectric  sus¬ 
ceptibility  in  enhancing  both  surface  and  volume  contribu¬ 
tions,  and  raising  the  magnitude  above  the  accepted  level  of 
eta  which  gives  a  low  value  of  KT^-KT11  C/m.  The 
enhancement  due  to  dielectric  susceptibility  is  of  cardinal 
importance  for  the  high  permittivity  incipient  ferroelectrics 
of  interest  here.  Long  experience  with  electrostriction  in  per- 
ovskites  has  shown  a  remarkable  consistency  in  the  magni¬ 
tudes  of  the  Qijk i  constants.  It  may  be  noted  however,  that  in 
spite  of  excellent  potential  in  transducers,  the  relaxor  ferro¬ 
electrics  like  PMN  have  unusually  low  Qtjkl  values,  particu¬ 
larly  at  temperatures  close  to  Tm ,  the  temperature  of  the 
dielectric  maximum.12  Should  we  expect  similar  consistency 
in  the  normalized  flexoelectric  constant  fiij lx  and  will  the 
relaxors  again  appear  anomalous.  Unfortunately  there  is  cur¬ 
rently  almost  no  database  of  flexoelectric  constants  from 
which  to  judge. 

Recent  measurements  on  the  Frequency  Agile  Materials 
for  Electronics  (FAME)  DARPA  program13  suggest  that  lead 
strontium  titanate  (PST)  solid  solutions  can  be  fabricated 
with  Curie  points  close  to  ambient,  classic  Curie-Weiss  be¬ 
havior,  no  dispersion  and  veiy  low  dielectric  loss,  testifying 
to  the  absence  of  relaxor  effects.  Large  samples  are  difficult 
to  procure  but  will  be  measured.  Similarly,  flexoelectric  mea¬ 
surements  on  a  nonlead  based  relaxor  ferroelectric  such  as 
barium  zirconium  titanate  could  also  give  important  addi¬ 
tional  information. 
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Abstract 

Strain-gradient-induced  polarization  or  flexoelectricity  was  investigated  in  unpoled 
soft  lead  zirconate  titanate  (PZT)  ceramic  whore  the  texture  symmetry  ooco m  forbids 
macro-piezoelectricity.  Even  under  high  strain  gradient  (1  m  *)  the  induced  polarization  is 
small  (1.6  pC/m2)  at  20  °C.  Higher  strain  gradients  induce  ferroelastic  poling  and  an 
additional  extrinsic  contribution  to  the  flexoelectric  coefficient  p12  raising  the  value  from 
0.5  to  2.0  pC/m.  Cooling  through  the  Curie  Point  ( Tc )  under  maximum  stress  (80  MPa) 
where  the  peak  permittivity  (-20,000)  could  raise  yn  to  20  pC/m,  the  equivalent  electric 
field  is  still  only  —1  kV/m,  inadequate  to  achieve  significant  ferroelectric  poling  The 
situation  may  be  different  in  thin  PZT  films  where  much  larger  strain  gradients  can  occur. 
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The  well-known  elasto-electric  coupling  effects  include  piezoelectric  effect, 
electrostrictive  effect,  and  Maxwell  stress  effect.  Piezoelectric  ceramics  [1]  and  more 
recently  single  crystals  [2]  have  demonstrated  the  capability  for  broad  applications  in 
sensor,  actuator  and  transducer  devices.  In  soft  polymers,  it  is  shown  recently  that 
Maxwell  stress  effect  (attractive  forces  between  opposite  charges  on  the  electrodes)  can 
generate  ultrahigh  strain  responses  and  exhibit  great  potential  for  a  variety  of 
electromechanical  device  applications  [3].  All  the  above-mentioned  physical  effects, 
however,  generally  assume  the  situations  of  uniform  stress  or  strain.  In  nature,  there  is 
elasto-electric  coupling  caused  by  inhomogeneous  deformation  where  stress  or  strain 
gradient  associated  polarization  effects  (flexoelectric  effects)  need  to  be  considered. 
Overall,  mechanical  stress  or  strain  can  generate  electric  polarization  in  a  deformable 
dielectric  material  through  the  following  two  mechanisms, 

dsjk 

~^ijk^jk  Mijkl  •  (^) 

In  Eq.  (1)  and  thereafter  Einstein  summation  convention  is  assumed  (i,j,  k,l=  1,  2,  3). 
The  first  term  on  the  right-hand  side  refers  to  the  well-known  direct  piezoelectric  effect, 
where  <7Jk  is  the  stress  uniformly  distributed  across  the  sample  and  dijk  is  the 

piezoelectric  coefficient,  a  third-rank  polar  tensor.  The  second  term  on  the  right-hand  side 

dSjk 

refers  to  the  strain  gradient  ( — — )  induced  polarization  and  juijkI  is  the  flexoelectric 

dxt 

coefficient,  a  fourth-rank  polar  tensor.  In  centrosymmetric  materials,  dijk  =  0 ,  so  the 

piezoelectric  term  in  Eq.  (1)  can  be  eliminated,  therefore, 
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where  Pt  is  the  electric  polarization  induced  solely  by  strain  gradient. 

When  reviewing  the  history  of  flexoelectric  investigations,  it  is  noted  that  while  the 
concept  was  originally  formed  in  early  1960s  [4,  5],  not  until  1981  was  such  effect  in 
crystalline  solids  given  the  name  “flexoelectric”  [6].  Based  upon  an  ionic  model  [7], 
Tagantsev  analyzed  the  flexoelectric  effect  and  suggested  possible  larger  effects  in 
ferroelectrics.  In  soft  polymers  Marvan  et  al.  [8]  observed  flexoelectric  coefficients  of  the 
order  of  10"11  to  10'10  C/m.  Up  to  now  very  little  attention  has  been  paid  to  test  the 
magnitude  of  flexoelectric  coefficients  and  the  mechanism  remains  unclear.  Recently  we 
measured  the  flexoelectric  coefficients  (jil2)  in  lead  magnesium  niobate  (PMN)  [9,  10]  (a 
well-known  relaxor  ferroelectric  material)  and  barium  strontium  titanate  (BST)  [11]  (a 
normal  ferroelectric  material).  Both  materials  were  tested  in  the  phase  region  with 
macroscopic  cubic  symmetry,  moreover  the  samples  were  measured  in  the  form  of  a 
cantilevered  beam  so  that  any  remnant  piezoelectric  contributions  from  the  top  and 
bottom  halves  of  the  beam  would  cancel.  In  this  letter,  we  investigate  the  flexoelectric 
effect  in  a  well-known  PZT  piezoelectric  ceramic  in  the  ferroelectric  phase  by  using  a 
four-point  bend  fixture  to  generate  a  uniform  strain  gradient. 

La  this  work  several  interesting  questions  were  explored:  (i)  Using  a  four  point  bending 
fixture  is  it  possible  to  induce  and  measure  a  quasi-static  flexoelectric  polarization 
generated  by  y.n\  (ii)  Is  the  flexoelectric  polarization  enhanced  or  inhibited  by  the  onset 
of  ferroelastic  domain  wall  motion  which  will  be  evidenced  by  change  of  the  flexural 
stiffness  and  the  development  of  remnant  curvature  in  the  sample;  (iii)  Is  it  possible  in 
soft  PZT  to  reach  levels  of  flexoelectric  induced  field  sufficient  to  pole  the  ceramic  into  a 
piezoelectric  form. 
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The  samples  used  were  unpoled  PZT-5H  ceramics  (doped  with  La  and  Sn)  fabricated 
by  TRS  ceramics  company,  State  College,  Pennsylvania.  Dielectric  measurements 
performed  using  an  HP4284A  LCR  meter  show  a  weak  field  permittivity  of  2,200  at  20 
°C  and  a  strong  but  rounded  dielectric  maximum  (f3  =20,000)  without  little  dispersion 

over  the  frequency  range  of  1  kHz  to  1MHz  [Fig.  1(a)],  suggesting  a  diffuse  phase 
transition  without  strong  relaxor  character.  Polarization  hysteresis  loops  measured  using  a 
modified  Sawyer-Tower  circuit  show  remnant  polarization  of  35  pC/cm2  and  coercive 
field  of  7  kV/cm  [Fig.  1(b)]  with  no  discemable  bias.  Young’s  modulus  was  measured  to 
be  70  GPa  by  a  dynamic  resonance  method. 

A  uniform  strain  gradient  was  generated  using  a  4-point  bend  fixture  schematically 
illustrated  in  Fig.2(a).  Samples  for  measurement  had  dimensions  60  mm  long,  7  mm  wide 
and  3  mm  thick.  Surfaces  were  carefully  polished  and  the  samples  were  annealed  at 
700°C  to  relieve  surface  stresses.  Sputtered  gold  electrodes  with  dimensions  10x7  mm2 
were  applied  to  upper  and  lower  surfaces  in  the  center  of  the  beam.  No  residual 
piezoelectricity  could  be  detected  by  Berlincourt  d33  meter,  and  the  impedance  trajectory 
from  1  kHz  to  1MHz  was  free  from  evidence  of  piezoelectric  resonance.  The  bend  test 
was  carried  out  according  to  the  ASTM  C-l  161-94  by  using  an  Instron  machine  Model 
4202  with  a  10  kN  load  cell.  Outer  and  inner  spans  of  the  fixture  are  40  mm  and  20  mm 
respectively.  The  generated  electric  charge  was  detected  by  a  Keithley  6517  electrometer 
that  can  resolve  10  fC  and  will  measure  up  to  2.1  pC.  Before  measurement  the 
electrometer  was  carefully  calibrated  for  voltage  burden  and  input  offset  current. 

The  stress  distribution  along  the  sample  length  direction  (xi)  is  shown  in  Fig.2(b). 
Within  the  inner  span  the  stress  <7n(x,)  is  uniform  along  the  length,  while  along  the 
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thickness  direction  (x3)  the  stress  <rn(x3)  varies  and  there  is  a  stress  or  strain  gradient  as 

shown  by  Fig.  2(c),  which  can  generate  electric  polarization  through  the  flexoelectric 
effect.  Because  the  bar  length  is  much  greater  than  the  bar  thickness,  we  omit  the  shear 
stress  and  strain  and  only  consider  the  principal  stress  and  strain.  Therefore,  for 
simplification  here  only  one  suffix  was  used  for  describing  the  stress  and  strain  tensors. 
The  absolute  value  of  surface  stress  was  calculated  using  the  following  equation 


<7i(x3) 


|  -  3FL 
*s=±§  4  wd2 


(3) 


where  F  is  the  load,  w  the  width,  d  the  thickness,  and  L  the  outer  span  of  the  bend  fixture. 
The  strain  gradient  in  the  thickness  direction  is  given  by 

dsx  (x3)  12  st 
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where  s  is  the  crosshead  speed  and  t  the  time  gone  by. 

Figure  3  presents  data  on  the  surface  stress  vs.  strain  curves  for  tests  carried  out  at  1 
mm/min  and  0.2  mm/min  crosshead  speed.  As  shown  in  Fig.3(a),  the  flexure  strength  is 
measured  to  be  85  MPa.  The  softening  of  the  samples  at  certain  stress/speed  levels  (e.g., 
around  30  MPa  for  1  mm/min  crosshead  speed)  we  believe  corresponds  to  the  onset  of 
ferroelastie  domain  motion.  As  expected  in  PZT  the  sample  becomes  softer  at  lower 
loading  rates  due  to  the  relaxational  nature  of  ferroelastie  domain  switching.  That 
ferroelastie  switching  has  occurred  at  the  higher  loading  levels  is  evident  by  a  static 
remnant  curvature  in  the  sample  after  testing  (Fig.4). 

For  an  unpoled  PZT  ferroelectric  ceramic  at  the  morphotropic  phase  boundary, 
although  individual  grains  may  have  lower  symmetry  (tetragonal  or  rombohedral)  which 
permit  piezoelectricity,  in  the  volume  of  the  ceramic  one  may  expect  a  macroscopic 
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symmetry  of  oooo m ,  so  the  non-zero  components  for  the  flexoelectric  coefficient  //,y# 
should  be //nil,  ^1122,  and// 1212,  or  in  matrix  notation  // n,  // 12  and  /M-  Since  in  the  current 
investigation  any  remnant  piezoelectric  effects  from  the  top  and  bottom  halves  are  well 
balanced  and  there  is  only  strain  gradient  in  the  thickness  or  x?  direction,  so  the  generated 
electric  polarization  can  be  composed  of  only  flexoelectric  polarization  as  shown  below, 


PZ  (*3  )<*]  +  Mn  d-e'-Xl>  =  ft, 2 


Flexoelectric  polarization  vs.  strain  gradient  is  displayed  in  Fig.5.  It  is  clear  that  the 
behavior  is  not  linear,  showing  a  low-gradient  and  high-gradient  linear  behavior  but  of 
different  slope.  The  changeover  occurs  at  the  gradient  ~  0.3  m'1  corresponding  to  the 
onset  of  ferroelastic  switching.  At  very  small  strain  level,  the  sample  is  only  subjected  to 
elastic  deformation,  so  the  measured  electric  polarization  response  represents  the  intrinsic 
flexoelectric  effect.  The  low-gradient  slope  corresponds  to  a  flexoelectric  coefficient  //12 
~  0.5  pC/m,  and  the  high-gradient  slope  takes  //12  ~  2.0  pC/m  showing  that  ferroelastic 
domain  wall  motion  aids  the  response. 

The  measured  flexoelectric  polarization  is  very  small  compared  to  the  remnant 
polarization,  only  1.6  pC/m2  at  a  strain  gradient  of  1  m'1.  By  using  a  relative  permittivity 
value  of  er  =2,200  (at  1  kHz),  we  figure  out  that  a  strain  gradient  of  1  m‘l  is  equivalent  to 
an  electric  field  of  100  V/m,  which  is  obviously  too  small  compared  to  the  coercive  field. 
On  a  simple  point  charge  model,  flexoelectric  coefficients  has  the  magnitude  el  a 

where  e  is  the  electron  charge  and  a  the  dimension  of  the  unit  cell  [5].  Thus  in  high- 
permittivity  ferroelectrics  we  may  expect  from  elementary  theory  of  flexoelectricity  [7] 
that  /Uy  be  proportional  to  dielectric  susceptibility  Xy  following  a  relation  [11]: 
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(6) 


M9  ~  YXtj  ~ 
a 

where  y  is  a  constant  of  value  close  to  unity.  For  the  PZT  unpoled  ceramic  at  low 
gradient  levels  the  normalized  flexoelectric  coefficients  %n  at  20  °C  is  0.23  nC/m 
giving  7  « 0.57.  In  lead  magnesium  niobate  (PMN)  from  our  earlier  study  [9,  10], 
/^i2  •  Z 22  is  0.26  nC/m  giving  y  ~  0.65,  both  in  reasonable  accord  with  the  elementary 
theory.  For  Bao^Sro^TiCb  (BST)  ceramic  at  25°C  however //i2/j22 -3.72  nC/m  yielding 
a  value  y  =9.3  [1 1]  much  higher  than  those  in  the  lead-based  systems. 

It  may  be  noted  that  in  earlier  investigations  of  the  thermo-polarization  effects  Strukov 
et.  al.  [13]  found  exceedingly  high  values  for  the  normalized  thermopolarization 
coefficient  (bjj)  in  triglycine  sulfate  (TGS)  which  they  attributed  to  the  order-disorder 

nature  of  the  ferroelectric  phase  change  in  this  compound.  We  note  that  in  BaxSri_xTi03 
[13],  there  is  strong  evidence  of  a  local  order-disorder  component  between  polarization 
vectors  in  the  unit  cell.  It  will  be  interesting  to  measure  juj2  in  potassium  niobate  tantate 
(KTN),  where  again  order-disorder  has  been  identified  [14],  to  see  if  the  coefficients  are 
again  anomalously  large. 

In  conclusion,  the  PZT  ceramic  does  exhibit  modest  flexoelectricity  and  ferroelastic 
domain  wall  motion  enhances  the  response,  but  it  is  not  possible  in  our  samples  using 
stress  levels  up  to  the  full  fracture  strength  to  induce  ferroelectric  poling,  although 
ferroelastic  poling  was  patently  obvious. 

In  the  PZT  thin  films  epitaxially  grown  on  lattice-mismatched  substrates  the  strain 
gradient  between  surfaces  and  strain-relieving  dislocations  can  be  exceedingly  large  and 
the  flexoelectric  effects  could  be  of  major  importance  in  these  systems. 
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Figure  Captions 


Fig.  1  (a)  Weak  field  dielectric  permittivity  (E  ~  10  V/cm)  as  a  function  of  frequency  and 
temperature  in  the  soft  PZT-5H  sample  measured;  (b)  dielectric  hysteresis  (P  vs.  E) 
in  the  soft  PZT-5H  composition. 

Fig-2  Schematic  illustration  of  strain  gradient  induced  quasi-static  polarization 
measurement,  (a)  typical  4-point  bending  fixture;  (b)  stress  distribution; 

(c)  schematic  of  strain  gradient  and  the  induced  polarization  along  the  thickness 
of  the  sample. 

Fig.3  Surface  stress  as  a  function  of  strain  measured  at  two  different  crosshead  speeds  on 
PZT-5H  bar,  (a)  1  mm/min  cross  head  speed  (the  arrow  indicates  the  level  of  stress 
at  which  the  sample  broke);  (b)  0.2  mm/min  crosshead  speed. 

Fig-4  Obvious  static  curvature  induced  in  the  PZT-5H  sample  by  ferroelastic  domain  wall 
motion  at  the  higher  stress  levels. 

Fig.5  Polarization  versus  strain  gradient  for  an  unpoled  PZT-5H  sample  induced  during  a 
4-point  bend  test  carried  out  at  a  crosshead  speed  of  0.2  mm/min 
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List  of  the  changes  on  Ms.  #L02-3586 


1 .  The  discussion  of  physical  nature  for  the  experimental  result  is  added,  please  see  the 
1st  and  2nd  paragraphs  in  Page  7.  It  appears  that  the  lead-based  systems  (PMN,  PZT) 
are  in  reasonable  accord  with  the  elementary  theoiy,  while  in  BST  the  coefficients  are 
exceedingly  large.  The  order-disorder  components  are  supposed  to  be  related  to  this 
phenomenon,  but  future  efforts  are  needed  to  clarify  it. 

2.  Ref.  1 1  was  completed,  ref.  12  was  changed,  and  refs.  13  and  14  were  added. 

3.  Dielectric  measurements  and  ferroelectric  hysteresis  loop  were  added  and  listed  as 
Fig.  1(a)  and  Fig.  1(b),  respectively. 

A  picture  showing  the  remnant  curvature  after  4-point  bend  test  was  added  and  listed 
as  Fig.4. 

4.  Eqs  4, 5  and  9  in  the  old  manuscript  were  removed  in  order  to  make  the  paper  more 
concise. 

5.  In  the  abstract,  the  last  paragraph  in  Page  3,  and  the  last  paragraph  in  Page  7,  the 
significance  of  the  investigations  on  flexoelectric  effect  were  further  elaborated  and  it 
is  suggested  that  in  thin  films  the  effect  can  be  of  major  importance. 
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The  macroscopic  model  proposed  fay  Kittel  for  antiferroelectrics  assumes  two  interpenetrat¬ 
ing  sublattices  with  antiparallel  polarizations.  Although  the  model  can  explain  many  proper¬ 
ties,  it  does  not  provide  microscopic  lattice  displacement  patterns  associated  with  these 
sublattices  since  the  order  parameter  in  the  model  is  a  volume-averaged  quantity.  We  propose 
a  model  based  on  microscopic  positioning  of  symmetry  elements  and  use  an  order  parameter 
field  directly  associated  with  the  local  dipoles,  which  can  provide  a  more  transparent  defini¬ 
tion  for  the  sublattices  used  in  the  Kittel  model.  We  also  include  nonlocal  gradient  interac¬ 
tions  to  model  multidomain  structures  and  domain  walls  in  ferroeieetrics. 

Keywords:  antiferroelectric;  Landau  theory  of  phase  transitions;  sublattice  model;  micro¬ 
scopic  symmetry 

INTRODUCTION 

The  antiferroelectric  (AFE)  state  was  modeled  by  Kittel  based  on  die 
Landau  theory  and  use  two  interpenetrating  sublattices  with  equal  but 
opposite  polarizations.1  This  model  can  successful  describe  the  double 

A 

hysteresis  loops  observed  in  antiferroelectrics,  but  contains  no 
mechanism  to  describe  the  spatial  relationship  between  the  two 
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sublattices.  In  other  words,  there  is  no  mechanism  to  guarantee  the  two 
sublattices  to  exist  alternately  in  space  to  form  the  antiferroelectric 
state. 

We  propose  a  microscopic  model  for  the  AFE  state,  which  not  only 
provides  the  origin  of  the  Kittel  sublattices,  but  also  gives  the  fixed 
spatial  relationship  between  them.  This  model  is  based  on 
microscopic  symmetry  and  group  theoretical  methods.  By  adding  the 
gradient  energy  associated  with  the  AFE  state,  it  can  also  be  used  to 
describe  inhomogeneous  structures,  i.e.,  orientation  twins  and 
antiphase  configurations. 

In  describing  this  model,  we  will  use  ammonium  dihydrogen 
phosphate,  or  ADP,  as  a  prototype.  The  method  itself  is  quite  general, 
and  can  be  applied  to  any  other  antiferroelectric  species. 


MICROSCOPIC  SYMMETRY 


The  high-temperature  paraelectric  phase3"6  of  ADP  is  tetragonal,  with 
space  group  I42d .  At  Tc  *  -125°C,.a  transition  to  an  orthorhombic 
AFE  phase4"7  with  space  group  P2j  2j  2\  occurs,  with  a  doubling  of  the 
primitive  unit  cell. 


Z=1/4^  -o 


z  =  1/2  + 
2  =  3/4 


Figure  1  Dipole  moments  within  the  AFE  unit  cell. 

The  transition  is  driven  by  an  M-point  soft  mode  corresponding  to 
the  M3M4  physically  irreducible  representation.8  The  antiferroelectric 
state  consists  of  four  individual  molecular  dipole  moments,  parallel  to 
the  a-b  plane,  within  the  unit  cell  of  the  AFE  phase  as  shown  in  Fig.  1. 
The  transition  can  be  described  by  a  two-component  order  parameter 
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Figure  2  Directions  of  (000)  dipole  moment  in  the  four  AFE  domains. 

that  corresponds  to  the  components  of  the  molecular  dipole  at  the 
(000)  position.  There  are  four  energetically  equivalent  orientations  for 
this  dipole,  as  shown  in  Figure  2,  and  so  there  are  four  distinct 
domains  in  the  AFE  phase,  with  the  two-component  order  parameters 
as  listed  in  Table  I.  For  convenience,  we  assume  pa>Pb  in  the  low 
temperature  orthorhombic  structure  with  space  group  P2j2]  2;. 


TABLE  I.  Values  of  the  order  parameter  for  the  four  different  domain 
states. 


Domain 

Order  Parameter 

Si 

(p*Pb ) 

s2 

(Pk-Pa) 

Si 

(-Pc-Pb) 

s. 

(" Pb,  Pa) 

Since  the  space  groups  of  the  two  phases  and  the  irreducible 
representation  are  known,  the  symmetry  allowed  distortions,  based  on 
the  positions  of  the  microscopic  symmetry  elements,  can  be 
calculated.  This  is  done  using  die  ISOTROPY  computer  program9,  and 
the  results  for  the  Wyckoff  a-sites  are  shown  in  Table  II.  Using  the 
distortions  in  this  table,  the  dipole  arrangement  can  be  reproduced  for 
any  single  domain  state  given  in  Table  I  (similar  to  the  arrangement 
shown  in  Fig.  1  for  Si).  The  relationship  between  these  allowed 
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distortions  and  the  positions  of  the  microscopic  symmetry  operations 
is  illustrated  for  domain  Si  in  Fig.  3 


TABLE  n.  Group-theoretically  allowed  distortions  at  Wyckoff  fi¬ 
shes  in  domain  Si. 


Position 

Distortions 

(0,0,0) 

pa  (1,0,0)  +pb  (0,-1, 0) 

Pa  (-1,0,0)  +pb  (0,1,0) 

(0>~2>i) 

Pa  (-1,0,0)  +p*(0,-l,0) 

Pa  (1,0,0)  +Pb  (0,1,0) 

RELATIONSHIP  TO  THE  KITTEL  SUBLATTICES 

Our  model  has  focused  on  the  microscopic  structure  of  the  AFE  state. 
In  order  to  understand  the  macroscopic  behavior  (in  particular,  the 
double  hysteresis  loops  observed  when  the  antiferroelectric  material  is 
placed  in  an  external  electric  field),  we  must  make  a  connection 


fig.  3  Actions  of  symmetry  elements  generate  dipole  arrangement 
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between  our  microscopic  model  and  the  Kittel  macroscopic  model. 
Since  polarization  is  an  average  of  molecular  dipole  moments  over  a 
given  volume,  the  sublattice  polarizations  Pi  and  P2  in  the  one- 
dimensional  Kittel  model  can  be  defined,  in  terms  of  the  molecular 
dipole  moments,  as 

<P.*»  (pj< 0) 

In  other  words,  in  the  AFE  state  those  dipole  moments  with  a  positive 
component  in  the  a-direction  form  one  sublattice  with  polarization  Pj, 
mid  those  dipole  moments  with  a  negative  component  in  the  a- 
direction  form  the  other  sublattice  with  polarization  P2,  as  shown  in 
Fig.  4. 

On  a  unit  cell  level,  we  can  further  define  the  sublattices  in  terms  of 
symmetry  operations.  Referring  to  die  positioning  of  symmetry 
operations  in  Fig.  3,  define  p+  and  p.  as 

P+=P  +  {Cu\j[0i}p 

P-=fC2lliii)p +{C2,|0ii)p  (2a’b) 


Figure  4  Defining  the  macroscopic  Kittel  sublattices  in  terms  of  the 
microscopic  dipole  moments. 


134/[346] 


R.A.  HATT  and  W.  CAO 


These  two  quantities,  when  volume  averaged  Over  all  unit  cells, 
naturally  recreate  the  Kittel  sublattices  on  the  microscopic  scale. 

Our  microscopic  theory  presented  here  not  only  can  derive  the 
macroscopic  Kittel  model,  but  also  can  account  for  the  spatial 
relationships  among  dipoles  of  adjacent  cells.  More  importantly,  it 
also  gives  the  dipole  tilt,  which  occurs  in  the  ADP  system,  leading  to 
the  observed  unit  cell  doubling  associated  with  the  antiferroelectric 
phase  transition. 


FREE  ENERGY  DENSITY 


The  free  energy  originally  incorporated  in  the  Kittel  model  included 
only  local  terms: 

F  =  aj(P;  +  p|j+a2  PjP2+cc3( Pl+ P42)  (3) 


In  our  model,  however,  knowledge  of  the  irreducible  representation 
involved  in  the  transition  allows  non-local  terms  to  be  added  in  the 
form  of  gradient  invariants.  For  the  ADP  case  being  considered,  this 
free  energy  density  can  be  written  as 


F  =  A(p?  f  Pj)  +  B,(P.J  +  p\)7  +  Bt(pt  +  Pi) 

+  Bj(pIPi  -  P,pI)+Ci(P?  +  pl)3  +  Cjipfpl  +  pIp\) 

+  C3(p?p,  -  P, )  +  Dx 

+  D2 
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where  pi  and  P2  refer  to  the  two  components  of  the  order  parameter, 
i.e.,  the  a-  and  b-  components  of  the  (000)  molecular  dipole  moment. 
All  coefficients  are  temperature  independent  except A  =  A^(T-T0). 

With  the  addition  of  the  order  parameter  gradient  terms,  this  free 
energy  can  be  used  to  describe  the  change  in  order  parameter  across  a 
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domain  wall  in  an  inhomogeneous  structure.  For  example,  the  change 
in  order  parameter  across  the  domain  wall  in  an  orientation  twin 
formed  by  domains  Si  and  S2  can  be  calculated  by  normalizing  the 
order  parameter  using  the  following  change  of  variables: 


Pi  =  P?  fi  +  PTfi  =  pT  (ft  +Kfi)  (5a) 

Pi  =  pTfi-pfh  =  Pfitfi  -/2).  (5b) 

This  leads  to  a  coupled  set  of  differential  equations  upon  energy 
minimization, 

+/,s  +  (6a) 

fu v— T/a  ~ft  +ft  +Wft  +2/.V,3),  (6b) 


with  the  boundary  condition 


/i=jvi+Vi+T »  x 


t*  *  , 

x  — >+<x> 


.  /*- 


// 

X  -»-oo 


ft+Ji+T ,  *"-»+«= 


where  t  describes  the  temperature  dependence.  These  aquations  can 
be  solved  numerically,  and  solutions  for  three  temperatures  are  shown 
in  Fig.  5.  Further  details  of  this  process  are  reported  elsewhere.10*11 


*3-2-10123 
rescaled  x"  axis 


Fig.  5  Calculated  profiles  for  a  representative  orientation  twin 
for  three  different  temperatures. 
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CONCLUSION 

We  proposed  a  Landau-Ginzburg  type  model  for  the  antiferroelectric 
state  based  on  microscopic  symmetry  and  group  theory.  The  order 
parameter  corresponds  to  the  local  dipole  moment  within  the  unit  cell, 
instead  of  the  macroscopically  averaged  polarization  conventionally 
used.  Our  model  can  rigorously  derive  the  sublattices  of  the  model 
proposed  by  Kittel  and  address  the  dipole  tilt  and  cell  doubling  in  the 
antiferroelectric  transition.  The  so  derived  sublattices  are  intrinsically 
linked,  which  resolved  the  confusion  encountered  in  the  Kittle  model. 

The  model  is  also  extended  to  include  order  parameter  gradient 
terms.  These  nonlocal  interaction  terms  allow  inhomogeneous 
structures,  such  as,  orientation  twins  and  antiphase  walls,  to  be 
modeled. 
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The  only  Landau-type  model  for  antiferroelectric  phase  transitions  was  proposed  by  Kittel,  in  which  two 
interpenetrating  sublattices  with  opposite  polarizations  of  equal  amplitude  were  assumed.  The  theory,  however, 
did  not  include  any  mechanism  to  specify  the  relative  spatial  positions  of  the  two  sublattices,  and  therefore 
could  not  address  the  cell  doubling  during  antiferroelectric  phase  transitions.  We  propose  a  Landau-Ginzburg- 
type  model  based  on  microscopic  symmetry  and  group  theory,  which  can,  without  having  to  assume  sublat¬ 
tices,  account  for  all  aspects  of  antiferroelectric  states,  including  local  dipole  orientation  and  cell  doubling.  The 
average  of  these  dipoles  naturally  leads  to  the  Kittel  model.  The  inclusion  of  gradient  terms  in  the  free  energy 
allows  the  modeling  of  multidomain  structures  and  domain  walls  in  antiferroelectric  states. 


INTRODUCTION 

Kittel1  proposed  a  macroscopic  Landau-type  model  for 
the  antiferroelectric  (AFE)  state  by  introducing  two  interpen¬ 
etrating  sublattices  with  opposite  polarizations.  His  model 
describes  a  second-order  transition  from  the  paraeleetric  state 
to  the  antiferroelectric  state,  by  truncating  the  free  energy  at 
the  fourth  order: 

F=a1(P2l+Pl)+a2P1P2+an(Pi+P42),  (1) 

where  Pt  and  P2  are  the  polarizations  of  the  sublattices.  If 
a2>0,  the  transition  will  favor  P1  and  P2  being  antiparallel, 
making  the  low-temperature  phase  antiferroelectric.  On  the 
other  hand,  if  a2< 0,  the  transition  will  favor  Pl  and  P2 
being  parallel,  and  the  transition  will  lead  to  a  ferroelectric 
state.1”2 

This  model  has  the  intrinsic  limitation  that  it  contains 
only  local  interactions.  In  other  words,  there  is  no  mecha¬ 
nism  to  fix  the  spatial  relationship  between  the  two  sublattice 
polarizations  within  the  crystal.  This  local  model  creates  un¬ 
certainty  in  the  antiferroelectric  state.  For  example,  the  sub¬ 
lattice  polarizations  Px  and  P2  are  assumed  at  the  same  lo¬ 
cation  in  space  (or  can  be  anywhere  in  space),  which  leads  to 
the  cancellation  of  Px  and  P2.  Such  a  situation  does  not 
fully  describe  the  antiferroelectric  state  in  which  adjacent 
primitive  cells  acquire  opposite  dipole  moments  and  the  lo¬ 
cal  polarization  at  any  space  point  is  actually  nonzero.  Using 
group-theoretical  techniques,  we  propose  a  continuum  model 
for  the  antiferroelectric  state  built  upon  microscopic  symme¬ 
try.  The  symmetry  allowed  distortions  associated  with  the 
soft  mode  are  given  in  our  model,  which  accounts  for  the 
formation  of  antiparallel  dipoles  in  adjacent  cells.  The  free 
energy  has  also  been  expanded  to  include  gradient  terms  of 
the  order  parameter  so  that  multidomain  structures,  such  as 
orientation  twins  and  antiphase  walls,  can  be  modeled  using 
the  same  formulation.  Specifically,  in  this  paper  we  will  use 
ammonium  dihydrogen  phosphate  (NH4)H2(P04)  (commonly 


referred  to  as  ADP)  as  a  prototype  system  to  illustrate  the 
procedure  and  the  characteristics  of  the  antiferroelectric 
phase  transition. 


ANTIFERROELECTRIC  PHASE  TRANSITION 

The  high-temperature  phase  of  ADP  is  tetragonal  with 
space  group  /42<L3"6  Following  the  notation  of  the  Interna¬ 
tional  Tables  for  Crystallography ,7  the  conventional  (non¬ 
primitive)  unit  cell  contains  four  formula  units.  Two  of  the 
formula  units  are  in  the  primitive  cell  while  the  other  two  are 
related  by  the  centering  translation  (555).  At  rc*= 
~  125  °C,  the  material  undergoes  an  antiferroelectric  transi¬ 
tion  to  an  orthorhombic  phase  with  space  group  P232121 ,6,8 
The  transition  is  driven  by  a  zone-boundary  M-point  soft 
mode  and  results  in  ionic  displacements  that  create  one  di¬ 
pole  moment  in  each  formula  unit  (four  dipoles  per  AFE  unit 
cell)  as  shown  in  Fig.  1.  These  dipole  moments  form  a  net 
antiferroelectric  polarization  along  the  [100]  or  [010]  of  the 
parent  phase.  This  means  that  the  fourfold  rotational  inver¬ 
sion  axis  of  the  parent  paraeleetric  phase  is  lost  during  the 
transition,  which  leads  to  two  rotationally  related,  energeti¬ 
cally  equivalent  orientation  domain  states  in  the  low- 
temperature  phase.  Moreover,  the  transition  causes  a 
primitive-cell  doubling  (equivalently,  the  centering  point  in 
the  conventional  parent  cell  is  lost),  leading  to  two  additional 


z=1/4^ — q 

z  =  3/4  a  — 


FIG.  1.  Dipoles  formed  in  the  antiferroelectric  phase  of  ADP  in 
a  conventional  cell.  The  lattice  displacement  pattern  can  be  gener¬ 
ated  from  group  theory. 
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2  =  0 


z=1/4'*"'-o 


Z=V2-±-'* 


2  =  3/4 


FIG.  2.  Microscopic  positions  of  the  symmetry  elements,  creat¬ 
ing  the  antiferroelectric  polarization.  The  AFE  unit  cell  is  shown  by 
the  dashed  line,  with  origin  at  (0,  J,  f).  The  screw  axes  parallel  to  a 
are  at  z=  |,  and  the  screw  axes  parallel  to  b  are  at  z=  4 . 


antiphase  states  that  are  translationally  related  to  the  two 
orientation  states.  Overall,  there  are  four  possible  domains  in 
the  AFE  state  of  ADP. 


GROUP  THEORETICAL  DESCRIPTION  OF  THE 
TRANSITION 

Once  the  structures  and  space  groups  of  the  high-  and 
low-temperature  phases  have  been  determined,  the  micro¬ 
scopic  positions  of  the  symmetry  elements  are  fixed  (we  will 
use  the  settings  of  Ref.  7),  as  shown  in  Fig.  2.  Specifically, 
we  note  that  the  three  mutually  perpendicular  twofold  screw 
axes  in  the  AFE  phase  do  not  allow  a  net  dipole  moment  for 
the  unit  cell  of  the  low-temperature  phase,  although  they  do 
allow  antiferrolectric  polarizations. 

The  M-point  soft  mode  driving  the  transition  in  ADP  may 
be  described  by  the  physically  irreducible  M3Af4 
representation,4  which  is  simply  the  direct  sum 
(here  we  use  the  labeling  of  Miller  and  Love9).  This  repre¬ 
sentation  carries  a  two-component  order  parameter  (OP), 
(p  1  ,p2)>  whose  values  are  listed  in  Table  I.  Physically,  the 
OP  corresponds  to  the  molecular  dipole  moments  within  the 
a-b  plane;  px  represents  the  component  along  the  a  direction, 
and  p2  represents  the  component  along  the  b  direction.  This 
OP  is  a  continuum  field,  which  is  equal  to  the  dipole  moment 
at  the  sites  of  each  formula  unit. 

Knowing  the  irreducible  representation  allows  the  lattice 
distortions  that  arise  in  the  transition  to  be  calculated  using 
the  ISOTROPY10  software  package.  For  simplicity  we  will  re¬ 
strict  our  discussion  to  distortions  arising  at  Wyckoff  a  sites 
only,  which  are  listed  in  Table  II.  These  distortions  can  be 
used  to  construct  the  dipole  arrangement  within  the  low- 
temperature  unit  cell.  As  an  example,  consider  the  distortions 
that  arise  in  domain  Slf  for  which  the  OP  is  (/?j,p2) 
=  (pa>Pb)-  ^ Pa>Pb>  ^e  distortions  result  in  the  configu- 


TABLE  I.  Values  of  the  order  parameter  for  different  domain 
states. 


Domain 

Order  parameter 

Si 

(pfl  -Pi,) 

s2 

(Pb’-Pa) 

S3 

(~ Pa~Pb ) 

S4 

(“Pi.Pa) 

TABLE  II.  Group  theoretically  allowed  distortions  at  Wyckoff  a 
sites  in  domain  Sx . 


Position 

Distortions 

Relationship  to  (0,0,0)  site 

(0,0,0) 

p„(  1,0,0) +p4(0,- 1,0) 

{£|000} 

(Hj) 

Pa(~  l,0,0)+p4(0,l,0) 

{^2zl  2  2  l} 

(0.U) 

Pa(~  1.0,0) +P»(0,-  1,0) 

(i.l.l) 

pfl(  1,0,0)  +/>„(0,1,0) 

{^-2x1  2  1  1} 

ration  shown  in  Fig.  1,  which  agree  well  with  the  observed 
dipolar  distribution  by  Blinc  et  al} 

Figure  1  was  constructed  based  only  on  the  microscopic 
positioning  of  the  symmetry  elements  and  the  allowed  dis¬ 
tortions  from  group  theoretical  considerations,  as  in  Fig.  2. 
For  example,  a  dipole  moment  at  position  (0,  0,  0)  necessar¬ 
ily  means,  by  action  of  the  screw  axis  parallel  to  the  z  axis, 
an  oppositely  directed  dipole  moment  at  the  parent  cell  cen¬ 
tering  point  of  (555).  The  screw  axis  parallel  to  the  x  axis 
further  implies,  from  these  two  dipoles,  that  the  other  two 
antialigned  dipoles  within  the  AFE  unit  cell  must  be  present. 
Hence,  there  is  no  need  for  defining  sublattices  and  no  need 
to  explain  why  such  sublattices  would  adopt  exactly  equal 
magnitudes  but  opposite  orientation.  More  importantly,  the 
spatial  relationship  between  the  two  sublattices  naturally 
comes  out  of  our  model.  If  a  single  dipole  moment  is  formed 
within  a  parent  formula  unit,  one  can  generate  its  counter¬ 
parts  by  symmetry  requirements  to  form  the  antiferroelectric 
polarization  of  the  whole  unit  cell  of  the  low-temperature 
phase.  Also,  since  the  dipole  moment  at  the  centering  point 
(555)  is  opposite  to  the  dipole  moment  at  (0,0,0),  the  dou¬ 
bling  of  the  unit  cell  size  is  naturally  explained. 

HYSTERESIS  AND  THE  RELATION  TO  THE 
SUBLATTICES  PROPOSED  BY  KITTEL 

A  well-known  characteristic  of  antiferroelectrics  is  the 
double  hysteresis  loop,  as  shown  in  Fig.  3.  In  the  one¬ 
dimensional  antiferroelectric  model  proposed  by  Kittel,  there 
is  no  net  polarization.  However,  if  an  external  electric  field  is 
applied,  die  sublattice  polarization  parallel  to  the  field  grows 
and  the  other  sublattice  polarization  opposite  to  the  field 
shrinks,  resulting  in  a  net  polarization.  When  the  field 
strength  becomes  sufficiently  large,  the  polarization  in  the 
direction  opposite  to  the  field  abruptly  switches  orientation 
to  become  parallel  to  the  field,  resulting  in  a  ferroelectric 
state.  Because  of  symmetry  constraints,  a  double  hysteresis 


FIG.  3.  Double  hysteresis  loop,  a  characteristic  of  the  antiferro¬ 
electric  state. 
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FIG.  4.  Molecular  dipole  moments  in  the  presence  of  an  external 
electric  field,  (a)  before  switching  and  (b)  after  switching. 


will  be  produced  during  such  a  switching  process.2 

In  our  model,  we  describe  the  hysteresis  in  terms  of  the 
dipole  moments  associated  with  each  of  the  four  formula 
units  in  an  AFE  unit  cell.  Consider  an  external  electric  field 
applied  to  domain  Sx ,  parallel  to  the  -ha  direction.  This  field 
will  interact  with  the  dipoles,  causing  pa> 0  to  increase  in 
magnitude  and  pa< 0  to  decrease  in  magnitude,  as  shown  in 
Fig.  4(a).  Note  that  pb  is  not  affected  in  our  model.  At  suf¬ 
ficiently  large  field  strength,  the  dipoles  with  pa  antiparallel 
to  the  external  field  will  flip  so  that  pa  becomes  parallel  to 
4-a,  as  shown  in  Fig.  4(b),  This  results  in  a  state  in  which  the 
dipole  components  are  parallel  along  the  a  direction  and  an¬ 
tiparallel  along  the  b  direction.  It  is  a  ferroelectric  state  with 
die  polarization  along  a,  but  with  a  cell  size  twice  that  of  the 
Kittel  model. 

In  order  to  make  the  connection  between  our  microscopic 
model  and  the  Kittel  macroscopic  model,  we  use  the  defini¬ 
tion  of  the  polarization,  which  is  an  average  of  molecular 
dipole  moments  over  a  given  volume  P=(l/V)Xpi-.  The 
sublattice  polarizations  Px  and  P2  in  the  Kittel  one¬ 
dimensional  model  can  then  be  defined,  in  terms  of  the  mo¬ 
lecular  dipole  moments  in  our  model,  as 


■^1  y  2  Pat »  y  E  Pal *  (2) 

iPa>  0}  iPa<  0) 

In  other  words,  in  the  AFE  state  the  formula  unite  containing 
a  positive  dipole-moment  component  in  the  a  direction  form 
one  sublattice  with  polarization  Pt ,  and  the  formula  units 
containing  a  negative  dipole-moment  component  in  the  a 
direction  form  the  other  sublattice  with  polarization  P2.  The 


TABLE  IE,  Invariants  for  the  antiferroelectrie  transition  driven 
by  a  soft  mode  corresponding  to  the  M3M4  irreducible  representa¬ 
tion. 

^““(Pi^pf)2  Bpi  dp2  dp\  dp2 

^  dx  dx  dy  dy 

MW 

3)=p\P2~PlPl 


difference  between  the  microscopic  theory  presented  here 
and  macroscopic  continuum  theory  of  Kittel  is  that  the  latter 
cannot  provide  the  spatial  relationships  among  dipoles  of 
adjacent  cells,  and  cannot  account  for  the  dipole  tilt  which 
occurs  in  the  ADP  system. 

LANDAU-GINZBURG  FREE  ENERGY  BASED  ON 
MICROSCOPIC  SYMMETRY 

Adding  OP  gradient  terms  to  the  free  energy,  i.e.,  using 
the  Landau-Qinzburg-type  free  energy,  will  allow  us  to  de¬ 
scribe  inhomogeneous  structures,  such  as  orientational  twins 
and  antiphase  walls.  The  invariant  polynomials  of  the  OP 
and  its  derivatives  can  also  be  obtained  using  the  ISOTROPY 
program.  In  this  paper,  we  have  truncated  the  free  energy  at 
the  fourth  power  of  the  OP  to  limit  our  discussions  to  a 
second-order  transition.  We  also  assume  that  the  OP  field 
varies  slowly  in  space  so  that  only  the  first  derivatives  of  the 
OP  are  included  in  the  free  energy.  We  write 

F^AipW+B&W+DA  (»=  1,2,3),  (3) 

where  file  invariant  polynomials  <p{2\  (p\4) ,  and  are  given 
in  Table  EH,  The  coefficients  A  -  A0(T —  Tc),  Bt ,  and  Di  are 
constants. 

For  single  domain  states,  all  derivatives  of  the  OP  must 
vanish;  Le.,  <5f=0.  Considering  the  single  domain  state  S1 , 
the  two  component  OP  is  (pa  ,pb)  with  both  pa  and  pb  non¬ 
zero.  Since  the  symmetry  of  the  low-temperature  phase  is 
known,  the  ionic  displacements  and  dipole  moments  in  each 
unit  cell  can  be  determined  experimentally  from  neutron 
scattering.  In  other  words,  the  ratio  of  the  amplitudes  k 
-pJPb  can  be  measured  in  the  antiferroelectrie  state 

Pb=f<Pa*  (4) 

Because  all  derivatives  of  the  OP  vanished  in  the  single 
domain  state,  energy  minimization  of  Eq.  (3)  leads  to 


=  2Apa+[4Bl(l  +  K1)+4B2 
+£3(3/c-/c3)]p3=0,  (5a) 


dF 

dPl 
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dF 

dp2 


=  2  Axpa+ 

P~{pa,Kpa) 


4B1{\  +  k1)  +  4B2k2 


■+■  B2 


xp\= 0. 


(5b) 


In  order  for  Eqs.  (5a)  and  (5b)  to  give  the  same  solution  for 
pa ,  the  sum  involving  the  B2  and  B3  terms  must  be  equal, 


4B2+B2(3k- k3)  =  4B2k2  +  B3 


(6) 


which  can  be  written  as  a  fourth-degree  polynomial  equation 
in  k: 


K4+i>/c3-6/r2-fc/c+l=0  (7) 

with  b=4B2/B3.  This  equation  has  four  roots, 

(8) 

where 

*1  =  —  T—  7  Vl6+h2+  -^r(l  6+b2)U4ylb+  Vl6+h2, 

4  4  2V2 

(9a) 

b  1  , -  1  - l- 

k2=--+_V  i6+h2+ - (16+h2)1/4V-h+Vl6+h2. 

4  4  2V2 

(9b) 

One  of  the  above  solutions,  Eq.  (9a),  will  be  less  than  1  and 
greater  than  0  (i.e.,  0</c<  1),  which  is  the  solution  we  seek 
since  it  corresponds  to  the  OP  in  domain  Su  with  pa>pb 
(see  Fig.  1).  Because  Eq.  (7)  does  not  contain  any 
temperature-dependent  terms,  the  solution  k  will  be  tempera¬ 
ture  independent. 

Once  k  is  determined,  Eq.  5(a)  allows  the  value  of  pa  to 
be  determined  as  a  function  of  the  expansion  coefficients  in 
Eq.  (3): 


2  — 

Pa=451(l  +  K2)  +  452  +  fi3(3/c-K3)‘  (10) 

We  note  that  the  expansion  coefficient  A  is  temperature  de¬ 
pendent,  A=A0(r-rc),  so  that  the  amplitude  of  the  order 
parameter  satisfies  the  universal  relation  for  a  second  order 
phase  transition,  i.e., 


TABLE  IV.  Equivalence  classes  of  two-domain  (twins  and  an¬ 
tiphase)  structures. 


Domain  pair 

Type 

Relation  to  Sx 

(Si,S2) 

orientation  twin 

{°mI°2  4} 

(Si,s4) 

orientation  twin 

Wbdll  1  1} 

antiphase  structure 

{EWW} 

theory,  the  pairing  set  can  be  much  reduced.  We  find  three 
classes  of  two-domain  structures  to  be  considered,  as  listed 
in  Table  IV  along  with  the  type  of  structure  represented  (the 
trivial  degenerate  structure  that  has  the  same  domain  on  both 
sides  of  the  wall  is  actually  a  single  domain  structure):  two 
orientation  twins,  and  one  antiphase  structure. 

We  will  first  consider  the  orientation  twin  formed  by  do¬ 
mains  Sj  and  S2.  The  domain  wall  orientation  for  such  a 
twin  can  be  determined  by  strain  matching,11,12  noting  that 
this  transition  belongs  to  the  Aizu13  species  42mF222.  The 
strain  allowed  domain  walls  are  the  planes  x=±y  (see  Fig. 
5).  Since  properties  of  the  crystal  will  change  only  along  the 
direction  normal  to  this  wall,  we  will  rotate  to  a  new  coor¬ 
dinate  system  so  that  x 9  is  normal  to  the  wall.  The  OP  com¬ 
ponents  (and  any  other  properties  of  the  twin)  then  will  be 
functions  of  x(  only.  This  means  that  the  derivatives  in  the 
gradient  part  of  the  free  energy,  Eq.  (3),  can  be  rewritten  as 


dp  l  dp  dp  1  dp 

dx  ~v5  dx"  dy  ~^dxf  '  K  J 

For  such  a  twin  structure,  the  OP  varies  across  the  wall 
from  its  value  in  Si  to  its  value  in  S2, 


P  i  = 


Pb> 


x 

xf 


►  +  oo’ 


P  2  = 


Pb>  * 
-Pa,  *' 


►  +  oo 


.  (12) 


In  order  to  match  the  boundary  conditions  in  Eq.  (12),  it  is 
convenient  to  change  the  dependent  variables 


Pi  ~Paf  Pbf 2* 
Pl^Pbfl-pJl- 


(13) 


Here,  the  functions  fx  and  f2  are  normalized  order-parameter 
components  and  have  simpler  boundary  conditions 


Pb^Pa^^Tc-T. 

ORDER-PARAMETER  PROFILES  FOR  TWINS  AND 
ANTIPHASE  STRUCTURES 

The  Landau-Ginzburg  free  energy,  Eq.  (3),  allows  us  to 
describe  inhomogeneous  structures.  We  will  study  both  ori¬ 
entation  twins  and  antiphase  structures  in  the  AFE  phase, 
each  composed  of  two  different  domains  separated  by  a  do¬ 
main  wall.  Since  there  are  four  possible  domains  in  the  AFE 
state,  there  will  be  4X4=  16  possible  pairings  of  domains 
from  which  we  can  construct  orientation  twins  or  antiphase 
walls.  However,  using  the  idea  of  equivalence  in  group 


FIG.  5.  Twin  structure  with  a  domain  wall  oriented  parallel  to 
the  ( 1  TO)  lattice  plane. 
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FIG.  6.  Numerical  solution  for  the  order  parameter  profile  of  an 
orientation  twin  with  7'  =  3.0. 

1,  *'-►-«>  JO, 

0,  >  +  co’  f2  1?  ‘x*l+  +  ocG* 

We  now  apply  the  Euler-Lagrange  equations 


FIG.  7.  Dipole  arrangement  in  the  vicinity  of  an  antiphase  wall 
between  Sx  and  S3. 

For  the  second  type  of  orientation  twin  involving  domains 
Si  and  S4,  the  boundary  conditions  are 


—  00 
+  OD? 


[Pb*  *'-* 

.Pa  * 


—  00 

+  CO* 


(21) 


d  dF  dF 


0-U) 


and  obtain  a  coupled  set  of  ordinary  differential  equations 

d2f\ 

DfaT5=afi+Pf\+yfifl, 


‘7777 -  0/2+  r/l/s-. . 


D=<1  +  k2)(D,+D3), 

a=(l+-*2)A,  ( 

/8=^[4(  1  +  /c2)2^ +4(  1  +  k4)B2+4k(  1  -  k2)B3]  , 


We  make  the  following  change  of  variables 
Pl~Paf\~~Pbf2* 

Pi^Pbfi+PJi*  "  (22) 

Hie  same  set  of  equations,  Eqs.  (20),  for  the  functions /j  and 
ft  can  be  derived  and  similar  solutions  can  be  obtained. 

For  the  antiphase  structures,  we  cannot  use  the  strain 
compatibility  relations  to  predict  the  orientation  of  the  wall 
since  the  strain  compatibility  relations  are  automatically  sat¬ 
isfied.  However,  we  can  simplify  the  equations  to  quasi-one- 
dimensional  using  the  same  procedure  described  above. 
Again,  we  rotate  the  coordinate  system  such  that  x*  is  per¬ 
pendicular  to  the  wall.  The  derivatives  in  the  gradient  part  of 
the  free  energy  can  again  be  expressed  in  terms  of  deriva¬ 
tives  in  the  x'  direction.  By  defining  the  new  normalized 
order  parameter  g. 


?=Pa[4(  1  +  k2)2Bi+24k2B2—  12k(  1  —  *c2}2?3].  Pi=Pag,  P2=Pbg  (23) 

In  order  to  meet  the  boundary  conditions  in  Eq.  (14),  we  and  using  a  rescaling  of  the  space  variable  similar  to  Eq.  (19) 
must  have  (but  dependent  upon  the  specific  orientation  of  the  wall),  we 

can  simplify  the  system  to  a  single  differential  equation  for 
0.  (18)  the  function  g. 

The  xr  axis  can  be  rescaled  to  ^2^ 

IP2  =  -g+g3-  ■  (24) 

For  the  antiphase  structure  formed  by  domains  Si  and  S3 
to  arrive  at  a  simplified  set  of  dimensionless  equations  (see  7)>  the  boundary  conditions  are 


drfi 

^a—fi+A+y’fifl 

d2f2 

=  ~h +fl  +  y'flfi  >  (20) 

where  y*  -  yf [3  is  temperature  independent  Equations  (20) 
have  been  solved  for  other  systems*4  and  the  solutions  for 
the  choice  of  yr  =  3.0  are  illustrated  in  Fig.  6. 
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g  =  tanh| 


(27) 


SUMMARY  AND  CONCLUSIONS 

We  have  proposed  a  Landau-Ginzburg  model  based  on 
microscopic  symmetry  to  describe  a  second-order  antiferro- 
electric  phase  transition.  The  OP  field  used  is  directly  corre¬ 
lated  to  the  local  dipole  moments  of  the  formula  units,  in¬ 
stead  of  the  macroscopically  averaged  polarization.  Our 
model  not  only  can  derive  the  macroscopic  model  proposed 
by  Kittel  but  also  can  address  the  dipole  tilt  and  cell  doubling 
in  the  antiferroelectric  transition.  There  is  no  need  to  assume 
separate  sublattices  in  our  model. 

When  the  antiferroelectric  state  is  switched  to  a  ferroelec¬ 
tric  state  by  an  external  electric  field,  our  model  predicts  a 
unit  cell  size  twice  as  large  as  that  of  the  macroscopic  picture 
obtained  from  previous  investigations.2  This  picture  is  more 
consistent  with  the  microscopic  picture  observed  experimen¬ 


tally  that  the  dipoles  are  actually  tilted  from  the  a  and  b 
directions. 

There  are  four  domain  states  in  the  antiferroelectric  phase 
and  they  form  two  distinct  orientation  twins  and  one  type  of 
antiphase  structure.  The  addition  of  gradient  energy  terms  in 
the  free  energy  allows  us  to  model  multidomain  inhomoge¬ 
neous  structures.  Numerical  solutions  for  the  inhomogeneous 
OP  profiles  describe  the  gradual  change  of  the  dipole  ampli¬ 
tude  and  orientation  across  the  domain  walls.  We  find  that  if, 
in  the  free  energy,  only  the  coefficient  of  the  quadratic  term 
is  assumed  to  be  temperature  dependent,  the  local  dipoles 
have  a  fixed  orientation  independent  of  temperature,  while 
the  amplitude  of  the  dipole  moment  is  a  function  of  tempera¬ 
ture. 
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